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Abstract 

Switch-level simulators model a metal oxide semiconductor (MOS) large scale integrated (LSI) 
circuit as a network of transistor "switches". They can simulate many aspects of MOS circuits which 
cannot be expressed in the Boolean logic gate model, such as bidirectional pass transistors, dynamic 
storage, and charge sharing. Furthermore, the logic network can be extracted directly from the mask 
specification by a relatively straightforward computer program. Unlike analog circuit simulators, 
however, the nodes are assigned discrete states 0, 1, and X (for unknown), and the transistors arc assigned 
discrete states "open", "closed", and "unknown". As a consequence, switch-level simulators operate at 
speeds comparable to logic gate-simulatora. 

In this thesis, a formal model of switch-level networks is developed. The networks in this model 
may contain transistors of different strengths and types, as well as nodes of different sizes and types, and 
hence the logical behavior of a wide variety of ratioed, complementary, and ratioless designs can be 
expressed. In keeping with the concept of a hgfc model, however, both the transistor strengths and the 
node sizes may take on only discrete values, and electrical behavior is modeled in a highly idealized way. 
The operation of a network is characterized by its target stale ft ; etion, which for a particular state of die 
network yields the logic states which the nodes would eventual!) reach if all transistors were held fixed in 
their initial states. This characterization abstracts away the rate at which nodes approach their target 
states and the voltages through which they pass but provides adequate detail for many simulation and 
analysis technique* The target state function can be defined in teum of an abstraction called hgk 
signals, where a logic signal gives a composite description of 'the network at some node much as a 
Thevenin equivalent network gives aeomposfee description of a trocar network at some pott 

Logic signals can be formalized into a simple, discrete algebra with operations describing the effects 
of performing some elementary network transformations. A technique for finding the target state of an 
arbitrary switch-level network can be derived by utilizing concepts from abstract algebra and lattice 
theory. This technique leads to an algorithm for a switch-level simulator which improves on previous 
algorithms in its generality, speed; and fiimplteity. Furthermore, the mathematical formulation provides a 
means for proving useful properties about toe simulajwft »etJao4 and opens up further areas of 
application for the switch-level modeL 
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1. Introduction 

Recently, a new class of logic simulator has emerged specifically for simulating metal oxide 
semiconductor (MOS) large scale integrated (LSI) circuits. These switch-level simulators model an MOS 
design as a set of nodes connected by transistor "switches" with each node assuming a state 0, 1, or X 
(unknown) and each switch a state "open", "closed", or "unknown". Programs such as the author's 
MOSSIM [8, 9] and others [5] show remarkable accuracy and versatility in simulating such logic elements 
as logic gates, pass transistor logic, busses, and both static and dynamic memory. The accuracy results 
because the logic network closely matches the actual circuit, while the versatility results because 
transistors form a common denominator for all LSI design techniques. Furthermore these simulators 
operate at sufficient speeds to test entire LSI systems, because behavior is modeled at a logical rather than 
a detailed electrical level. Unlike previous attempts at developing MOS logic simulators by adding ad hoc 
extensions to gate-level simulators, switch-level simulators are based on a uniform and consistent model 
which provides a powerful level of abstraction for viewing MOS designs. In this thesis the concept of 
switch-level simulation is developed into a mathematical model of MOS logic networks from which 
simulation algorithms and other analytic tools can be derived. 

The ability to implement digital logic has progressed gready in the past decade with circuits of 
increasing size and complexity being fabricated at decreasing cost. Metal oxide semiconductor (MOS) 
technology has played a major role in this "integrated circuit revolution" due to its relative simplicity in 
both design and fabrication. In more recent years MOS design has become part of the university 
curriculum in both electrical engineering and computer science. By using simplified design rules and 
conservative clocking schemes, and by following systematic methodologies such as those presented by 
Mead and Conway [28], the basics of MOS design can be learned in one semester. As this training 
becomes widespread, we will see a new form of integrated circuit revolution in which nonspecialists 
design their own custom integrated circuits rather than relying on the limited variety of commercially 
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available products. 

With the increasing number of custom-designed integrated circuits, and with the growing size and 
complexity of commercial LSI products, the inadequacy of current LSI design techniques has become 
apparent. The semiconductor industry has traditionally relied on humans to design, lay out, and verify 
LSI systems. Typically many man hours are spent, and several prototype chip designs are fabricated in 
developing a single IC design. Industry analysts have extrapolated the current design techniques and 
estimated that a 100,000 device microprocessor (the expected state of the art in 1982) would take 60 
man-years to lay out and another 60 to debug [41]. Rawer than accepting such predictions as inevitable, a 
change in design techniques is called for. 

Computerized tools have been applied to commercial LSI design, but most of these can be viewed 
as extensions of manual techniques <such as graphical .layout systems), or as expert in a specialized 
domain (such as circuit simulators.) , JBoth kinds require close roppetatjonwitha human whauoderstands 
the exact capabilities and limitations of the program. In addition, humans are cequued to bridge the gaps 
between programs with expertise in different domains. For example, before a designiQan be testedby a 
logic simulator, the actual design typically must be translated by hand into a description which can be 
understood by the simulator. This translation process wastes manpower in performing a rather tedious 
task and also decreases the level of confidence provided bjr the simulation. 

Logic simulators form an important class of computerized tools for LSI design. TJieu* utility has 
long been recognized for analyzing designs which by their size and complexity exceed the capability of 
humans to fully understand. The usefulness of a logic simulator dependsgreatly on the consistency and 
accuracy with which it can model the mil range of design techniques available to the designer. Of course 
no logic simulator can model all designs with complete accuracy, because it does not simulate the detailed 
analog behavior. Nonetheless, it should provide as close a model as possible within a set of well-defined 
limitations. As a further requirement, a logic simulator for LSI must be efficient enough to simulate 
entire systems with reasonable speed. The size of single-chip, very large scale integrated (VLSI) systems 
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will far exceed the small scale integrated (SSI) systems for which conyenjipnal logic simulators were 
designed. 

A logic simulator has as its basis an abstract model of hbw digital systems function. This logical 
model describes both the structure and the behavior of a system in terms of a set of primitive elements, a 
set of interconnections, and a set of rates for operation: For a simulator to accurately arid reliably 
simulate a system, the logical model must reflect its actual structure and Operation. 

Unfortunately, the development of togiestmutetorsfiasnbt Kept pace with LSI technology. This 
inadequacy stems hv part from the lack of ftirnidA logic models for describing the behavior of MOS 
circuits. Instead, systems are designed and simulated using an ad hoc combination of Boolean gate 
models, relay models, and electronic models. Such a representation may be appropriate for human 
designers, who can combine different modes of operation and* resotvelhe coftfflcts between these models. 
Computers; however, lack the intuition required to simultaneously view a system at several different 
levels. Computerized tools must be based on models which can describe a [large class^>f systems in a more 
uniform way. 

Most logic simulators; are based on the Boolean gate model which adequately models systems built 
from SSI components but fails to support the wide variety of techniques available to me LSf designer, 
especially for MOS LSI. Many extensions of the Boolean gate model have been attempted with die usual 
result that only a slightly larger class of designs can be simulated and'rrMy 50trfces ! of unpredictable of 
inaccurate behavior are introduced. This claim wifl b^sup^bited by briefty surveying the development of 
logfc simulators with respect to their support ft>rM08LSf design. ' 



1. Brzozowski and Yoeli [10] present a logic model in which "T-etements" provide a-simplified model of 
field-effect transistors. This model, however, only expresses the operation of static logic ^ gates. 
Furthermore, it has not received widespread attention. 
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1.1 The Boolean Gate Model 

The Boolean logic gate model has formed the theoretical basis for logic design ever since die advent 
of electronic logic, to this model a system consists of a set of logic gates connected by unidirectional, 
memoryless wires. The logic gates compute Boolean functions of their input signals and transmit these 
values along the wires to the inputs of other gates. Each gate input has a unique signal source. 
Information is stored only in the feedback paths of sequential circuits. This model directly implements 
Boolean algebra [20] and hence has a well-defined specification which can. guide the simulator 
implementation. 

The Boolean gate model cannot describe many; of the techniques available to the logic designer, 
especially the MOS LSI designer. MOS pass transistor networks can impfcanent cotnbinatiorul logic in 
ways which more closely resemble relay contact networks than logic gate networks (see [2S] or [16] for 
several examples.) Dynamic memory can store information without feedback paths by exploiting the 
capacitances of the wires and the gates of the transistors attached to diem. A variety of bus structures can 
provide multidirectional, multipoint communication. A logic simulator which implements only the 
Boolean gate model provides limited support to the MOS LSI designer. Most existing logic simufcKpa, 
however, extend the Boolean gate model in various ways> Hence, any evaluation of logic gate simulators 
must consider these extensions as wefl. 

Many simulators extend the two-valued logic of Boolean algebra wiw a third value to represent an 
unknown or undefined logic level. This "X" level can mdicate an, uninitiali«d state variable, a signal 
held between the two logic thresholds, or a signal in transition between and 1 or between 1 and 0. The 
X logic level can be handled algebraically by changing the two-valued Boolean algebra to a three-valued 
DeMoigan's algebra [3, 11, 23, 46]. Thus, even with this extension many of the desirable mathematical 



1. A DeMorgan's Algebra satisfies all postulates of Boolean Algebra except for the Law of Excluded 
Middle(A+ -iA= 1). 
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propcrties of the Boolean gate model are preserved. Alternatively, some simulatorsimpjeinent the X level 
by an enumeration technique in which the simulation is repeated with the nodes at the X level set to all 
possible combinations of O's and l's [6, 45]. Nodes which remain at a; unique level iter all combinations 
are set to this level, while all others are set X. Still ether simulators |44] use ad. hoc techniques to 
implement the X level often resulting in inconsistent or a riomalous 'faetevior: The Xiogic level is useful in 
simulating all forms of digital logic including MOS. 

To model the behavior oftros structures, some logic simulators have a fourth or "high-impedance" 
logic level $2]. This H level corresponds to the third state of tri-state logic. To simulate a bus structure, 
the outputs of a number of gates are connected to a common node. Typfcalf y all but one output will be at 
the ft level, and the level of "mat output win* dominate. Unlike the X level which can be viewed as an 
extension of Boolean algebra, the H level violates a basic principle of' me Boolean model, in that a logic 
gate input no longer has a unique signal source. Because the simulator is not based on a well-defined 
mathematical model, it becomes difficult to implement consistently and accurately. The H state may be 
adequate for simulating SSf designs in which only hmited fonra of^trf-state bussescan be implemented. 
The MOS LSI designer, on the other hand, can select from a wide variety of bus designs, such as 
pre-charge/discharge and multiple driver designs. The H state i only partially captures me behavior of 
these hus structures. Nonetheless, me H logic level M seen in simulators Kir both M6S and other forms 
of logic. ■"•■■■■■■•■ -"■-' >■■ 

Some simulators allow a special logic gate to represent the MOS pass transistor [44]. this logic gate 
models a field-effect transistor (FET) as a unidlrtttfonal device with' two inputs and one output as shown 
in Figure 1.1. The simulator cannot hielp in cases m" which ^ &e6i^reciiohal property of the FET r& 
important, such as in circuits where information may flow in eiiher direction, or where the design has a 



1. In actual fact, the bidirectional property of the FET is rarely u|ed intentionally. The ! author has seen 
only a few designs which have signals propagating in bdth directions tijjrQUjglba pass transistor. 
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Fig.M. The FET Logic Gate 
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malfunction due to a sneak path. More recent gate model simulators have implemented bidirectional 
transistor models, but these transistors usually entail a much higher computational cost, and hence their, 
use must be minimized. Furthermore, few of these simulators can simulate arbitrary combinational 
networks of pass transistors, because they cannot model the action of p^JJupresis^ors. Thus th^extensian 
does not fully capture the behavior of the MOSFET. like the h^ : impedance logic level, U also lacks the 
algebraic properties of the Boolean gate model and hence forces mad hoc implementation. 

Finally, some simulators model dynamic memory in a limited fashion [44]. A node is allowed to 
remain at a previous logic level if the outputs of all logic gates connected to the node are at die H level 
This extension is very limited in its generality and its accuracy. 

As new types of MOS logic circuits are developed, designers add more extensions of the Boolean 
gate model to their simulators. These extensions are doomed to failure, because they cannot correct the 
fundamental mismatch between the Boolean gate model and MOS logic circuits, MOS circuits consist of 
bidirectional switching elements connected by bidirectional wires with memory (considering the 
capacitive effects of the transistor gates as contributing to a wire's memory.) Instead, die simulators 
become increasingly cumbersome and unreliable, because they only partially capture the behavior of the 
logic technology. 



1. SIMULOG introduces even greater inaccuracy by failing to differentiate between the undefined (X) 
level and the high-impedance (H) level. 
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Using a gate model simulator requires both intuition about how the design is supposed to function, 
and detailed knowledge of the simulator implementation. The user must explicitly identify the logic 
gates, the Signal directions through pass transistors, the locations of busses, the sites of dynamic memory, 
and sometimes even the feedback paths. Often the actual togic design must be transformed into one 
compatible with the simulator which may not display tie exact same behavior. This transformation 
process not only decreases the level of confidence provided by the simulation, it virtually eliminates the 
possibility of automatically generating the simtilatk)Bo«etwoifk from some specification of the actual 
design. Unless we restrict our attention to a lirrrited class c^ designs, a very sophisticated program would 
be required to analyze the mask patterns for an M(>S layout and tSJnvert this to a gate-level description, 
performing the necessary transformations to provide cctopatibility with me Simulator. Without this 
capability; a togfc simulator cannot be used to help veifry the* obrrecliess bfa layout Gate-level logic 
simulators fit into the MOS LSI design process as shown in Figure 1.2. They provide mainly a 
verification ©f the high level functional description ^plaSliraltea Verification of the actual logfc design! 



Fig. 1.1 Role of Logic Gate Simulators in LSI Design 
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1.2 Analog and Hybrid Simulators 

LSI designers have recognized the limitations of conventional logic simulators for modeling MOS 
circuits and have at times resorted to analog or hybrid simulators. Analog simulators treat me entire 
design as network of analog circuit elements and try to model the detailed waveform at every node over 
time. Simulators such as SPICE [301 and even those which use fester (and nwre approximate) numerical 
techniques such as MOTIS [U\ require very large amounts of computation. Some reports claim the 
amount of computation scales as the square of the network, size [2]. Thus they are practical onry for small 
designs or for small sections of larger designs. Analog simulators, however, are based on a uniform and 
general abstract model and hence have been well received because of thetf consistency and accuracy. 
Furthermore, computer programs exist for deriving me simulation network automaticaUy from the layout 
descriptions (2L 

The amount of computation can be reduced significantly by hybrid techniques such as in SPLICE 
[31] in which some sections of the design are simulated as logic gates and others are simulated as analog 
circuits. Hybrid simulation works well as long as only small, isolated sections of the design need be 
simulated as analog circuits. Unfortunately, a human must decide which portions of the network can be 
modeled as logic gates, and which portions require analog simulation. Furthermore, trying to combine 
analog and logic models in a single program requires rather unsatisfactory approximations at the 
interfaces. For example, if an output of a section modeled as tegk gates is to be interfaced to an input of 
a section modeled as an analog circuit, die program must convert the logic signal into a voltage waveform. 
This, of course, cannot be done with any accuracy, because much of the necessary information is lacking. 
The resultant outputs of the analog section must then be viewed with skepticism. Similarly, if the logic 
simulator were extended to include the X state, it could not be interfaced*) an analog simulator because 
mis state docs not represent a single voltage. Unless great care is exercised, a hybrid simulation could 
well provide die accuracy of a logic simulator at the speed of an analog simulator, rather than vice-versa. 
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For this reason, a human must monitor the simulation very carefully. 
1.3 Switch-Level Simulators 

As an alternative to conventional logic simulators, the author has developed the simulator MOSSIM 
[8, 9] specifically for the logical simulation of MOS LSI. With MOSSIM the Boolean gate concept is 
discarded altogether and replaced with a logical model which closely matches the structure and behavior 
of MOS circuits. A logic network consists of a set of nodes connected by a set of FET "switches". 
MOSSIM uses three logic levels: 0, 1, and X (undefined.) There are three types of nodes: 

1. Input nodes provide a strong, externally generated signal (eg. power lines, 
clock drivers, data inputs, etc.) 

2. Pullup nodes are connected via a pullup resistor to a high voltage. They will 
generatea 1 signaluntess grounded. Theoutpotof an riMOS tdgiegate isan 
example of a pullup node. 

3. Normal nodes cannot generate a signal but can store a signal dynamically. 

Only two types of network elements are allowed: p-type and n-type field-effect transistors. A 
transistor is a three node device which acts as a voltage-controlled switch with no assumed direction of 
signal flow as shown in Figure 1.3. No distinction is made between the labels "source" and "drain". 
When the gate node of a transistor is in die X state, the switch status is unknown: it may be open, closed, 
or somewhere between. The user interface of MOSSIM allows the user to describe the network in terms 



Fig. 1.3. The MOSSIM Transistor Model 
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of transistors, logic gates, and user-defined macros, but these are all translated into a transistor-level 
representation for the simulation. C. Terman has developed a switeh-^i^ simulator patterned after 
MOSSIM [5] but differing in several respects, as is discussed at several points in this thesis. Researchers 
at Caltech [34] also developed a switch-level MOS logic simulator, but not to a degree of accuracy or 
generality required in a serious design tool. Researchers at other laboratories have developed their own 
switch-level simulators based on these earlier designs. 

Switch-level simulators can simulate almost the full range of circuit designs available to the MOS 
designer without any special logic levels or poorly-defined logic elements. Both logic gate and pass 
transistor combinational logic are simulated without difficulty, as aw,both s stafc and dynamic memory. A 
wide variety of bus structures can be simulated including tri-state busses, pre-charged busses, etc. Most 
significantly, theuser need not te» the simulator what type«f^je structure &ji»teodc4only the actual 
physical structure of the design. 

Switch-level simulators have been tested on a wide variety of MOS designs 'ranging from student 
homework problems to a LISP microprocessor chip containing over 10,000 transistors (21]. They have 
proved remarkably general and accurate, correctly simulating logic design techniques which were not 
even anticipated in die simulator design. The confidence in die simulation results is greatly enhanced by 
the fact that the user can see an exact correspondence between the actual design and the simulation 
network. Moreover, the simulation is fast enough that entire designs can be simulated. For the LISP 
microprocessor chfp, MOSSIM requires between 5 and 12 seconds of CPU time on a DEC20/60 to 
simulate each clock cycle. The designers were able to fuDy test the systemoy simulating 700 clock cycks; 
Experience has shown drat simulation inevitably uncovers fatal errors in die design. 

C. Baker [4] has written a program which can Tate layouts specified in the Gdtech Intermediate 
Form {28] and generate the equivalent transistoi'-tevel network. Unlike a program which generates a logic 
gate description, this program needs no special intuition about logic design. It need only look for 
electrical connectivity and transistors in the mask patterns. The simulation network for the LISP 
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microprocessor chip can be generated with 30 minutes of CPU time on a PDP 11/70. This technique has 
proved extremely valuable, uncovering errors both in the layout and the logic design. Furthermore, it 
saves the duplicative effort of entering the design by hand for the two different representations. 
Switch-level simulators can fit into the MOS LSI design cycle at several levels, being applied 
independently to the high level description, the actual design, and the layout, as shown in Figure 1.4. 

1.4 An Abstract Model for MOS LSI 

The generality and accuracy of switch- level simulators suggest that a formal model of MOS based 
on the switch-level concept can be developed. The uniformity and consistency of the switch-level 
approach are precisely those properties which make a concept amenable to a mathematical treatment. 
This model would serve not only as a basis for verifying the correctness of a logic simulator, but also as 
the foundation for new computer tools for MOS design. One need only look at the many advances made 
possible by Shannon's development of logic models based on Boolean algebra [38, 39] to understand the 
value of abstract logic models. While the expected benefits of a MOS logic model are much more 



Fig. 1.4. Role of Switch-Level Simulation in LSI Design 
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modest, its utility could be significant. Unlike traditional switching theory which was developed to help 
humans analyze and synthesize networks containing small numbers of elements, we now want models 
which can form the basis of computer programs to be applied to networks containing thousands of 
elements. This places more emphasis on the generality and uniformity of the model and the algorithmic 
complexity with which it can be implemented. 

In this thesis, a formal switch-level model of MOS logic nctworte.i&dey*lope(L The network model 
closely resembles the network model of MOSSIM but generalizes it in several respects. As with 
MOSSIM, a logic network consists of a set of nodes interconnected by a set of transistors. Unlike 
MOSSIM, however, only two types of nodes: input and normal are allowed. To model ratioed circuits, 
transistors may have different strengths, with a stronger transistor (such as an inverter pulldown) being 
able to override a weaker one (such as a pullup load transistor). A third type of transistor, d-type (for 
"depletion") is introduced which is closed regardless of the gate signal, this new model more closely 
matches the actual structure of MOS networks, because in MOS networks the relative sizes of transistors 
determine the logical behavior. The pullup node used in MOSSIM is a rather ad hoc way of representing 
this. The new model can also describe a wider variety of networks, mcfuding circuits which rely on 
multiple levels of ratioing. 

In MOSSIM, each normal node is modeled as having a capacitance of unknown value which can 
store a signal dynamically but dannot a^ve^'mis signal onto another 'WSk'% "'a different state. 
Unfortunately this model cannot describe the behavior of many bus designs in; which a relatively high 
capacitance bus node is connected to a lower capacitance node (such as the sttfage-dode of a 3-transistor 
dynamic RAM cell) resulting in both nodes obtaining the same logic state # was originally on the bus. 
Our new switch-level model can model this effect by assigning each normal node a size, where the signal 
on a larger node w#! predominate when connected to a smalleruode. 
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Our abstract model describes both the time and electrical behavior of a network in a highly 
idealized way. The time behavior is described by the target slate function giving the logic states which die 
normal nodes would reach for a particular set of input node, transistor, and initial normal node states. 
For designs containing no critical races, the logical behavior can be modeled by repeated application of 
the target state function. To model the electrical behavior, the target state is denned in terms of the set of 
steady state voltages in an "order of magnitude" electrical network. This class of networks models the 
conducting transistors by linear resistors, where the conductances of the resistors for different strength 
transistors differ by orders of magnitude. As a consequence, any path to an input node containing only 
transistors with strength greater than or equal to some value is modeled as overriding any path containing 
a transistor with strength less than this value. Similarly, the normal nodes are modeled by capacitors 
where the capacitances of the capacitors for different size nodes differ by orders of magnitude. As a 
result, the target states formed on a set of nodes throygh charge sharing depends only on the state(s) of 
the largest node(s) in the set : Furthejrmore, no attempt is made to accurately compute, the node voltages. 
Instead, they are classified into the three logic states 0, 1, and X. This model provides a simplified view of 
ratioed circuits and charge sharing which adequately describes die logical behavior of most MOS circuits. 

Although the target state is defined in terms of an electrical model, we will find that the target state 
of an arbitrary switch-level network can be computed without evaluating any, electrical networks. Instead, 
by introducing an abstraction called logic signals, an iterative method for computing the target state can 
be developed which uses only operations in a simple, discrete algebra. A logic signal provides a 
composite description of a switch-level network at some node Jbr a particular set of node and transistor 
states, much as aThevenin network [15] provides a composite description of a linear network at some port 
for a particular set of network parameters. However, whereas finding the Thevenui equivalent generally 
requires solving a set of simultaneous linear equations, finding the lqgic signal requires much less effort 
A simple set of rules describes the logic signals created by the input and normal nodes in their initial 
states and the effects on a node of other nodes connected ftrough conducting transistors. With- these 
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rules we can develop an equation expressing a set of constraints which must be satisfied by the signal for 
each normal node in terms of the initial node signal and the signals for input nodes and for other normal 
nodes connected through conducting transistors. It can then be shown mat die minimum solution of this 
set of equations equals the set of logic signals describing me network at each node, and the set of steady 
logic states can easily be derived from mis solution. Logic signals can be formalized into simple, discrete 
algebra with a domain corresponding to die signal values and operations describing the effects of the rules 
for logic signals. This algebra allows us to apply elementary concepts from abstract algebra and lattice 
theory to develop techniques for computing the target state. These techniques can be further developed 
into efficient simulation algorithms. 

US Relation to Relay Networks 

The switch-level MOS model can be viewed as an extension of Shannon's relay network model p8, 
395. The algebra of signal strengths introduced in Chapter 5 bean many similarities to Boolean algebra. 
As Shannon observed, a relay can be viewed as a switch with conductance 1 when closed and when 
open. 1 The rules for connecting relays in series and in parallel and the methods of analyzing relay 
networks are special cases of those for transistor networks. 

MOS networks, however, have several characteristics which are not found in relay networks. Fust, 
relay networks are used as a current-driven logic, in which the logic state of a node is determined by the 
connection between the node and the current source. Thus a simple characterization of the connection to 
the signal source determines die state Of a node. MOS networks, in contrast, are used as a voltage-driven 
logic, where die state of a node is determined by both its connection to die supply voltage and its 
connection to ground. One must characterize bom the states of the signal sources and the connections to 
them to determine the state of a node. Furthermore, in a voltage-driven logic erroneous behavior may 



1. Shannon actually described die state of a relay by its hindrance, die complement of its conductance. 
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result due to a short circuit between signal sources, and hence we also require the state X for logical 
completeness. Second, relay networks do not allow ratioing in which one closed switch can override 
another. Thus, a Boolean characterization of a conductance path suffices. Third, relay networks can only 
store information in feedback paths. No modeling of dynamic memory or charge sharing is needed. 
Finally, most theoretical work on relay networks was, conducted before the widespread availability of 
digital computers. Thus, most techniques were developed to aid die hand design of small circuits. The 
standards by which ideas are measured change greatly when they are to, be incorporated into 
computerized tools to aid the design of very large systems. 

1.6 Outline of Thesis 

In the next chapter, the details of how the switch-level model describes die structure and operation 
of MOS netwprks is presented. By mcnldir^ the network structure at a transistor level and by allowing 
transistors of different strengths and nodes of different sizes, this model covers a large variety, of MOS 
design techniques in a way which closely matches the actoaji circuit, d.e^jns., The switch-level, model can 
be viewed as either a simplification of analog network models or an extension of relay network models. 
The time behavior of a network is described by the target state function, giving the logic states toward 
which the nodes are driven or charged given the current node and^ansisior states. The value of this 
function is defined in terms of the steady state voltages, in ,$ .ltycar^ete&r^sd. network that models die 
transistor network. The logical behavior of many MOS networks is described by repeated applications of 
their target state functions. In computing the target state* the transistors are modeled with time-invariant 
elements, thereby simplifying the analysis considerably. 

In Chapters 3 and 4 the electrical circuit-oriented view of Ih^jarget state function provided by the 
definition given in Chapter 2 is transformed into a more abstract and logical view. It is shown that the 
target state can be defined in terms of the steady states of a set of logical conductance networks, where a 
logical conductance network represents a switch-level network in which each transistor is either 
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nonconducting or fully conducting. The concept of logic signals is then developed to express the 
behavior of logical conductance networks. With the logic signal abstraction we can derive an equation 
which gives the steady state of a logical conductance network, and consequently die target state of a 
switch-level network, which does not require evaluating any electrical networks. 

In Chapter 5, an algebra of logic signals is developed With operations describing die effects of a set 
of network tramformations. This algebra allows us to apply elementary concepts from abstract algebra 
and lattice theory to the study of switch-level networks. 

In Chapter 6 the mathematical formalism presented in Chapter 5 fe used to derive a technique for 
computing the target state of a switch-level network. This development utilizes, only the logic signal 
abstraction as expressed by the algebra of Chapter 5 and two equations which are derived in Chapters 3 
and 4 from the analysis of the electrical model. Although we could arrive at the desired results more 
directly by utilizing additional properties of the ekctrkal makf, this approach demonstrates the power of 
our abstract approach. 

In Chapter 7, the abstract solution technique of Chapter 6 is developed into an efficient algorithm 
for a swhch-tevel simulator. By exploiting the sparseness of the network the Smiuiator requires at most 
linear time to simulate one clock cycle for almost all networks. This algorithm improves on previous 
swhch-Ievel simulation ahjorithms in several respects. Some performance data for MOSSIM is presented 
to demonstrate die performance characteristics of switch-level simulation and how it compares to logic 
gate simulation. 

In Chapter S, the simplified timing model of MOSSIM is investigated more closely to see for what 
classes of systems it is valid. Possible methods of implementing logic simulators with other timing models 
are presented. In addition, a tern/ay simulation algorithm is developed which uses die X state to detect 
potential races in MOS networks. This algorimm is a straightforward extension of Brzozowski and Yoeffs 
algorithm for logic gate networks [HJ. Ternary simulation requires a much more accurate and efficient 
implementation of the X state than is required for functional simulation, because the X state will become 
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the most prevalent state in the network. The algorithm presented in Chapter 7 provides this accuracy and 
efficiency. 

Finally, in Chapter 9 some ideas for further improvements of logic simulators and for future 
applications of the switch-level model are described. 

1.7 Notation 

In the remainder of this presentation, the following notauonal conventions are observed. Scalar 
values are denoted with lower case letters (e,g. a, b); vectors with boldface, lower case letters (e.g. a, b); 
and matrices with boldface, upper case letters (e.g. A, B). Mathematical domains, i.e. sets of values with 
particular mathematical properties are denoted with script, upper case letters (e.g, JL, ft), while ordinary 
sets are written with italic, upper case letters (e.g. A, B). The extension of a domain 9 to vectors of size n 
is denoted 9 n , and the extension to matrices with n rowsandm columns is denoted 9 n - m . 
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2. The Switch-Level Network Model 

11 Introduction 

In this chapter a model of the logical structure and operation of MOS networks wilt he presented. 
This model attempts to capture those aspects of an MOS circuit which affect its logical behavior, while 
ignoring many of the detailed electrical properties. This network model extends the network model of 
MOSSIM but in a way which provides a consistent level of abstraction. Like MOSSIM the network can 
be extracted direcdy from a specification of the layout. The time behavior of the network is also 
described in a simplified way by the target state function. Given the current network state, this function 
yields the state toward which the nodes move without considering the rate at which these changes occur. 
This function bears a strong resemblance to the excitation function used in fogfc gate and relay models. It 
is assumed that the reader has a background in MOS logic design comparable to that provided by Mead 
and Conway [28]. 

12 Network Structure 

A logic network contains a set of input nodes / = { ij ^ }, a set of normal nodes 

# = {«],...,«„ }, and a set of transistors T = { fj /^ }. 

Input nodes provide strong signals from sources external to the network, much like voltage sources 
in electrical networks. Examples of input nodes include the supply (VDD) and ground (GND) 
connections as well as signals supplied through input pads. Normal nodes have states determined by the 
operation of the network. Each normal node n^ has a size capj, where capj is an element in the set 
K = { Kp . . . , k }. A normal node can store charge to provide dynamic memory. The size of a node 
gives an approximate characterization of the amount of charge it can store, where sizes are ordered 



«1<«2 < "- < V 
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When normal nodes are connected together they will share charge and settle in a state dependent only on 
the state(s) of the largest node(s). The values in K have no properties other than their ordering; they only 
indirectly represent actual physical capacitances. This model provides a simplified view of charge sharing 
which is valid for most actual circuit designs. The number of node sizes q depends on the kinds of MOS 
networks to be modeled. For most networks, q = 1 will suffice. For those networks which rely on a 
sharing of charge between a high capacitance node and a low capacitance node for their logical behavior, 
q must equal 2 or more. For example, Figure 2.1 shows a three-transistor dynamic RAM circuit which 
relies on the high capacitance of the bus (size k-j) to override the charge on the storage node of the cell 
during a write and on the drain node of the storage transistor during a read. 



Fig. 2.1. Ratioless MOS Design 
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A transistor is a three terminal device as shown bdow. 

t drain 
gate —|[~ 

• source 

No distinction is made between the source and drain connections. Associated with each transistor 4 is a 
strength str^ where stTj is in the set T = { yj, . . . , y }. The strength of a transistor gives an approximate 
characterization of its conductance when turned-on, with strength values ordered 

0<y 1 Cy i <...< 1f;r 

The strength values in T have no properties other than their ordering; they only indirectly represent 
actual conductances. The number of allowable strengths p depends on die kinds of networks to be 
modeled. For networks which do not rely on ratioed resistances such as CMOS designs, all transistors are 
of equal strength and p = 1. Most ratioed nMOS or pMOS designs can be modeled with p = 2, where 
pullup and pulldown loads have strength y^ and all other transistors have strength y 2 - Some designs, 
including certain static RAM cells rely on multiple levels of ratioing and hence require a model with p 
equal to 3 or more. A transistor can be either n-type, p-type, or d-type. AD act as voltage-controlled 
switches as follows: 

n-type p-type d-type 

gate signal eJBxt yatesiynal gfigsfit* eate signal effect 

open "dasit closed 

1 closed l open l closed 
X unknown X unknown X closed 

A d-type (for "depiction") transistor can serve as either a toad resistor for a depletion mode nMOS logic 

gate, or a polysilicon-difrusion crossover such as is seen in some designs. When a transistor is in a 

"closed" state it provides a conductance between the source and drain nodes with value characterized by 

the transistor strength. When a transistor is in an "unknown" state it provides a conductance of unknown 
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value between (inclusively) the conductance when "open" (i.e. 0.0) and when "closed". This model 
provides a simplified view of ratioed circuits in which a connection through a stronger transistor will 
always override a connection through a weaker, as will be defined more rigorously later in this chapter. 

Examples of a variety of MOS circuits are shown in Figure 12, The first three show common 
nMOS and CMOS logic gates. The fourth one shows a foroeable inverter which, when connected in a 
ring with another inverter, can statically store 1 bit 



Fig. 2.2. Examples of ^fQS Networks 
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Transistors with strength less than y (called weak tr an s is tors) are used in ratioed circuits, where a 
stronger transistor may override a weaker one. Generally these weak transistors are configured only in 
limited ways, such as pullup loads on logic gate outputs. At nates we wul want to exploit the 
characteristics of networks restricted in the use of weak transistors to simpHfy the mathematical 
development or to improve the efficiency of an algorithm. In particular, we define a restricted network 
follows: 



In a restricted network every transistor of strength less man y has either its 
source or drain connected to an input node. 



All of the circuits in Figure 11 are restricted networks. An example of an unrestricted network is shown 
in Figure 13. In this example the pass transistor is assigned a strength y| to* titi&cM o^twhen a sneak 
path forms through the "kill" aiid pasr transistors with the inverterinput e^juaf to 0^ me4nverter output 
will go to X, while die output of the pass transistor will go to 0. Almost afiiactual MOS designs can be 

--•■•-.• r , 

represented as restricted networks, because weak transistors are used almost inclusively as pullup (for 

nMOS) or pulldown (for pMOS) loads^ with one side connected ^eitijer V3DD loe GND. Unrestricted 

networks seem plausible, however, so, we wfll develop results for this mcwgefterri case. We shall find 

f 
mat although unrestricted networks require a more complex mathematical development, dieir study will 



bcitii"- hhd .( ;f5omv;-'ii;n.-;.' 



Fig. 13. Unrestricted Network Example 
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provide much insight into networks containing transistors with gate nodes in die X state. 

2.3 Network Representation 

The structure of a network is represented by the sets /, W.andT.the vectors cap and str (indicating 
the node sizes and transistor strengths), and the following functions: 



TTYPE: r-»{n,p,d} the transistor type 

GATE: T-*IUN the gate node 

SOURCE: T-+IUN the source node 

DRAIN: T-+IUN thedrainnode 



2.4 Logic States 

Each normal node flj has a iogie state yjG {O.l.X }• In Watts of me actual voltage Vjatnode^: 

y. = ■ «fr - OJ9<^:^ T 

where V~ and V + are die logic thresholds. Note that our logic model makes no attempt to accurately 
model these logic thresholds, and they are used here only to aid our Mfoimal discussion. Each input node 
jj has a logic state Xj € { 0,1, X } widi the same interpretation. The values U and 1 correspond to the 
Boolean logie levels. The value X indicates either an unknown or Hmfe/w«/ logic level. An unknown logic 
level arises when an ambiguity in the network condition prevents a unique determination of a node's logic 
level, such as from uninitialized state variables. For example, when power is applied to a bistable device 
such as a Nor gate latch, the output will obtain a valid, but unknown logic state. An unknown level 
corresponds to a voltage either below V - or above V . An undefined logic level arises when the network 
operation creates a voltage which could lie between the two logic thresholds, such as due to a short circuit 
or improper charge sharing. Hence an undefined level corresponds to a voltage which may be anywhere 
between 0.0 and V^. 
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These two concepts differ slightly in that an "unknown" value satisfies the Law of Excluded 
Middle, while an "undefined" value does not For example, if y represents die state of an uninitialized 
bistable device, then y + -iy = 1. On the other hand, if y represents the state of a node along a shorting 
path from VDD to GND, then y + -iy = X. Some circuit designs exploit the Law of Excluded Middle 
during the initial power-up sequence to assure that all feedback paths will be initialized to valid logic 
levels. No known algorithm, however, can utilize information about "unknown" logic values in a 
completely general way, except by enumerating over all possible combinations of Boolean values. Thus, 
to avoid an exponential algorithm, we shall not attempt to distinguish between the "unknown" and 
"undefined" logic levels but instead use the single value X. 

Each transistor /j also has a logic state Zj € { 0, 1, X }, where indicates "open", 1 indicates 
"closed", awl X indicates "unknown". Although transistor states and node states «re different physical 
phenomena, we will use the same mathematical objects to represent both. 

15 Network State 

At any instant in time the state of a network is given by the the logic states of the input nodes 
xE{e,l,X> m , the states of the nonnal nodes j 6 { 0, 1, X } n , and the states of the transistors 
z€{0, 1, X} k . Under stable conditions, the transistor states x are factions of the node states. Suppose, 
for example, that node nj is the gate node for transistor /j (ie. a^ = GAT^t)). Then the value of z* is 
given by the following table 

riTPE(t) 

Yj B p d 

1 1 

11 1 

XX X 1 

A similar table would result if (he gate node were an input node. The function trans (x, y) denotes the 

transistor state z resulting from the node states x and y. During the actual circuit operation, some t"»e 
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may elapse between a change in node state and the resulting change in transistor state. Hence the total 
state of the network must include the transistor states as well as the node states. 

2.6 The Target State Function 

With the network in state (x, y, z) each normal node n { will move toward a target state y { given by 
the function 

y 4 = target j(x, y, z). 

The target state of a node equals the state the node would eventually reach if the input nodes were held 
fixed in state x, the transistors were held fixed in state z, and the normal nodes were initialized to state y. 
As shall be seen, the target state y is the steady state solution of a time invariant network with parameters 
x and z and initial conditions y. We can view target as a vector-valued function 

y = target (xy,z). (2.1) 

giving the target states for all normal nodes. The state y may never actually be reached, however, 
because the transistor states will change in response to the changing node states, and the input node states 
may be changed externally. The function target only describes a tendency in the network and not a 
definite reality. However, it provides a basic characterization of the logical behavior of a switch-level 
network. Much of the development of this thesis will be directed toward a mathematical formulation of 
the target state function. 

Define the function slep x as 

step x (y) = target (x, y, trans (x,y)). (2.2) 

For a particular state of the input nodes step x gives the target states of the normal nodes as a function of 
their initial states, assuming the transistor states are functions of the initial node states. 
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During actual aeration, the network may not move through me succession of states predicted by 
the function step x due to changing transistor states and changing input conditions. However, for a large 
class of networks, the ultimate behavior is equivalent to the behavior modeled by successive applications 
of the function siep x as long as the input nodes remain unchanged. That is, if the normal nodes initially 
have state y and the input nodes are held fixed in state x, the network will eventually reach a state 
phase(x, y) defined as 

Phasdx^) = ^ step*(j) (2.3) 

where the superscript k indicates k applications of the Junction step^. Furthermore, the networks of 
interest will be guaranteed to stabilize after a bounded number of «tep& Thus for some k<oo, 
siep^Hs) = siep x k+ %) = pkautx, y). Once the network arrives at this state, it will remain there until 
some input node is changed. This ignores the possibility that nodes may lose their charge due to leakage. 
With current technology, in which clock speeds are measured in megahertz while leakage times are 
measured in milliseconds, this assumption is appropriate. Examples of systems which can be modeled by 
repeated applications of step % include combinational logic, any system free of critical races, and a variety 
of systems which can be modeled with unit delay logic elements. The limitations of this assumption wul 
be discussed in Chapter 8. 

For a large class of MOS systems an equation for the function target will lead to a description of 
the complete functional behavior of a network. In other words, the behavior of a system can be modeled 
by repeatedly freezing the states of the nodes and transistors, computing the target state for each normal 
node, and then updating the node and transistor states accordingly. "~ This technique has obvious 
advantages over modeling the detailed voltage waveforms on each node. 
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Concepts similar to the target state function have been applied to the study of both relay networks 
and logic gate networks. The target state function describes the excitation of a switch-level network, i.e. 
the node states created in response to the current states. Similarly, the excitation of a relay network is 
defined as the states which would appear on the relay coils if aM relay contacts ace held fixed in their 
current states, while title excitation of a logic gate networkis defined as the outputs of me logic gates as 
functions of their current inputs. In general, the excitation of any logic network can be defined as the 
logic states which would form on toe aodefc if all active elements (e>g. transistors,; relays* m logic gates) 
were held fixed |n their current states, ' 

Huffman [22] first recognized toe importance 0f Ae network excitation as providing a basic 
characterization <)f the dynamic behavior of a logic network, although he expressed H in the focm of a 
flow table rather than a function. Wito Ihis (^aracteruatieiUiB«jci»pfjm# physical behavior is abstracted 
away, such as toe rate at which the nodes move toward through their exoitations and the analog values 
through which they pass. While such a characterization caniKA O^rtcct certam error conditions dependent 
on the detailed voltages or toning, U provides a useful level of abstraction. Although logic states are 
formed in switch-level networks in much different ways than in logic gate or relay networks, these three 
kinds of networks share much in common when viewed as systems computing logical functions. 

2.7 Specification of the Target State 

For both relay and logic gate networks, the excitation function for a node can be defined with a set 
of Boolean functions in a relatively straightforward way. With switch-level networks, on the other hand, a 
more complex formulation is required, because the logic elements are bidirectional, and because toe state 
depends on the relative conductances of toe pulhip and pulldown paths acting on it or on the relative 
capacitances of nodes with winch it may share charge. We will define toe target state in terms of a linear 
electrical model called toe "order of magnitude" model With this model toe concepts behind toe 
switch-level model can be expressed in relatively conventional mathematical terms. In later chapters, a 
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rather unconventional algebra is presented to express these concepts more directly and to allow the 
development of efficient simulation algorithms. The order of magnitude electrical model serves only as a 
means by which these more abstract concepts can be lriotivated&d derived. 

Order of magnitude networks contain voltage sources; resteers, and capacitors where the resistor 
conductances and capacitor c^r^acitaiKes are spedfedaipewefs of a rttio parameter p. Transistorsm toe 
X state form conductances of unknown value bounded by powers of p, and nodes M me X state form 
unknown voltages; Henceywemusteonsder toesetofposs^^ 

order of magnitude network when the parameters and initial condffiom range over sefe of values. The 
target state of a node in a switch-level network is defined in terms of Ae set of possibk steady state 
voltages on me corresponding node hi the order of magnitude network asp is made very large; such oat 
the conductances of difleneM strength transistors in me 1 state and die^eapacitancels ^different sbse 
nodes differ by orders of magnitude; Thfe provides a simplified view of ratioed logic and charge sharing 
in which only me dominathig effects are considered. 

The order of magnitude network corresponding tea swfcch-tevel network in a particular state 
(vy,*) isconstructed with dements shown in Figure 2A Each input node t is modeled by a voltage 
source with negative terminal connected to GND and with voltage x.-, where X:= 0.0 if x. = 0, 
x j = V dd tf *j = *• and x j ran 8 es over me xt of voltages { v 1 0.0<v<V dd } if x, = X. A normal node 
/ij of size k^ is modeled by a capacitor with one side connected to GND and with capacitance ap* where 
*is an arbitrary positive constant This capacitor is initially charged tea voltage y= defined in terms of 
toe logic state yj in thesameway Xj is defined in tcrmsofX:. A trasststorwkh strength y k and state 1 is 
modeled by a resistor of conductance mp* where a is ata artwtrary positive constant A transistor wife 
strength y^ and state X is modeled by a resistor with conductance ranging over the set 
{ 9 |00<3<<»P }. where a is the same constant used when the transistor state is 1. The values of a 
can be different for different resistors and capacitors. It wHI be shown later that fee values of feese 
coefficients do not affect the value ofthc target state. 
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Fig. 2.4. The Order of Magnitude Network Model 
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For a particular value of p, the conductance parameters in an order of magnitude network can be 
described by two matrices G G * n x n and E G * n X m , where % denotes the set of real numbers. Each 
element gj: of G equals the sum of all conductances between nodes n- Y and n-, while each element ey of E 
equals the sum of all conductances between node «j and the voltage source corresponding to input node 
L When the switch-level network contains transistors in the X state, these matrices may range over 
infinite sets: 



G min (p) < G < Gr*(p) 
E min (p) < E < E^Xp) 



with the restriction that G be symmetric. Note that the partial order < is defined between matrices in the 
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usual way, i.e. A < B if and only if a^ < by for all i and j. The elements of G min , G m ", E min , and E raax are 
polynomial functions of p, and hence the elements of the conductance matrices G and E are bounded by 
polynomial functions of p, although they may take on arbitrary values within these ranges. With this 
formulation, we are assuming that transistors with the same gate node behave independently when that 
node is in the X state. This models the possibility that transistors may have slightly different threshold 
voltages, and hence when the gate node has a voltage close to one of the thresholds, the transistors may 
behave quite differently. Furthermore, as shall be shown in Chapter 3, this assumption allows us to look 
only at the minimum and maximum conductance values for each transistor when computing the target 
state. 

The remaining parameters of an order of magnitude network are described by a vector c(p) G 3> n , 
giving the capacitances corresponding to the normal nodes, and the vector x 6 { v 1 0.0<v<V dd } m 
giving the settings of the voltage sources corresponding to the input nodes. When the switch-level 
network contains input nodes in the X state, this vector may range over an infinite set 
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< x < x"", 



where if xj = X, x mi \ = 0.0 and x™^ = V dd . The initial conditions of the order of magnitude network 
are described by the vector y G { v 1 0.0<v<V dd } n giving the initial voltages on the capacitors 
corresponding to the normal nodes. When the switch-level network contains normal nodes initially in 
state X, this vector may range over an infinite set 



-rata 



< y < y™". 



Let Vj(G, E, c(p), x, y) denote the steady state voltage on node n y for the network with parameters 
G, E, c(p), and x, and initial conditions y. Since the network contains only passive, linear elements, this 
voltage must be unique. Furthermore, since the network contains no floating capacitors, all node voltages 
must lie between 0.0 and V dd> When nodes or transistors in the X state are present in a switch-level 
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network, this voltage can range over a set Kj(p) where 

Kj(p) = { Vj(G, E, c(p), x, y) 1 6 = 6 T , 6 mta (p)^6^G m », E^^E^E" 1 "^), 

x mte <x<x m ", f^<y<f^}. 



This set is uniquely determined by the structure and state of die switeh-fevel network, the constant 
coefficients a (which will he seen to be unimportantX and the ratio parameter p. 

The target state of a node «j?is defined in terms of the set Vfp) as me ratio parameter p is made 
very large. As this occurs, ar# conductance paths to input nodes formed by transistors with state 1 and 
strength greater than onequalto y* will dominate over paths containing transistors of, strength less than 

i" ' : ■■■>%■! 

y^. Similarly, the charge on capacitances formed by normal nodes of size kj, will dominate over the 
charge on capacitances formed by normal nodes of lesser size. To form a proper logic state (i.e. or 1), 
there can be no conflict between the pullup and pulldown paths or between the high and low charges. As 
p is made very large V^p) should approach eimer-the%et|M }loMh4set{ V^ }, If we take the limit as 
p approaches infinity, the set V(p) should converge to one of these twjo sets. Thus, if we define the set 



K,°°as 



'" -,*•»»» 



then the target state on node «j is defined as 



{ 



uY^MttOT (14) 

y i = i 1. ^°°«{ 

X, else. 



In this formulation, the X state jmaavBtm'^IKm &e set Kj(p) converges to some value hot equal to 0.0 or 
V44, indicating erroneous behavior due to a short circuit or improper charge sharing, or when the set 
K|(p) fails to converge to a single value, indicating an ambiguity in me steady state voltage caused by 
nodes or transistors in the X state. This definition of the target state in terms of a limiting process 
expresses in a mathematical way the concept that the switch-level model considers only the dominating 
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effects acting on each node, either through conductance paths formed by transistors in the 1 state to input 
nodes or by charge sharing between normal nodes. When the dominating effects conflict or when nodes 
or transistors in the X state create uncertainties, the target state equals X. 

For example, Figure 2.5 shows a switch-level model of an nMOS Nor gate and the corresponding 
order of magnitude network. Let us see how the target states would be defined for several sets of inputs. 
If k»j = 1 and mj = X, feen flj = ejp, g^ = a-jp , and Oj0^g*<o*p'', where a j, a^ and «i are 
arbitrary positive constants. Tim gives a set of steady state voltages ftraodcitj 

1 I I «1 + <«2 + «3*» ~ ~ «1 + "7* J* 



.00 



and therefore V y = { 0.0 } and y 4 = 0. Now suppose that taj = and uv^ = X. Then g A = a^p, 
9 2 = OJ), and 0.0<g3<o 3 ^ 2 . This gives a set 



w = C f¥ |J^4L<.».<!|lSfcl;l 



.00 



and therefore Vf = { v [&.0<¥<V dd }aiid yj = ». 



Fig. 25. Example of aa Order of Magnitude Network 
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2.8 Relation to Actual Circuits 

We have defined the target state in terms of an electrical network model as the ratio between the 
conductances of different strength transistors and between the capacitances of different size nodes is 
made very large. It then follows that the switch-level model would correctly describe an MOS circuit 
containing only transistors in which the length-width ratios differ by orders of magnitude along with 
nodes with capacitances that also differ by orders of magnitude. 

In designing actual MOS circuits, of course, one does not use transistors with length-width ratios 
which differ by orders of magnitude. Instead, the relative sizes of transistors along a conducting path 
from VDD to GND are set so that the node voltages will lie sufficiently within the logic thresholds. 
Furthermore, transistors are sized according to their required speed, power, and driving capabilities. As a 
result, the pullup load in a clock driver may have much greater conductance than the pullup in an 
ordinary inverter, even though both perform the same logical function ~ to provide a 1 signal in the 
absence of a stronger signal. This transistor sizing can be viewed as an optimization of our order of 
magnitude model to improve the electrical charact^|tics^ii^^reta3«^:^,sa^ logical behavior, lac 
almost all MOS circuits the logic value formed on a node depends only on the dominating effects. 

Similarly, node capacitances span a wide range of values, but the logical behavior is affected by 
node sizes only in a few isolated locations, such as in pre-charged bus circuits. Typically the large 
capacitance node (e.g. the bus) greatly exceeds me, smaller capacitance node, and hence our order of 
magnitude model more nearly approximates the actual circuit in this case. 

To model the logical behavior of a correctly designed MOS circuit we need only characterize 
transistor and node sizes according to their logical function in the networjt. This correctness can be tested 
prior to the simulation by a computer program which compares \ the conductance and capacitance 
parameters of the circuit against the proposed logic network. Thus the simulation model can assume that 
the circuit is correctly designed in this respect and take a more abstract view <of ratioed circuits and charge 
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sharing. 

Not all MOS digital circuits can be modeled by a switch-level network. For example, the 
switch-level model cannot describe circuits in which slight variations in voltages can represent different 
logic values, such as one-transistor dynamic RAM designs using sense amplifiers to detect these 
variations. Furthermore, our model can only describe the behavior of ratioed circuits or charge sharing 
when there is a clear precedence between the different transistor conductances and between the different 
node capacitances. Other forms of MOS circuits can only be modeled with partial accuracy. For 
example, the switch-level model ignores the effects of floating capacitors in "bootstrapping" node 
voltages above W^ or below 0.0. This technique, however, is used primarily to overcome the saturation 
effects of the field-effect transistor, but our model ignores these effects anyhow. Thus a circuit which 
utilizes bootstrapping for this purpose can be simulated with a switch-level model, but the simulation will 
not check whether the bootstrapping actually occurs. Except for these limitations, switch- level networks 
provide an accurate and simple way to describe the behavior of MOS logic circuits. 

2.9 Comparison to Other Switch-Level Models 

The logic networks allowed by the simulator MOSSIM [8] can be described in the model using only 
one node size (q = 1) and two transistor strengths (p = 2). Input and normal nodes in MOSSIM 
networks are modeled as input and normal nodes, while pullup nodes are modeled by one of the 
following circuits: 



T 
n 






The new model corresponds more closely to the actual circuit implementation, because transistors of 
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different sizes define the behavior of ratioed circuits. The pnllup node of MOSSIM provides a rather 
inelegant model of this and also lacks generality. All transistors in a MOSSIM network are modeled by 
transistors of strength yj- A MOSSIM network always corresponds to a restricted network in the new 
model. 

The logic networks allowed by C. Term^n's switch-level simulator {5} are for the most part identical 
to MOSSIM networks. A more general model of charge sharing is allowed, however, in which each 
normal node has a real-valued capacitance. This technique is applied when a set of normal nodes is 
interconnected only by transistors in the 1 state, none are connected to input nodes, and none are in the X 
state. In such a case, the program sums the capacitances of those nodes in state 1 as well as the 
capacitances of nodes in state 0. If one sum exceeds the other by at least a factor ofi, the nodes are all set 
to the corresponding state, and otherwise they are set to X. No attempt is made to perform this 
calculation when transistors or nodes m the X state are present Instead, the aodesare all set to X. 

At first this approach might seem superior to our model of charge sharing in which nodes are 
assigned a size in the set K and the value on a larger node always-overrides the value on a smaller. Using 
real-valued capacitances more closely matches the actual electrical behavior and utilizes only information 
easily calculated from the layout specification. However, it seems as tf transistors MMhe X state cannot be 
dealt with in a consistent, way with this model. Tennaa has chosen not ft> simulate the effects of 
transistors in the X state with great accuracy. Instead, nodes are sometimes set to X even when they would 
have state or 1 regardless of the conductances of transistors in the X state. For example, if a high 
capacitance node in state is connected to a low capacitance node in state 1 by a transistor in state X, 
both nodes are set to X even though the high capacitance node would remain in state even if it shared 
charge with the other node. 
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Suppose that we were to utilize real-valued capacitances in our logic model but tried to provide a 
more accurate model of transistors in the X state. That is, a node would be set to or 1 if it has this 
unique state regardless of the conductances of transistors in state X, and otherwise it would be set to X. 
Consider, for example, the network shown in Figure 2.7 containing a set of nodes of increasing 
capacitance connected by transistors ha the X state. Suppose initially that node n l is set to 1 and all others 
are set to 0. With mis model, die target states of nodes n^ and nj equal X, because these nodes would 
have undefined states if only they share charge. Nodes »j and n^, on the other hand, would have target 
states 0, because no setting of the transistor conductances could cause them to be charged above the logic 
threshold. We have defined the target state as the steady state solution of die network, and hence when 
the nodes are set to their target states, the network should remain stable until a transistor or input node 
changes state, if we set nodes n± and n^ to X, however, the target stale of node /13 will become X, because 
by our naive approach, we must consider the case in which the X states on nodes nj and n^ actually 
represent high voltages, and tfiese nodes share charge with just ^3. Hence, die original target state is not a 
true steady state solution. Similarly, if we set node /13 toils new target state, me target state of node ^ 
becomes X, indicating that the previous target state was also not stable. In me final steady state solution, 
the initial charge on node n± has created an undefined state on a node with 3© times greater eapacitance. 
Thus, this model of charge sharing yields unstable solutions and also lack* accuracy. 



Fig.2£. Charge Sharing Anomaly 



XXX 

In . 2 n . 3 r u 4 

1.0 2J> 8.0 30.0 





1 


"2 


"3 


"4 


initial 


1 











stepl 


X 


X 








step2 


X 


X 


X 





step3 


X 


X 


X 


X 



43- 



More sophisticated approaches could be devised such as storing the range of possible voltages for 
each node, but for any scheme there seems always to be a circuit which would be modeled incorrectly. 
The problem occurs because we are trying to combine exact electrical concepts such as real-valued 
capacitance with abstract logical concepts such as the X state for both nodes and transistors. With the 
logic model mapping many possible conditions of the electrical network into the state X, it cannot model 
behavior which depends on detailed electrical properties with sufficient accuracy or consistency. By 
adopting a more absolute view of charge sharing in which a larger node can always override a smaller, we 
obtain a more uniform level of abstraction leading to a more consistent modeling. Of course to describe a 
design in terms of our logical model the process translating the electrical design into the logic model must 
decide how the design should be viewed logically. For the network shown in Figure 2.7 this would 
require some rather unsatisfying decisions. In fact this network really should be modeled as an analog 
circuit, because its behavior depends too much on exact electrical properties. 

This aspect of the logic model design indicates a fundamental trade-off with abstractness and 
consistency on one hand and a desire to combine concepts from several different models on the other. 
Since we are concerned at the moment with developing the mathematical aspects of the switch-level 
model, the more consistent and abstract approach will be chosen. For other applications, a different 
choice may be appropriate. 

2.10 Derivation of the Switch-Level Network 

Switch-level simulators have proved quite successful in simulating networks extracted by a 
computer program directly from a description of the mask layouts. The combination of network 
extraction and simulation provides an important check of a design. 
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For networks which can be expressed in the MOSSIM network model, this layout extraction is 
relatively straightforward. The program need only follow the electrical connectivity in the network and 
find the transistors and their types. A simple set of rules allow one to translate this node and transistor 
description into a MOSSIM network. For example, in depletion toad nMOS technology, any depletion 
mode transistor with YDD as the drain node can be assumed to be a pullup load, while all other 
transistors are generally strong transistors. If the extraction program also calculates the lengths and 
widths of the transistors, it can verify that the ratioed circuit wOl operate correctly by determining 
whether the highest resistance pulldown path can override the lowest resistance pullup path for nodes 
which have independent pullup and pulldown paths. If these worst case conditions are not satisfied, a 
warning message can be issued. By performing these static checks of the circuits, we avoid the need to 
check the pullup and pulldown ratios dynamically as the simulation proceeds. 

As we generalize the switch-level model to include multiple levels of rauoing, charge sharing, and 
arbitrary use of weak transistors, the layout extraction becomes more difficult The extraction program 
can calculate transistor resistances and node capacitances without great difficulty, but no general rule can 
take these parameters and assign transistor strengths and node sizes. For example, recognizing that the 
forceable inverter shown in Figure 12 requires three different transistor strengths would require a much 
more sophisticated algorithm than our rule for finding the MOSSIM network. Similarly, identifying 
which nodes win be sources of signals during charge sharing and which wifl be recipients cannot always ■ 
be done with complete accuracy. 

If we tailor the layout extraction program toward a particular class of designs and not try to utilize 
the full generality anowed by the model presented here, the extraction can be made reasonably 
straightforward and reliable. For example, circuits using more man two transistor strengths are relatively 
rare and even then appear only in limited configurations. The extraction program can look just for these 
configurations and for other portions of the design apply simple rules such as those used in deriving the 
MOSSIM network. Similarly, a bus node can usually be identified by ks large capacitance relative to the 
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nodes with which it may share charge. Such nodes can be assigned size «2 anc * a -^ otners s ' ze K l *° r most ' 
cases. Thus, layout extraction should still work well for this more general switch-level model. 

One can see the advantage of Terman's model from the standpoint of layout extraction, where the 
relation between the physical capacitances and the logical behavior is computed dynamically as the need 
arises. Similarly, a program could compute the relative conductances of the pullup and pulldown paths 
dynamically to model ratioed circuits in a very direct way, although this is more difficult than computing 
the relative capacitances in charge sharing. For some applications, one may be willing to sacrifice the 
accuracy and consistency with which the X state is modeled to gain this direct correspondence between 
the electrical parameters and the simulation model. 

2.11 Summary 

The switch-level model provides three major simplifications i over more detailed analog circuit 
models: ^ 

1. Timing is not modeled in great detail. Instead, the dynamic behavior of a 
network is modeled by. 9 se^uence^of tajgejustates, where each target state 
represents the steady state solution of a ume-invanant network. 

2. Node states are characterized by just the three logic levels 0, 1 arid X, where 
states and lajr^ duxing^i^ 

from an ambiguity in the network : or Atot enoneous operation. 

3. The effects of ratioed logic and charge sharing ; are modeled in a simplified way 
as if the conductances formed by transistors of different strength and the 
capacitances of nodes of different size differ byorf^gtfpi^tyclf. iThis, 
assumes the circuit correctly obeys the ratio rules ancl hence the exact voltages 
need not be computed during the simulation. 

These assumptions lead to a uniform and consistent view of MOS logic designs in which only those 

aspects which determine die logical behavior during normal operation are considered. As has been seen, 

these assumptions lie in a very delicate balance. If we try to introduce greater accuracy in one area, such 

as incorporating real-valued node capacitances, we can lose consistency in another. 
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3. Logical Conductance Networks 

3.1 Introduction 

As was discussed in Chapter 2, die target state function provides a basic characterization of the 
dynamic behavior of a switch-level network. The value of th& fiiriction was defined in terms of the set of 
possible steady state voltages for each node in a linear electrical network in the limit as the ratio 
parameter p approaches infinity. This definition helps clarify the relation between die switch-level model 
and MOS logic circuits but provides litde aid in devising efficient simulation algorithms. In mis chapter 
we start a transition from a circuit-oriented view of the switch-level model to a more abstract and logical 
view. A new form of network is introduced called logical conductance networks to focus attention on a 
key aspect of computing the target state. A logical conductance network represents a switch-level 
network with each transistor either nonconducting or ftuly conducting. Trie nock states have values 0, 1, 
and X, but unlike switch-level networks, the network elements are "logical" conductances which can take 
on values from a small discrete set The target state of a switch-level network can be defined in terms of 
the steady states of a set of logical conductance networks. The steady' state of a logical conductance 
network is in turn defined in terms of the limiting case steady state voltages in an order of magnitude 
network with a unique set of parameters and initial condition^ Lc^catco^hrctar^ networks provide an 
intermediate level of abstraction between electrical networks and switch-level networks, 

12 Properties ofthcHectricalModd 

We have defined the target state in terms of a class of linear time-invariant electrical networks 
containing voltage sources, capacitors with one side connected to ground,, and linear resistors. The 
settings of the voltage sources are given by the vector x, the capacitances are given by the vector c, and 
the conductances of the resistors are given by the matrices G and E. The initial conditions of the network 
arc defined by the vector y, giving the initial voltages on ate capacitors. The nodes in this network 
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correspond to both the normal nodes and input nodes in the switch-level network. Some properties of 
these networks will be given here. Formal derivations of these properties are given in books on electrical 
network theory such as Desoer and Kuh [15J. They rely only on KirchofFs' Current arid Voltage Laws 
and on Ohm's Law. 

In a linear network containing only passive, Bnear, time-iirvariant resistors and capacitors, any node 
connected by some conducting path to a voltage source (which has its negative terminal connected to 
GND) will have a steady state voltage uniquely determined by the voltage source settings and the resistor 
conductances. Such nodes are said to be driven. When node n- x is driven, its steady state voltage is a linear 

function of the voltage source settings x: 

m 
v 4 = 2 ayxj. 

The coefficients ay depend only on the resistor conductances and obey the following properties: 

ay > 0.0, foraflj 
m 

2 ay = 1.0. 
j *= 1 

These properties follow from the fact that for any possible Set of voltage source settings given by the 
vector x, vj is bounded below and above by the minimum and maximum elements of x, respectively. 
That is, if ay were negative for some value of j, then an out-of-range voUage would be obtained by setting 
Xj to a positive value and all other voltage sources to 0.0. Similarly, if the coefficients did not sum to 1.0, 
then an out-of-range voltage would be obtained by setting all voltage sources to the same nonzero value. 

Nodes which have no conducting paths to voltage sources are said to be charged. Since every node 
has a nonzero capacitance and the network contains no floating capacitors, the steady state voltage of a 
charged node will be uniquely determined by the node capacitances, the resistor conductances (only 
whether each conductance is zero or nonzero), and the initial node voltages. When Hj is charged, its 
steady state voltage is a linear function of the initial node voltages y: 
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vj = 2 b ijyj . 



j = l 
The coefficients by depend only on the resistor conductance and node capacitances and obey the 

following properties: 

by > 0.0, for all j 
n 
2 by = 1A 

j = l 

These properties follow from the fact that for any possible set of initial node voltages y, the steady state 
voltage on any charged node is bounded below and above by the minimum and maximum elements of y, 
respectively. 

A node is either charged or driven to its steady state voltage, and therefore if we adopt the 
convention mat ay = 0.0 for all j when n y is charged and by = 0.0 for all j when \ is driven, the two 

modes can be described by a single equation: 

ffl n 

V; = 2 ayXj + 2 byyj. (3.1) 

j=l j=l 

The coefficients obey the following properties: 

ay > 0.0, for alii and j 
by > 0.0, for alii and j 



m 



2 ay + 2 by = 1.0, for alii 



m 



2 ay • 2 by = 0.0, fbraHL 

We will also be interested in networks where me elements off, 6, and c are given by continuous 
functions of the parameter p. In this case the coefficients ay and by will be given by continuous 
functions of p. For any positive value of p, the network described by the conductance matrices <%p) and 
Up) and the capacitance vector c(p) will be a passive, linear network with no floating capacitors, and 
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therefore the properties of the coefficients listed above must hold for all p > 0.0. Therefore, if we define 



ay and by as the values 






then 

m n 

p-,oo V i " .V' i t i * * yj 

and the coefficients ay and by obey the following properties: 

ay 00 > 0.0, for alii and j 

by 00 > 0.0, for alii and j 

m n 

2 ay 00 + 2 by 00 = 1.0, for alii 



j=l j F l 

m n 

2 ay 00 • 2 by 00 = a0» fpralli. 
j= 1 j=l 



These properties show that the limits used in the definition of th£ target state are mathematically 
well-defined. 

3.3 Simplification of the Target State Definition 

When a switch-level network contains transistors or nodes in the X state, the target state of a node is 
defined in terms of the set of possible steady state voltages for the node in an electrical network as the 
network parameters and initial conditions are varied over uncountably infinite sets. Fortunately, we need 
not evaluate all of these possibilities, because we only wish to know whether each set of voltages 
converges either to the set { 0.0 } or trie set { V^ } as p approaches infinity. If it can be determined that 
such a convergence does not occur, the target state is X. Thus we need not find the enure set of steady 
state voltages for all possible network parameters and initial' conditions for each node, but only certain 
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propcrties about this set 

33.1 Node Voltages 

First, let us consider the effects of having the voltage sources range over die settings x'^x^x"" 
and the initial node voltages range over y^y-Cy"". Let die vectors x' and y' be any vectors satisfying 
die following requirements: 

**<*'<*-. and x-'-^x—j => x-^^x'^x— j (3.2) 

y*Sr*£ir. «d y"Vy-*i =» y -, i <y' i </-r <3J) 

Hie following theorem shows that we need only evaluate die network for this single set of voltages. 



ill. 

For any conductance matrices 6 and E and raparih >nr«» vector c, if 

V x = {v f (6,£«,x,y)|x , *<x<x-, y"*<y<y-} 

then 

*i = {V^} if and only if v^fi, E, c, x', y*) = V^ 

V x = {0.0} if and only if Wj(€, £ ]c, x', y*) = QJ> 



for any x' and y' satisfying equations 32 and 3 J. 

Proof of Theorem 3.1: 

We have already seen that the voltage Vj is given by a linear function of die dements of x and y as 
shown in equation 3.1, including in the limit as p approaches infinity. Furthermore, all coefficients are 
nonnegative and sum to IjO. Clearly, vj equals V^ if and only if Xj = V^j for all j such that a- > 0.0, 
and yj = V^ for all j such that by > 0.0. Therefore, since x'j equals V^ if and only if 
x - ^ = x""j = V dd , and similarly for y j, Vj(G, E, c, x', y') equals V^ if and only if Vj(6, E, c, x, y) 
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equals V dd for all x such that x min <x<x max and y such that y min <y<y ma \ The proof for V { = { 0.0 } 
follows identically.! 

This theorem shows that even though the voltage formed by an input or normal node in the X state 
can range over an entire set of values { v | 0.0<v<V.jj }, we need only evaluate networks with this node 
having a unique voltage v such that 0.0 < v < Vjj. This single value can be used to test the sensitivity of 
nodes in the network to this X state, where sensitivity is indicated by the values of the coefficients a^ and 
bjj. Thus we can simplify the definition of the target state from one in terms of entire sets of voltage 
parameters and initial conditions to one in terms of a single set x' and y'. Furthermore, the exact values 
of the voltages are unimportant, only whether they equal 0.0, Vjj, or lie between (exclusively) 0.0 and 

V dd" 

3.3.2 Conductance Matrices 

When the switch-level network contains transistors in the X state, the target state is defined in terms 
of the steady state node voltages as the conductance matrices range over infinite sets, G min <G<G m " and 
E rain <E<E m ". Furthermore, while the elements of G min , G" 1 ", E rain , and E™* are functions of p, this 
property does not hold for all matrices in the set, which makes it difficult to establish the limiting case 
voltages. Observe, however, that rather than finding the entire set V^p) for each node, if we could 
determine the minimum and maximum elements of the set, we would have sufficient information to find 
the target state. The following lemma demonstrates an important property of electrical networks as the 
resistors vary between their minimum and maximum values. 
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Lemma 3.1. 

Suppose a passive, linear, time-invariant network contains a single variable resistor with conductance h. 
If vj(h) denotes the steady state voltage on node «j as a function of h, then for any h, h""", and h"" such 
that0J)< h"*<h<h''", either 

VjCh - ") < Vj(h) < Vjtfi!"*) 



or 



VjOO < v,<h) < VjCh-). 



Proof of Lemma 3.1: 

The proof must consider twocases. ; 

First, suppose that node n- is charged when h = OJO, le. there is no conducting path from n^ to a 
voltage source. Then for nonzero h, n^ could be connected to more charged nodes than when h = 0.0, or 
a conductance path could exist from a voltage source to nj and hence ru is driven.. Thus, .there can be a 
discontinuity in vj as a function of h at h = 0.0. For nonzero values of h, the steady state current through 
mis resistor must equal 0.0, because this resistor is not contained in any loop in the network. As a result, 
the value of h can have no effect on the steady state voltages, and Vj must remain constant for any h > 0.0. 
Therefore, either h = h** = O.Oand v^ft"") == Vj(fi),orh >0.0andv i (h) ■'=•*$?*). 

The more difficult case occurs when node /ij is dfivefi for all values Of h, Le. there is always a 
conducting path to a voltage source, m Appendix I an equation is derived for the node voltage V: as a 
function of h by an analysis of mufti-port networks. This equation has the form 

Vj(h) = Vi (0.0) + j^^ (3.4) 

where b is nonnegative. This equation indicates that vj is approximately linear with respect to h for small 
h but approaches a constant asymptote as h becomes very large. Taking the derivative with respect to h 
gives 

dh (U + b h) 2 " 
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which equals 0.0 only in the trivial case where a = 0.0. Therefore Vj has no minimum or maximum 
between h = h™ 1 " and h = n™ unless it is a constant function, and one of the two inequalities must hold. 
I 

With this lemma we can show mat the minimum and maximum steady state voltages for each node 
can be determined by considering only the minimum and maximum conductance -values for each 
transistor. 



Theorem 3.2. 

Let E and G be the following sets of conductance matrices: 

E = {E|ej k = e^j k or»| k = *' , ffj k } 

G = { 6 1 g jk = g- A or fljk = g^ and g jk = g kj }. 



If 



and 



V i = { Vi (G, E, c, x, y) | E^E^E"", 6 B,ta <6<6 m ", 6 = G T }, 



v{ = { vjCe, e, t>mtm € E, 6 € G }, 



then 



minimum Kj == mimmumy^, (3JS) 

maximum V- x = maximum Vf. (3.6) 
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Proof of Theorem 32: 

Suppose some pair of conductance matrices G and E give a minimum value for Vj(G, E, c, x, y). 
For any element of G such that g" ta j k < gj k < g""j k , Lemma 3.1 shows mat we could set g: k to one of the 
values g^jj, or g™"^ without affecting Vj, or else vj would not have been a minimum value eriginaliy. 
This process can be repeated for all such ^, until we have a matrix 6'€<J such that 
Vj(G, E, c, x, y) = Vj(6', E, c, x, y). A similar process would find a matrix E'€f such that 
vj(G, E, c, x, y) = Vj(G\ E', c, x, y), and therefore equation 3.5 must hold. A similar technique proves 
equation 3.6. 1 

This theorem shows that although transistors in the X stale may create arbitrary conductances 
within some range, we can find the minimum and maximum node voltages by considering each transistor 
to be either nonconducting or fully conducting In general, however, we must enumerate over all 
combinations of these two possibilities for all transistors in the X state, because the steady state voltage on 
a node can be minimized or maximized with some transistors nonconducting and others fully conducting 
Furthermore, the voltages on different nodes can be minimized or maximized under different conditions. 
Thus, unlike the node voltages, we must still evaluate, the network for a number of sets of conductance 
parameters, but mis evaluation is now finite. Later we will show that the target state can be computed by 
a more direct method. Note also that this proof assumes that the variable conductance parameters are 
independent of one another, which stems from our original assumption that transistors with the same gate 
node win behave independently when that node is in the X state. 
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3.3.3 Revised Definition of the Target State 

The results of Theorems 3.1 and 3.2 can be summarized by giving a new, but equivalent definition 
of the target state. Let x' and y' be any vectors satisfying equations 3.2 and 3.3 and define the sets £(p) 
and G(p) as 

Kp) = {E|e jk = e mi " jk (p)ore jk = e ,n » jk (p)} (3.7) 

Gip) = { G | g Jk = g rain jk (p) or g jk = g m " jk (p), and g jk = g kj }. (3.8) 

For any node n- v let Fj'(p) denote the set 

V{(p) = {v i (G,E,c(p),x'.y')|Ee£(p). G€G(p)}. 

The matrices in the sets £(p) and G(p) have elements which are continuous functions of p. Therefore, we 
can take the limit of this set of possible steady state voltages as p approaches infinity to get the set 

*T°° = '*L *Y( P ). 

1 p— » 00 ' 



The target state y j is then given as: 



{ 



0, V™ = {0S} (35) 



'1 " 1 I.>' 1 «- t V (id > 

X, else. 



3.4 Rational Functions 

In the formulation of the order of magnitude network model given in Chapter 2, a fully conducting 
transistor of strength -y k is modeled by a linear resistor of conductance ap*. Similarly a normal node of 
size « k is modeled by a linear capacitor of capacitance ap . In both cases, a is an arbitrary positive 
constant and can be different for different capacitors and resistors. Let us generalize this definition to 
model a fully conducting transistor of strength y k with a resistor of conductance g(p) where g is an 
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arbitrary rational Junction of degree k such that g(p) > OJ) for all p > 0.0. That is, g can be expressed by 
an equation of the form 

where N and D are both polynomial functions of p, and if N has degree n and D has degree d, then 
k = n - d. Similarly, we will model a normal node of size k^ by a capacitor of capacitance c(p) where c 
is an arbitrary rational function of degree k such that c(p) > 0.0 for all p > 0.0. 

Observe that this generalization has no effect on our definition of the target state, because any 
rational function a of degree k can be expressed in the form 

a(p) = ap k + HpX 

where a is a constant and b(p) is a rational function of degree less than or equal to k- 1. Therefore, for 
any e> 0.0, there exists a constant pq such that 

(«-*)p* < a(p) < (o+e)p k , for aD p>pQ. 

A rational function of degree k behaves like the function ap as p becomes large. This generalization 
wftl assist our mathematical development, because the dfmain of rational functions has many of the 
properties of the domain of real numbers, i.e. they both (okm fields [27]. We can use expressions such as 
"a+b" to denote "the function of p which gives the value a(p)+b(p)" and assume that the + operation 
in this expression satisfies the usual properties of addition. We can also replace an interconnection of 
resistors by a single equivalent resistor in the order of magnitude model much as one can replace an 
interconnection of resistors having real-valued conductances by an equivalent resistor in ordinary 
electrical networks. 
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The degree of a rational function a can alse be defined as the maximum value k such that 

p _oo p k 7 " 
This function generalizes the usual idea of the degree of a polynomial function. Observe that 
<feg(0.0) = - oo. All other functions usedliefe have integer-valued degrees. For any rational function a, 
we will use the notation a to denote the vawe 

a 00 = ^aO)). J 
p — » oo 

This notation will only be used when deg(a) < 0, and hence the limit is well-defined. 

Let us define the domain «F as the set consisting of the constant function 0.0 along with all rational 
functions a such that a(p) > 0.0 for all p > 0.0. We will deal mainly with rational functions in this domain. 

Ifa.bG 5 then 

defa + b) = maxidegHt), detfb)) (3.10) 

<feg(ab) = <feg(a) + <fe«(b) (3.11) 

<feg(L0/a) = -<feg(a) (3.12) 

defa - b) < nmidggi^ &#>)) (3.13) 

Note that when functions are subtracted, only a weak statement can be made about the degree of die 
resuldng function, and furthennore 5 is not closed uncter this proration. 

As a final notation regarding Rational (unctions, define thwe equivalence relation ~ on. rational 
functions as 

a ~ b if and only if <feg(a) = <feg(b). 

The relations < and > are defined as 

a < b if and only if deg(a) < deg(b) 
a > b if and only if <feg(a) > deg(b). 
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These relations are extended to vectors and matrices in Ae usual way, eg. a ~ b if and only if a^ ~ bj 
for all i. These relations characterize die limit of an expression which is often encountered in equations 
for node voltages in radoed circuits. For any a,bGf: 



{ 



, v _ 1,0, a > b (3.14) 

lim afp} I ^ ^ 

p _ooa(p)+b(p) P *«" 

r ■ a, a ~ d, 



where a is a constant such that 0.0 < o < 1.0. 

15 Equivalent Networks 

We are only interested in the steady state behavior of order of magnitude networks as p approaches 
infinity, and even men we need only a partial characterization of the node voltages. Thus many of die 
details of the electrical network can be ignored. This idea of ignoring certain details can be expressed in 
mathematical terms by defining equivalence relations such that networks which differ only in 
unimportant respects are equivalent 

The equivalence relation ca is defined on elements of die set { v | O-O^v^V^ } as v ~ v' if 
v = v' or if 0.0 < v < V^ and 0.0 < v' < V^. This relation defines three equivalence classes: { 0.0 }, 
{ V^ }, and { v ) 0.0 < v < V^ }, Which correspond cfosely to the logic states 0, 1, and X. This relation is 
extended to vectors in the usual way, Le. V ~V if v 4 ~ v'j for all L The ib^wmg lemma expresses the 
fact diat the exact voltages on die nodes in the order of magnitude model are unimportant, only Whether 
they equal 0.0, V^ , or lie between (exclusively) 0.0 and V^. 
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Lemma 3.2. 

For any conductance matrices 6 and E and capacitance vector c, if x =^ x' and y c^ y', then 

v(G,E,c,x,y) ~ v(G,E,c,x',y'). (3.15) 



The proof of this lemma closely follows the proof of Theorem 3.1. Let v denote the steady state node 
voltages for the vectors x and y, and v' denote the steady state node voltages for the vectors x' and y'. 
Referring to equation 3.1, Vj equals V dd if and only if x= = V dd for all j such that ay > 0.0 and y: = V dd 
for all j such that by > 0.0. Therefore, since x': (or y '=) equals V dd if and only if X: (or y:) equals V dd , v'j 
equals V dd if and only if Vj equals V dd . Similarly, v'j equals 0.0 if and only if Vj equals 0.0, and by a 
process of elimination 0.0 < \/\ < V dd if and only if 0.0 < V: < V dd . I 

We can also show that if the coefficients of the conductance and capacitance parameters in an order 
of magnitude network are varied, the steady state voltages will be equivalent under ~ as p approaches 
infinity. 

Lemma 3.3. 

Suppose an order of magnitude network contains a single variable resistor with conductance ijp . If 

vj(p, tj) denotes the steady state voltage on node «j as a function of p and tj, then for any constants 

ij 1; il2>0.0 

Proof of Lemma 3.3: 

We have already seen that for a particular (positive) value of p, if n-. is charged when tj = 0.0, its 
steady state voltage will be the same for any positive value of tj, and therefore Vj(p, tj^) = Vj(p, tj^) 
including in the limit as p approaches infinity. If «j is driven for tj = 0.0, then equation 3.4 can be 
applied to order of magnitude networks to give 

a(p) • Tjp k 

V:(p, Tj) = V:( P , 0.0) + T . (3.16) 

1.0 + b(p) • Tjp* 
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In the derivation of this equation in Appendix I it is also shown that b(p) > 0.0 for all p and that 
deg(&) < detfo) < 0. Define Vj°°(i») as 

*i*(n) = *** vjCp,,) 

* p— ♦ 00 ' 

and consider the following three 
L&x(b)<-k. 



1 a— »0O ' 



which is independent of if. 



v i°° (, »> = ^L^*^ ^fS 



p — * a» * a 

which is also independent of if. 

l*f(b)= -k. 
where 



p— »00 
fill I 
p-OO 



* = ^NpV*. 



In this case the voltage depends on the value of n. The derivative of this function with respect to if is 
given by: 



*■ (1.O+01J 2 * 

which equals 0.0 only if a = 0.0, m which case Vj 00 ^) is a constant function. We know that 
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O.O^Vj 00 ^))^^ for aH values of ij, and therefore if v^Ojj) - 0.0 (or V^) for some value ^^ then 

-vu = oo - 

in which case vj (ij) is a constant function. 

Combining these three cases, we can see that Vj (nj must either be a constant function or else 
0.0<v i oo (^)<V J y for all ij>Q.Q. Therefore, the possible values oJf v^Cfl) ibr pesJtiye q must be 
equivalent under =^.l .. -.-».-. v 

Lemma 3.4. . i 

Suppose an order of magnitude network contains a single variable capacitor with capacitance y\p . If 

Vj(p, n) denotes me steady state voltage on node ^ as a function of p and n, then for any constants 

q 1 ,ij 2 >0.0 

p _► 00 * x p — ♦ 00 * A 



■ f . » -- — ■■' \ ■■ ' ■! i-P—i n i.. - .— — ^ — 1 1 1 m i iiWf w ti » «■■» i»fi ■ +A**A*1imwt*m**m^ ^ i y..»^ifciiw**ifci. t .. l |i. i >h ii I 'M i ] - 1 - M i ri i> ii 



Proof of Lemma 3.4: 



Suppose the variable capacitor is associated with node ik. This parameter wiQ affect only me steady 
state voltages of nodes connected by some conducting path to n-. and then only if fter?: is, ho conducting 
pam from /i: to a voltage source. For any such node Up * 

, - c(p)vi<p,<he)^Vyi _ 

Vj(p, ij) = L — r » , (3.17) 

where c equals the sum of all capacitances connected by some path to a when ij— O.Q. By considering 
three cases: de&c) < k, de&c) > k, and de&c\ ^ k, we can, show ?thal Vj 00 ^) must f$»er be a constant 
function or else 0.0 < Vj°?{n) < V^for all positive values of ^muc^ as, in the proofof Lemma 3.3. I . ,• ; 
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These two lemmas can be extended to fte case in wakfa the conductance of the variable resistor or 
the capacitance of the variable capacitor is given by an arbitrary rational function of degree k. To see this, 
suppose the terms np k in equations 3.16 and 3.17 are replaced by terms ap k + d(p), where d is a rational 
function of degree less than or equal to k-1. Then as we take me Hmft as p approaches infinity, all 
instances of d will drop out, leavmg the oiigiiud equations. 

These two lemmas Can be combined into a single result by defining the equivalence relation s 
between two order of magnitude networks N and N" as N » N' if and only if 

« ~ € 

i ~ r 

where are uiiprimedsymbob refer to the nctwoitN and u^ 



i3L3L 

Suppose order of magnitude electrical networks N and N* have steady state node voltages v(p) and v*(pX 
respectively. ThenifNsiN\ 

p — ♦ 00 p— » 00 



The proof of this theorem follows directly from Lemmas 3.2, 33, and 3.4 (when extended to arbitrary 
rational functions) and the transitivity of me three equivalence relations, this theorem shows that 
aitiwugh network equivalence is defined in terms of trw structore of the network ^t also raphes a Ram of 
behavioral equivalence. It abo justifies our earner s tatem e nts mat the coefficients of me order of 
magnitude network elements can be arbitrary positive constants without affecting the value of die target 
state. 
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3.6 Logical Conductance 

As an aid in studying switch-level networks in which all transistors are either nonconducting or folly 
conducting, let us define a new type of network called logical conductance networks. Like a switch-level 
network, a logical conductance network contains input nodes and normal nodes, where each normal node 

«j has a size cap^ € K = { icj k }. Instead of containing transistors, however, these new networks 

contain logical conductances, which may have values in the set { } U T = { 0, y± y_ }., The 

correspondence between a switch-level network with all transistors either nonconducting or fully 
conducting and a logical conductance network is quite simple: a nonconducting transistor forms a logical 
conductance (i.e. an open circuit), while a fully conducting transistor forms a logical conductance equal 
to its strength. Examples of the logical conductance networks corresponding to the switch-level model of 
an nMOS inverter are shown in Figure 3.1. A switch-level network with no transistors in the X state has a 
unique corresponding logical conductance network. When a switch-level network contains transistors in 
the X state, however, we must consider the logical conductance networks formed for all possible 
combinations of these transistors being either nonconducting or fully conducting. By focusing our 
attention on logical conductance networks, we temporarily set aside issues concerning transistors in the X 
state and solve a simpler problem. Once methods for analyzing logical conductance networks have been 



Fig. 3.1. Examples of Logical Conductance Networks 

x = 1 

Switch-Level Network -p 




* 



x = 

I 



T 



* 
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developed, we will be able to generalize these methods to handle arbitrary switch-level networks. 

The steady state behavior of a logical conductance network is defined in terms of the steady state 
voltages in an order of magnitude electrical network with a single set of network parameters and initial 
conditions. That is, a logical conductance of strength y k is modeled by a resistor of conductance g(p), 
where g is an arbitrary rational function of degree k in the set 5. A normal node of size k^ is modeled by 
a capacitor of capacitance c(p) where c obeys the same requirements as g. Input node L in state Xj is 
modeled by a voltage source with voltage xj, where Xj = O.Oifx: = 0, Xj = V dd ifxj = 1, and X: equals 
some arbitrary voltage such that 0.0 < Xj < V^ if Xj = X. When normal node /t has an initial state y-„ the 
corresponding capacitor in the order of magnitude network is initially charged to a voltage y: defined 
simSariy. 

The steady slate of a logical conductance network is defined as me set of logic states (denoted with 
the vector y ) which the normal nodes would eventually reach when started in an initial state y. That is, if 
the network is modeled by an order of magnitude network with parameters 6(p), €(p), c(p), and x and 
initial conditions y, and Vj 00 is defined as 

v i°° = ^ v i(«P).E(p).c(p),x,y) 



the. 

f 0, v t °°=(LQ (3.19) 

V x , ao<v i 00 <v dd . 

The target state of a switch-level network can now be defined in terms of logical conductance 



networks rather than order of magnitude networks. From the matrices G and E in the order of magnitude 
model of a logical conductance network we can define two matrices G and E describing the logical 
conductances connecting the nodes in the network. That is, if degiq^d = t, then gy = y k , and if 
gy = 0.0 then gy = 0, and similarly for E. Let the sets { E} and { G } equal the sets of logical 
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conductance matrices corresponding to the sets of order of magnitude conductance matrices E and G 
given by equations 3.7 and 3.8. These sets describe the set of logical conductance networks for all 
possible combinations of nonconducting and fiilly conducting transistors. Let y (G, E) denote the steady 
state of the logical conductance network with logical conductance matrices G and E, with the remaining 
parameters cap and x, and initial conditions y assumed impheidy. The target state can then be defined in 
terms of these steady states as 



{ 



1, y ^G, E) = 1 for all G G { G } and E G { E } (3.20) 

y ! * 0, y i(G, E) = for all G G { G } and E € { E } 

X, else. 



That is, the target state of a node will equal or 1 if and only if it has this unique steady state regardless 
of variations in the conductances created by transistors in the X state. Note that if the switch-level 
network contains no transistors in the X state the sets { G } and { E } each contain one element and the 
target state equals the steady state of this unique logical conductance network. 

If a logical conductance network cante modeled by an order of magnitude network N then it can 
also be modeled by some other network N' if and only if N s N'. Thus, there is a 1-1 correspondence 
between logical conductance networks and equivalence classes of order of magnitude networks. Logical 
conductance networks can be viewed as abstractions of order of magnitude networks, where those aspects 
of the structure which have no important effects on the behavior are ignored, as was shown in Theorem 
3.3. The logic state y provides all of the information we require about the steady state voltages for the 
corresponding class oforder of magnitude networks. 
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3.7 Properties of Logical Conductance Networks 

With the analogy between logical conductance networks and classes of order of magnitude 
networks, we can derive some simple rules for interconnections of logical conductances. First, if logical 
conductances a and b are combined in series, then 

^series = ""/*>, bX 
where logical conductance values are ordered 

0< Yi <y 2 <...< V 

To show this, suppose logical conductances a and b are modeled by resistors of conductance a and b 
where a and b are rational functions in the domain 9. Then the net conductance g is given by 

9 a + b' 
and applying equations 3.10, 3.11, and 3.12 gives 

*S(g) = de&i) + <feg(b) - muidediXdetm « m^(AiCa),<fc«(b)). 

Similarly, if logical conductances a and b are combined in parallel, then equation 3.10 gives 

detto) = *s(a+b) = max(d^aX*8(b)X 

and therefore 

parallel = WKU£ ( a . b )- 

For more complex interconnections, the net logical conductance between two nodes equals the 
maximum logical conductance of all paths connecting them, where the logical conductance of a path is 
defined as the minimum logical conductance in the path. To see this, suppose that for a network of 
arbitrary linear resistors the net conductance between two terminals equals g. Let P denote some set of 
resistors forming a path between the two terminals and C denote some set of resistors which if removed 
would eliminate all paths between the two terminals. If g: equals the conductance of resistor j, then 
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• ifl ir ^ fl ^ 2 9j- 



2 g p c 

That is g is bounded below by the series conductance along any path and is bounded above by the 
parallel conductance through any cut-set These inequalities can be shown by simple applications of loop 
and cut-set analysis. If these resistor conductances are given by rational functions of p, then for 
sufficiently large p, all conductance values will be partially ordered by the degree of their rational 
functions. Therefore 

mindegigd < degHq) < maxde^qX 
P J '-■■€' 

This can be expressed in terms of the corresponding logical conductances as 

mwg: < g < max*. 
P J C * 

where g is the net logical conductance between the two terminals and & is the logical conductance 
corresponding to resistor j. These inequalities hold for any path P and any cut-set C. Suppose P' is the 
path of maximum logical conductance, and C" is constructed by repeatedly removing the minimum 
element in the maximum uncut path until all paths between the two terminals are eliminated. Then the 
minimum element of P ' equals the maximum element of C: 

/m'ngj = g = max*. 
P' J C" 

Therefore g equals both the minimum logical conductance in the maximum path and the maximum 
logical conductance in the minimum cut-set This rule illustrates how the switch-level model describes 
the operation of an MOS network in terms of only the dominant effects acting on each node. It considers 
only the strongest conductance path between two nodes and characterizes this path by a strength value 
from a small, discrete set Much as rules for combining networks of linear resistors may not apply when 
voltage sources are present we^ must be careful m applying these rules when input nodes are present 
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We can also think of the normal nodes in a logical conductance network as forming logical 
capacitances with properties much like logical conductances. Logical capacitances can take on values in 
the set K = { kj k } which are ordered 

0<«!<ic 2 <...<it^ 

When logical capacitances a and b are connected in parallel through a nonzero conductance, they form an 
effective logical capacitance: 

c paraflel = *"*fabfc 

Other forms of interconnection cannot occur, because the capacitors have one side tied to ground. This 
rule illustrates that the switch-level model considers only the largest capacitors connected to a node and 
characterizes the net capacitance by an element from a smafl, discrete set 

As a final set of rules, consider die logical co ndnrtaary network shnwq in Figure \\ rat^ iqigg fl 
normal node ^ connected by logical conductances g x , g 2 , andgj to input nodes in states t, 0, and X, 
respectively. The steady state y i will depend on the relative strengths of these conductance*: 



Fig.ll State Formation in Ratioed Orcato 



*3 

X MZ3 
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1, g 1 > trwigq, 83) 

y » = i °' g2 >Wfl *tel.83) 
X, else. 



{ 



To show this rule, we can model the network by an order of magnitude network containing resistors of 

conductance gj(pX nfo), and §3(0) connected to wofcages V^ &0, and some voltage V such that 

0< V < V dd , respectively. Node ^ will have a steady state voltage 

_ frfo) * v dd * fl^ ? Aft. *<3fr>' v 
Vi(p) " gi(p) + 9 2 (P)+9 3 (p) 

From equation 3.14 one can see that 

' v dd- 9l»92 + 93 



{ 



v i°° = i 0.0, g 2 >gi + 9 3 
V, 



where 0.0 < V < V^. This rule illustrates that the logic state (or 1) is formed on a driven node only 
when the connections to input nodes in state (or 1) clearly dominate over all other connections to input 
nodes. 

Similarly, if a set of normal nodes are connected such that the net logical capacitance of nodes 
initially in the 1 state equals cj, of nodes initially in the state equals Oj, and of nodes initially in the X 
state equals c 3 , then the steady state of any node n x in the set depends on the relative values of these 
logical capacitances: 



{ 



1, cj > max^, C3) 
y i = *| °> c 2 >inax(c 1 ,C3) 
X, else. 



This rule illustrates that the node state (or 1) is formed on a charged node only when the n« 
capacitance of nodes initially charged to (or 1) clearly exceed the capacitance of all other nodes. 
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3.8 Summary 

Logical conductance networks provide an abstraction of our electrical model of switch-level 
networks. The rather intractable definition of the target state in terms of electrical networks with 
parameters and initial conditions ranging over infinite sets can be reduced to one in terms of a finite set of 
logical conductance networks. The steady state of a logical conductance network is in turn defined in 
terms of an electrical network with a unique set of parameters of initial conditions. 
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Fig. 4.1. Thcvcnin-Norton Equivalent Theorem 



N 




L 








The Thcvenin-Norton Equivalent Theorem (15J, illustrated in Figure 4.1, states that for a linear 
network N connected to an arbitrary toad L (i.e. some other network), we can replace N by its Tbevenin 
(or Norton) equivalent network and obtain the same behavior. The Tbevenin equivalent can be found 
for N without considering die nature of the toad L. The voltage vg^ equate the voltage across the port 
when no toad is attached and hence is called the open-circuit voltage. The admittance Yg^ equals the 
admittance across the port with no toad attached and with all (independent) sources in N set to zero and 
hence is called dw zero-stale admittance. In our case setting all sources to zero involves short-circuiting 
die voltage sources corresponding to input nodes. Since we are concerned only with die steady state 
voltages we can consider the admittance to either be a conductance g^gy or a capacitance c^ ev , because 
the conductance values will have an effect only when die node is driven to ks steady $*afr vohage, and die 
capacitance values wifl have an effect only when die node is charged. One can also see that a voltage 
source connected in series with a discharged capacitor is equivalent to die capacitor charged to that 
voltage. 

Logic signals are related to Thevemn networks in two respects, as is shown in Figure 42. First, a 
logic signal describes the total behavior of a logical conductance at a particular node in terms of a single 
source of state, i.e. either an input node or a normal node, and a single network element, Le. either a 
logical conductance or a logical capacitance. Second, a logic signal provides a direct model of the 
Thevenin equivalent of an order of magnitude network. That is, the logic signal at node tu models the 
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Fig. 4.2. Relation Between Logic Signals and Thevenin Networks 



Driving Signal 
?k 



W\r»i 



Charging Signal 




. GND 




• GND 



Thevenin equivalent of the order of magnitude network as viewed at a port with positive terminaf ^ and 
negative terminal GND. Thus we caa prove properties aboot logic signals by demonstrating the 
analogous effect in the order of magnitude modeL 

The relation between a logic signal and the Thevenin equivalent of an order of magnitude network 
is defined as follows. For an order of magnitude network the open-circuit voltage v^y and the 
zero-state conductance g^y or capacitance c^ ey are given by rational functions of p. If we let 

v oo _ lim v ' 
v thev " p _oo v thev 

then the logic signal corresponding to an order of magnitude network will have state 1 if v^y = V^, 
state if vjjgy = 0.0, and state X if 0.0 < v^ v < V dd . The strength of the logic signal depends on the 

degree of g^y or c^y. A driving signal has strength in the set T = { y^ y p }. A signal of 

strength y^ corresponds to an order of magnitude network in which degHQfa ew ) = k. This strength also 
equals the strength of the strongest path in the logical conductance network from the node to some input 

node. A charging signal has strength in the set K = { kj k }. A signal of strength k^ corresponds 

to an order of magnitude network in which g^ = 0.0 and degic^y) = k. This strength also equals 
the size of the largest node connected by some path in the logical conductance network. A null signal 
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represents an open circuit That is, die corresponding order of magnitude network has zero-state 
conductance and capacitance equal to 0.0. As a consequence, the Thevenin voltage is indeterminate and 
is denoted by the logic state X. The null signal serves as an identity element when combining logic 
signals much as the number is used in other domains of mathematics. 

Just as a Thevenin network provides a composite description of a linear network at a particular port, 
a logic signal provides a composite description of a logical conductance network at a particular node. 
However, whereas the Thevenin equivalent of a linear network depends on die complete structure and 
exact parameters of all network elements, a logic signal depends only on the dominating effects at die 
node. The strength of a signal depends only on the strength of the maximum logical conductance path to 
an input node or on the size of die largest connected normal node, and fee state depends only on die 
states of input nodes connected by maximum conductance paais or the states of die largest connected 
normal nodes. As a consequence, finding the logic signal for a logical conductance network will prove 
much easier than finding die Thevenin equivalent of a linear network. 

43 Rules for Logic Signals 

A simple set of rules describe die logic signals resulting when a logical conductance network is 
constructed by a series of primitive steps. These rules wul first be fisted and (hen shown to describe 
analogous effects in order of magnitude networks. 
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Formation 



a. Input node L forms a logic signal of strength y_ and state X:. 

b. Normal node «= forms a logic signal of strength cap: and state y:. 



2. Coupling 

A logic signal coupled through a nonzero logical conductance forms a signal with the 
state of the original signal and with strength equal to the minimum of the original 
signal strength and the conductance strength. 

3. Combination 

Two logic signals can be combined into a single signal as follows: 

a. A stronger signal will override a weaker, and the weaker signal 
can be ignored completely. 

b. If the signals have the same strength and state, the resulting 
signal has this strength and state. 

c. If the signals have the same strength but different states, the 
resulting signal has this strength and state X. 



4.3.1 The Formation Rule 

The formation rule describes the logic signal formed ^hy an isolated mput or normal node. 
Comparing this rule to the Thevenin networks of Figure 4.2, one can see that the logic signal formed by 
input node /: corresponds to a TheveMn network with Vg^ set according to the logic state Xj, and with 
rfegtggjgy) = p, where y is the maximum allowable transistor strength. This resistor was not present in 
the formulation shown in Figure 2.4, but for the order of magnitude model it witt behave just like an 
infinite conductance. That is, it acts as an identity element when resistors are combined in series and as 
an annihilator when resistors are combined in parallel. This avoids the need to add a special strength to 
represent an infinite conductance. The logic signal formed by normal node «: corresponds to a Thevenin 
network with v^ ey set according to the logic state yj and if capj - *\ then de^c^ ey ) = k. This 
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combination is equivalent to a capacitor with one side connected to GND and charged to the voltage 
v thev 



43.2 The Coupling Rule 



The coupling rule defines the effect of connecting a network described by a logic signal through a 
logical conductance to a node. Figure 43 shows how mis rok describes the effect in the corresponding 
order of magnitude network. A driving signal of strength y k corresponds to a Thevenin network with the 
passive element having conductance Qj, where degig^ - k. As was shown in the derivation of the series 
rule for logical conductances, when this element is connected in series with a conductance g 2 , the net 
conductance has degree equal to miridefa^, defajft. By the ordering of signal strengths, the 
connection rule describes mis effect. A charging signal of strength « k corresponds to a Thevenin network 
with the passive element having capacitance c where de&c) = k. This capacitor in series with a resistor 
of nonzero conductance will have a net conductance of ao and a net capactomee of c. By our ordering of 



Fig. 43. The Coupling Rule 
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signal strengths, a nonzero logical conductance strength will always exceed the strength of a charging 
signal, and hence the minimum of the signal strength and the logical conductance strength equals the 
signal strength. 

4.3.3 Combination Rule 



Fig. 4.4. The Combination Rule for Acyclic Connections 
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The combination rule describes the effect of connecting two networks described by logic signals to 
form a single network. First, let us assume the two networks are independent Then this rule can be 
demonstrated by showing an analogous effect when ports of independent order of magnitude networks 
are connected. This involves three different cases as are shown in Figure 4,4, not counting the trivial 
cases where one of the signals is a null signaL 

When two driving signals are combined, mis corresponds to connecting a Thcvenin network with 

voltage vj and conductance gj to one with voltage v 2 and conductance g 2 - The zero-state conductance 
equals the effect of connecting the two resistors in parallel, and therefore 

The combination rule yields a signal with strength equal to the maximum of the two signal strengths and 
hence corratfy characterizes the value of g^. The open-circuit voltage is given by 

Therefore 



{ 



***< = 1 "2°°. ■l<i 2 

«*l°°+#*2 O0 '«2~*2 

where* and fi are positive constants. When9 1 ~g 2> vjj^ will equal 0.0 (or V^if and only if both 
vj and v 2 °° equal 0.0 (or V^ The combination rote yiete a signal with stare equal to the state of 
the stronger signal if the two signals have different strengths, and otherwise it yields state (or 1) if and 
only if both signals have state (or l.) Thus the combination rule correctly characterizes the value of 
v thev Wnen two eh 31 !"*? signals are combined, mis corresponds to connecting a Thcvenin network 
with voltage vj and capacitance Cj to one with voltage v 2 and capacitance c 2 . The analysis of this case 
proceeds much as with the previous case. When a driving signal is combined with a charging signal, this 
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corresponds to connecting a Thevenm network with voltage vj and conductance gj to one with voltage 
v 2 and capacitance C2- The resulting network has a zero-state conductance gj and an open-circuit 
voltage vi. Since a charging logic signal is always weaker than a driving signal, our rule correctly 
describes this effect 

We have shown that die combination rule holds when independent networks are combined. In fact 
a similar rule holds for arbitrary linear networks and hence we have not yet achieved a major 
simplification over more detailed electrical models. If the combination rule is applied only to 
independent networks, however, it can only be used toeonstract acychc networks. In general, networks 
may contam cycles, and to construct these we must create cycles by combining the logic signals describing 
two nodes in the same network. This corresponds to connecting together two ports of the same order of 
magnitude network to form a single port. With logical conductance networks, the effect of mis cyclic 
connection is #ven by simply applying the combination rule to the two signals. Ih other words, die port 
behavior of the order of magnitude network formed by connecting two ports of a network N is equivalent 
to die port behavior of die network formed by connecting die ports of two different copies of N, as is 
depicted in Figure 4.5. This can be motivated intuitively by noting that a logic signal describes only die 
dominating effects at a node and these effects involve offty simple chiBa^rizations of acyclic patiis. This 



Fig. 4.5. The Combination Sale TdrCycHc Connections 
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provides a major simplification over general linear network models, because the Thevenin-Norton 
Theorem applies only when the network N and the toad L are independent In general the Thevenin 
equivalent of a linear network can be found only by solving a system of linear equations. 

To show the combination rule holds for cyclic connections, suppose ports 1 and 2 of an order of 
magnitude network N are connected. If one of these ports is described by a driving signal while the other 
is described by a charging signal, there can be no path in N between the positive terminals of the two 
ports. Therefore this case is just like an acyclic connection. Similarly, if both ports are described by 
charging signals, either there is no path in N between the positive terminals of the two ports, and hence k 
is just like an acyclic connection, or there is a path and the new connection is redundant. The 
combination rule correctly describes bom of these possibilities. The more difficult case occurs when bom 
pom are described by driving signals. We must show mat the new network wiH have a zero-state 
conductance g^^ which obeys equation 42 and a limiting case open-circuit voltage v^ which obeys 
equation 4 J. 

We can show the zero-state conductance obeys equation 42 iising the rate derived in Chapter 3 for 
finding the net logical conductance between two nodes in a togkal conductance network. This rule states 
that the net logical conductance equals the strength of the strongest pam between the two nodes, where 
the strength of a path equals the weakest togkal conductance in die path. Given the correspondence 
between the logical conductance y^ and a resistor wiu conductance given by a tatio&d function of degree 
k, we can apply this rule to find the "degree" of me net conductance between two nodes in an order of 
magnitude network, i.e. die degree of the rational function which gives die net conductance. The degree 
of the net conductance between two nodes must equal the degree of die path with maximum degree, 
where the degree of a path equals the degree of die element in die nam with minimum degree. 
Furthermore, we need only consider acyclic paths, because for any cyclic path we can form a path with 
conductance of greater or equal degree by removing the cycles. In connecting the two ports of N, we do 
not form any new acyclic paths from the positive port terminals to GND, and hence the set of acyclic 
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paths across the new port equals the union of the sets of acyclic paths across ports 1 and 2 of N. 
Therefore 

^sCflthev) = *M*«tol« 92»- 

and the combination rule gives the correct signal strength for cyclic connections. 

We can show that the limiting value of the open-circuit vohage y^°, y obeys equation 4.3 by 
applying an equation for v^y derived in Appendix I by an analysis of multi-port networks: 

_ 9,00- Vn + gyj- fc^ 2 (44) 

v thevW>) - gi (1.0-k 2 ) + 92(1.0-^) X * A) 

All of the terms in the above equation are rational functions of p. The factors kj and k-^ describe the 
strength of the connection within N between the positive terminals of the two ports. If both equal 0.0, 
then there is no connection and the equation reduces to the one for an acyclic connection. If kj(p) equals 
1.0, then all paths from the positive terminal of port 2 to voltage sources pass through the positive 
terminal of port 1, and vice-versa. For values of kj(p)6%lweehi;exclusnfely) 0.0 and 1.0, some paths from 
the positive terminal of port 2 to voyage sources pass through |he positive terminal of port 1 and some do 
not These factors obey the following properties for any positive value of p: 

92<P)*l(p) = 9lO>) k 2<P> 
0.0 < k x (p) < 1.0 

Oj0 < kjfcO < L0 

k^VvjCp) ^ VjCp) 

k 2 (p)-V2(p) < vj0>). 

The second and third inequalities imply that kj and k 2 have degree less than or equal to 0. Let us 
consider 3 cases. 

l.g 1 »g 2 

Then de^k 2 ) = de&Q 2 k l /q l )<0, and hence <feg(1.0-k 2 ) = 0, while deg(l.Q- kj) < 0. Therefore 
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fll (1.0-k 2 )»g 2 (1.0-k 1 ),aiKlv^ ¥ = V! 00 . 
2. g x < g 2 



oo „oo 



An analysis similar to the previous case shows that vjjj^ = * 2 
We wish to show that for some positive constants a and fi 

00 00 * OO 

v thev = ^i+fiy 2 . 
Consider the following four possibilities: 

a).k 1 oo <1.0,k 2 QO <10 

Then 92(1.0- k 2 ) ~ g 2 (1.0- kj), and the desired result will clearly hold. 

b).k 1 oo = k 2 oo = 1.0 

This implies that vj 00 = v 2 °°. in which case 

*£ev - *1°° ■» 0.5v 1 oo ^0Jv 2 ^. 



c). k! 00 < k 2 °° = 1.0 



Then v°^ y = v 2 ° , and 



k i°° B ^-.-^«f«iW^K^>> >0 



P — *:QQ 

Furthermore k^-Vj 00 < vjj^ < vj 00 , and therefore 

O^k^-vj 00 + 0.5v 2 °° < v* v < O^ 00 + 0Jv 2 °°. 

For = 0.5 and for some a such that 0.0 < 0.5^°° < a < 0.5 the desired result must hoM. 

d).k 2 oo <k 1 oo = 1.0 

The analysis of mis case proceeds much like the previous one. 
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Therefore, when gj ~ g 2 , v^ ev must depend on both vj- and v 2 . This completes our proof that 
the combination rule correctly describes the effect of acyclic connection. 

4.4 The Steady State Signals 

Each node tu in a logical conductance network can be characterized by its steady state signal, 
denoted vj, analogous to the Thevenin equivalent of the corresponding order of magnitude network at 
this node once it reaches steady state. The state of this signal equals y j, the steady state of the node. The 
value of this signal can be found by constructing the network by a series of primitive steps according to 
the three rules. This task can be difficult and tedious, however, and must be done separately for each 
node. Instead, we will show that an equation can be derived which expresses the value of the steady state 
signal for each node in terms of the signals formed by the input nodes and normal nodes in their initial 
states and by the steady state signals on other nodes. 

Let us first introduce some notation, much of which will be replaced by more concise notation in 
Chapter 5 when the logic signal concept is formalized into an algebra of logic signals. For a nonzero 
strength value s, the expressions +s, -s, and xs denote signals with strength s and states 1, 0, and X, 
respectively. The null signal is denoted X. Signal-valued variables are denoted with italic characters. The 
vector x denotes the signals formed by the input nodei in state x. That is x % has state xj and strength y^ 
Similarly, the vector y denotes the signals formed by the normal nodes in their initial state y. That is >| 
has state y i and strength capj. The signal resulting when logic signal a is coupled through logical 
conductance g is denoted cple(a, g). According to the coupling rule this signal will have the same state as 
a and strength equal to the minimum of g and the strength of a. We will adopt a convention that 
cple(a, 0) = X, i.e. any signal coupled through a zero conductance yields a null signal 
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The rule for combining logic signals imposes the follow partial ordering on signal values: 



/\ 



+ y, 
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x Vi 
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XK, 



XKj 



/\ 
\/ 



-K\ +K 



That is, for signals a and b, a < b if and only if the effect of combining a and fc equals b. Thus, die result 
of combining a set of signals equals the least upper bound (abbreviated l.u.b.) of the set for this partial 
ordering. This partial ordering will prove important in our mathematical development It provides a 
concise statement of the concept that the logic model considers only the dominating effects acting on a 
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node. 

4.5 Constraints on the Steady Steady State Signals 

We can use the rules of logic signals to derive a set of constraints which the steady state signals must 
satisfy. Let N be an order of magnitude network in which node «j has capacitance c^, initial voltage y v 
and steady state voltage vj. Suppose we were to connect an additional capacitor to this node with 
capacitance c where c ~ Cj, and charged to an initial voltage y, where y c* yj, giving a network N\ 
Then the new network would be equivalent to the old one, i.e. N« N'. Equivalent order of magnitude 
networks are represented by die same logical conductance network, and therefore if an analogous process 
is performed with a logical conductance network, the logic signals should remain unchanged. This new 
capacitor is described by & logic signal 
application of the combination rule to >j and v- t giving 

Similarly, suppose mat N has a total conductance ey connecting node n^ to the voltage source 
corresponding to input node L Then if a conductance e is added between n- x and this voltage source such 
that e ~ ey, the new network N' will be equivalent, Le. N^N*. therefore, if we perform an analogous 
process with a logical conductance network, the logic signals should remain unchanged. A resistor with 
order of magnitude conductance e is described by the logical conductance ey, and hence the effect of ins 
new conductance on «j is described by the signal cpl^x-., ey)» giving 

^ = lu.b.{cpHx-y ey), *,} > cpKx-y ey). 

This result holds even if ey = 0, because cple(x^) - X, and vj^ X, By a similar line of reasoning 

Vj = Lu.b.{cp/£(Vj, gy), Vj} > Cp/<Vj,gy). 
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If vj is greater than or equal to all of these signals, it must be greater than or equal to their least 
upper bound: 

vj > Lu.b.({^} U {c/ > /e(x j ,e ij )|l<i<m} U { c/rfe^. gy) | l<j<n } ) . 

If we define die function/: as 

/j(«) = Lu.b.Jf^} U {cp/^e^ll^j^m} U {cpfc(^gy)(l<i<n}), 

then these constraints are expressed by an equation *|>/j(i$. If we then let / denote the function 
yielding a vector with ith element equal to the application of f x to the argument, then the set of 
constraints for all node signals can be expressed by a single equation f>/(r). 

We can show in feet that r = /(i* In other words, the steady state signal on each node equals the 
least upper bound of the signal formed initially on the node and tf^ signals on ao>^nt input and normal 
nodes coupled through the logical conductances between them. To show this, let /denote the vector 
given by v =f(r\ and let r and i denote the vectors of signal strengths for the signal vectors v and r'. 
That is, for any value of i, rj equals the strength of vj and r'j equals die strength of v'j. Suppose first that 
for some nj, v- ( is a charging signal, Le.it has strength rjGK. Then for any rt such that fc>0; 

which implies that »j=/|. If.onteotrKThano^j = Oforallj.ttenmjde^ 

theBefow ^ ~ *'i = J\- Ne3st ' «W"»» *i » a dr*wnfrj%n*4 ie. it has strength r^e T. If *j>V| men 
either ^ > r\ or *, has state X and v'j hat state or 1. IT^>^ then for any normal node nj. the 
constraints on the logic signals imply the following constraints on die signal strengths: 



r i > r 'i > 



OTMCgyfj) 



Therefore, if ij > rj, the first inequality can hold only if gy < r'| < rj. and if n < rj, the second inequality 
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can hold only if g^ < n < r y Similarly, for any input node j, 



r s > r' : > *M^, yj } >'«■ 



"i ' *r =- 



In other words, all logical conductances connected to «j have strength less than rj, the strength of die 
steady state signal. Therefore, since any path frriin «$ to anttputaode must pass through one of these 
condoctaaces it cannot havestrtagfli equal to *he steady stale signal stmigth. This contradicts the Act 
that 4he strength of the steady stale signal equals the maximum of these path strengths, and dierefee rj 
cannot be greater than r' iv Thus we can assume that Tj =tLr*jfor:att4 I*K^«upf»ie r^ 3 r'j but the state 
of ^equals X and me state of v'j equate & For any node 1^ if f^»* r ulen *j > w&^gg, rj>*= r- r 
Furthermore, v'j > cple{y y gy), and therefore vj must have state y j » &■ By a«mflar Hue of reasoning. 
any node L for whiehey equals r$ (it cannot be greater^ = 0. In oo^ words, any node connected to 14 
by a logical conductance of strength greater than or equal to rj must have state 0. Let us look at what this 
implies in the corresponding order of magnitude network. By KirchorTs Current Law, the net current 
flowing out from node nj equate 0.0: 

Kvj-xpey + SKvvpgy = 0.0. 

In the above.equation-au^termsare rational functions of p. Assuming rj . a y^ we can divide through by 

p k and take the limit as p approaches infinity: 

m q.. n n... 

P^»( j f 1 (v '' x ^ + j !i te "^H il0 ' ' 

Let us look at the individual terms in this equation. For any j such that <feg(gy) < k, 

p-,89 l Jy * ,,■■■■* ~ p~+,M p* 

For any j such that de&%£> k, yYs= 0, and therefore vy° = OjO. Furthermore, by our assumption 
7j s X, and therefore Vj 00 >OJ0. This implies that 

■"„(v«j)?-"(«i*' : »i^ ft "oo* >(u '- 
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A similar line of reasoning shows that for any j if de^e^)< k, 



and if *«(ejj) > k, 



*V (vj-x^ > ao. 



These results would imply that die above summatioa is greater than 0.0, which is not possible. Therefore 
Vj cannot have state X when v'j has state 0, and a similar line of reasoning shows that v-. cannot have state 
X when Vj has state 1. This completes our proof that? = /(»)* 

We have shown mat die rules for logic signals imply a set of constraints which must be satisfied by 
the steady state signal for each node «j: 

vj = l.u.b. ( { n } U { c P Hx y ey) | l<j<m } U { epldy r gy) | l<j<n } J . (4.5) 

Note how we used the fact that the combination rule holds for cyclic as wen as acyclic connections to 
derive the constraints on the steady state signal for each normal node in terms of the steady state signals 
on other normal nodes. No such set of constraints can be given for the Thevenin equivalents at the 
different ports of an arbitrary linear network, because cyclic combinations cannot be dealt with and 
because mere is no analogous partial ordering. Thus, we have taken a major step away from electrical 
network concepts. 

4.6 Specification of the Steady State Signals 

The constraints on the steady state signals given in equation 4.5 can be expressed by a recurrence 
equation of the form v = f(v). This equation in itself, however, does not provide a unique specification 
of Ae steady state signals, because it may have multiple solutions. For example, the networks k Figure ' 
4.6 contain no input nodes, and assuming the nodes have initial state X, me steady state signals of me 
nodes in both cases should equal xjcj. Suppose in the- first example, however, that we let 4 equal +iy, 
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Fig. 44. Networks with Extraneoas Solutions 
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Then, since ^ equals y 2 , and +y 2 = cpH+y 2 , y 2 )> *« signal a * tisfles ** recurrence reUtion a " /W. 
as would any signal such that xkj < a < xy 2 . The second example shows a similar case in which two 
nodes have mutually dependent signals, and hence any sehiliim such that x« x < a l = a 2 < xy 2 *<> uld 
satisfy the recurrence relation a = f(a). These examples demonstrate mat the equation a = f(a) may 
have solutions in which some nodes have signals greater than the steady state signals because of 
extraneous cyclic dependencies allowed by the recurrence relation. In the second case, a single logical 
conductance creates cycljc dependencies in both directions. These extraneous solutions have no physical 
significance; they are artifacts of the simplified view of electrical networks provided by logic signals. 

Fortunately, we can exclude these extraneous solutions by considering only the minimum vector 
which satisfies die recurrence relation. Let v denote the minimum solution of the equation a = /(a). 
That is, v = f(v) and for any a such that a = /(«), v\ < flj for all i. We will show in Chapter 6 that such 
a unique minimum solution must exist The minimum solution depends omy on the initial signals on the 
nodes and the consistency constraints impo»d by the recurrence relation. Therefore, it seems quite 
reasonable that the vector of steady state signals v should equal the minimum solution v. By a 
generalization of the technique used to show that v = f(v), we will show that v = v. This gives us a 
complete specification of the steady state signals in terms of a simple set of operations on logic signals. 
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First suppose for some node n^, vj is a charging signal. Let A equal the set of all nodes connected by 
some conducting path to /ij. Then Vj describes the result of charging sharing between these nodes: 

v- = l.u.b.{^|njG^}. 

We can show that v'j cannot be less than v- using this equation. For any th and n k in A such that gj k > 0, 
v'j > Q»/^gj k , v' k ) = v' k > cpljfy, v j) = Vj, 

and therefore v'j = v' k . For any node /t G A, we can show by induction on the length of the shortest 
path to «j mat v'j = v':. Furthermore v'j > y- for all th which implies that 

v'j > Lu.b.{>j|nG^} = + r 

Therefore Vj = v'j for any charged node. 

Next, let us consider driving signals. Let r and r' denote the vectors of signal strengths giving me 

strengths of the elements of v and r , respectively. For any node /ij, if V| > v'j, either Tj> r'j of v-, has state X 

and v'| has state or 1. Let s G T equal a strength value where for all j such that n > s, n = r'j, but for 

some i, t % = s>r'j. Let A and B denote the sets 

A = {« | |s = r i >r* i } 
B = N-A. 

Then for any i^ € A and * G & 

s > r'j > m/Vitejj. r j) 

Therefore, if Tj > s, then r'j = n and the first inequality can hold only if gy < s. Similarly, if n = s, and 
since «j g 4 then r'j = Tj and the first inequality can hold only if gy < s. Finally, if n < s, the second 
inequality can hold only if gy < s. Furthermore, for any input node t 
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s > r'j > mW*fryJ - ey 

and therefore ey < ij. In other words ait logical conductances connecting nodes in A to input or normal 

nodes outside of A have strength less than s, the strength of the steady state signals for all nodes in A. 

Any path from a node in A to an input node mast pass throygh one of these conductances and hence 

cannot have strength equal to the signal strength, but this contradicts the fact mat the strength of the 

steady state signal equals toe maximum of these path strengths. Iherefore, there can be no strength value 

s satisfying our requirement and r^ = r'j for all i. Next, le£S equal a strength value where for all j such 

mat n > s, v: = v ':, but for some i, q = sand ?j > Vj. Let A and 5 denote toe following sets: 

A = { «j | Tj = s, Vj has state X, and v'j has state } 
B - N-A. 

Then for any *j€ A and «: € *if gy > fj, then r-. > miiig^, rj) -i^S F&theimiore, v*- > ep$fci,t$, 
and therefore v= must have state y ■ = %. ' Hy a similar fine of reasoning, for any n^GA and any Mprit 
node t for which ey equals r^(it cannot be greater), . t = 0. Thisshowf that for* any normal node in Bot 
for any input node, if it is connected to a node in A by a logical conductance of strength greater than or 
equal to s, it must have state 0. Let us look at what this implies in toe corresponding order of magnitude 
network. By Kirchoff s Current Law, if we form a cut-set consisting of all branches connecting nodes in 
A to nodes in B and to voltage sources corresponding to input nodes, the net current through this cut-set 

must equal 0.0: 

m 
2 2(v i -x j )e ii + 2 2 (vj-vpn^ = 0A 
n^GAj=l n^GAmGB 

In the above equation, all terms are rational functions of j£ f8r s * - : f j*;^ can divide through by p K and 
take the limit as p approaches infinity, 

'*" / 2 2(v i -x i Aj +22 (v r v:)S{i\ = 0.0. 
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Let us look at the individual terms in this equation. For any i and j such that <feg(g ;:) < k, 

*V (vj-vpft = (Wj 00 -*: 00 ) *» H = 0.0. 

p -» oo ' r p K i j p _> oo p* 

For any j such that «j G B and tfegtey) > k, y j = 0, and therefore Vj°° = 0.0. Furthermore, we have 
assumed that for any nj 6 A, y { = X, and therefore Vj 00 > 0.0^ This implies bat 



(Vj-v:)^ = (v: 00 -^ 00 ) ■*" $f > 0.0. 
J p k * J p-»o©pk 



p — ♦ 00 



A similar line of reasoning shows that if <feg(ey) < k, 



/on 



P -» °° - p 

andif<feg(eji)>k, 



(VjXjhf = 00. 



p-+oo J ^pi 
These results would imply that the above summation is greater than 0.0* which is not possible. Thereibre 
the set A must be empty, but a similar result holds if A k defined as 

A = {^(fj = s, vjhas state X, and ?'jhas stale 1}, 

and hence there can be no such strength value s. For all driven nodes ^ = v'j, which completes our proof 
that the vector of steady state signals r will equal die unique minimum value satisfying the recurrence 
relation r = J\f). 

The constraints on die steady state signals given by equation 45 give us a specification of die set of 
steady state signals without any reference to an electrical model, as long as we consider only die minimum 
solution. This will allow us to develop a method of computing the steady state of a logical conductance 
netwoA without evaluating any electrical networks. 
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4.7 Signal Blocking 

Informally, we can view a driving logic signal as describing the combined effect of those input 
nodes connected by a path of maximum strength. Under certain conditions, however, this informal view 
may not be entirely accurate, because of a phenomenon we call signal blocking. Consider, for example, 
node «2 in the two networks shown in Figure 4.7. In both networks mere is a path of strength 
miiirfY- Yi) = 7i to an input node in state 1, and a path of strength mhiy^ Yj) = yj to an input node 
in state 0. Our informal view would then suggest that in both networks ^ = l.u.b.{ +Yj, -y± } = xy^, 
and hence node /^ has a steady state X. While this analysis yields the correct result for the second 
network, it fails for the first In the first network node «j will be driven to by the signal -y 2 and hem* 
« 2 will also be driven to by the signal cpH-yj, yj) = -y^ This example demonstrates that when the 
paws to input nodes are not independent, some signals may be blocked along a path by a signal of greaier 
strength. For example the signal +yj is blocked at nj by the signal -y 2 , and hence node * 2 is unaffected 
by tim signal Our iaformal view does not take such a possib*ty into account, although our more format 
method does. Note that this phenomenon can be ignored in restricted logical conductance networks, 
which are defined as networks in which any logical conductance between two normal nodes must have 
strength y . This class of networkscan be used to analyzing resfricted switth-level networks as were 

Fig. 4.7. Signal Blocking Example 
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defined in Chapter 2. In restricted networks, the strength of any path to an input node wHl be 
determined by the strength of the final logical conductance in the path, and hence the informal rule will 
correctly describe the precedence between signal paths. 

Signal blocking creates difficulties in the formulation of a method for computing the steady state 
signals, but a resolution of this problem leads to an efficient method for computing the target state of a 
switch-level network. 

48 CMrinsiM 

The simple set of rules for logic signals lead to a specification of the steady state signals: 

n = Luaj.({^} U {q*4x^e£\l<i<m} U {g^g^Hr^n}), (45) 

whereto elements of r must be the mmtt^^ Since 

the steady state y j equals die state of the sigi^ v^ this gives m a specifka<k» of the steady state for any 
logical conductance network in terms of a staple set of operations on logic signals rather than in terms of 
an electrical model Thus we have completed our Lansraon form a dicuk-oriented view of me 
switch-level model to a more abstract logical view. From trus point on ward we will deal only with these 
logical concepts. Tlieonter of inagnn^ 

The logic signal formalism takes advantage of our simplified model of me electrical behavior of 
switch-level networks. In general linear networks, the state variables (Le. node voltages) are to some 
degree dependent on the complete network structure and the exact parameters of all connected network 
elements. As a consequence, computing the steady state involves solving a set of simultaneous linear 
equations. With logical conductance networks, on the other hand, the signal on a node generally depends 
only on the network parameters along a small number of paths, and consequently die state can be 
computed much more easfly. 
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Many logic simulation programs use values (generally caBed "sW«C) which descrfte both logic 
state and relative precedence. Such values correspond closely to logic signals. For example the high 
impedance or M state described in Chapter 1 corresponds closely to the null signal X, because both 
represent open circuits and are overridden by any #her sta|e. .Most MPS logic simulator* |9. 34, 51 
encode both Jog* state and strength into a single va|u$ and sim?^.^ 

and propagating these values according to some set of rules. These rules, however, often do not display 
me consistency, accuracy, and gejjeia^^ 

outlined in mis chapter, on the other hand, caif - ; Jh«g 3 'fioMpim^K^i. .«|*>, .J». algebra witii operations 
corresponding to signal oBi&uulfa_.wi,&n&ak and algorithms can be dew^ped which perform 
simulatic^byso]^ , !■; 
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5. An Algebra of Logic Signals 
5.1 Introduction 

Shannon showed originally thai Boolean algebra could be applied to die study of relay networks, 
and later this algebra was applied to logic gate networks. With KIOS networks, however, Boolean algebra 
does not suffice. Although the nodes assume Boolean (of ternary) logic states, the network conditions 
which create these states depend on die relative sizes of conductances and capacitances. Therefore, the 
conductances and capacitances cannot simply be characterized by Boolean values. Furthermore, in a 
voltage-driven logic, one must consider the state of a signal source as well as die strength of the 
conductance path to it We will develop our own algebra based en the rules oflogic signals which retains 
much of the simplicity of Boolean algebra while allowing a more detailed description of the network. 

52 General Definitions 

Let us start by defining some terminology, most of which is consistent with standard mathematical 
practice. 

For domains 9 and S' and a function/^ -» 9' define die poinlwise extension to vectors of size n, 
Le. /:9 n -+g' n as the vector resulting from the application of the function to each component of the 
argument The pointwise extension of a function with more than one argument is defined as die vector 
resulting from the application of die function to the corresponding components of each argument The 
pointwise extension to matrices of size nXm is defined similarly. For example, in linear algebra matrix 
addition is the pointwise extension of scalar addition. Whenever a scalar function is shown applied to 
vector or matrix arguments, its pointwise extension is implied. 



-97 



For a domain 9 apatiaLonteriag < is extended to vectors ki 9° and mairieeSin9 n><in by saying 
that one vector (or matrix) is < anoAer if tlii entering holds for«ach pair of corresponding elements. 
Obsero that the pointwise extension of the feast upper bound operation tot some partial ordering equals 
me feast upper bound operation for the extension of the partial ordering. 

For a domain 9 with partial ordering < and a domain 9' with partial ordering <\ a function 
/.H — » 9' is monotonic if 

a<b -* /(a></M 

A function of more than one argument is monotonic if it is monotonic for each argument, Qne can easily 
see that the composition of monotonic functions must be monotonic, as is the pointwise extension of a 
monotonic function. 

5.3 The Algebra of Signal Strengths 

Logic signals have strengths in the finite set 

which are totally ordered: 

0<« 1 < ... <k_<Yj< ... <y- . 

When two signals akribtae, the ! resulmig signal rta* streh^ 

strengths. When a signal is coupled through a logical conductance, the resuhang signal has strength equal 

to the minimum of the signal and conductance strengths. The binary operations T and 1 are defined to 

give the maximum and minimum of their arguments, Le. 

a|b = max(a,b) 
a i b = mirfa b). 
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These operations have properties similar to addition and multiplication, respectively. We will refer to the 
minimum of a set of values as the product, and me maximum of a set of values as the sum. Signal 
strengths can be described by as algebraic system frtlbyp) whkh satisfies the following (not 
exhaustive) list of properties: 

la. (£1,0) is a monoid: 



a. at(btc) = (atb)Tc 
m. a|0 = ©t« =* 


(closed) 

(associative) 

(0 is the identity) 


b- (% i. J J » a monoid: 




i. a*(bic) = (a^b)ic 
»• *iT^ = -r>|a = a 


(dosed) 

(associative) 

(^Btbekfcntity) 


c. OisanaaamflatorfbrX: 
ai0 = 




L | is commutative and ktempofeent: 




L atbs=bt« 

i. afa = a 





3. i distributes over |: 

ai(bTc) = (alb)T(aicX 

4 - tf *!• *2- • ' • • *fr * • • * * countable sequence qf elements in 3", then 
a lT»2t--taiT--- exists and is unique. Moreover, associativity, 
commutabvity.arikieBTpotenttar^ 
(sequences of values combined with f .) 



5. i distributes over countably infinite sums as wefi as finite 
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Hopcroft. and UUman $1 H does not jqm a rii« tew v«, because most elements lack inverses under 
the sum operation ,f, Futifaetmm, &* mtMti&mtettm* a*** not t$ jdempotent 

We can define vectors and matrices of strength j^pe^f^ch mSk^pm us to oescribe the network 
structure and signal strengths in terms of matrix equations. If the pointwise extension of T (A) is applied 
to two vectors, the resulting vector will have elements equal to the componentwise maximum (minimum) 
of the two vectors. The pointwise extension oft then has properties similar to vector and matrix addition. 
We can define a matrix product • for strength vaftKS where f is analogous to addition and J. is analogous 
tomuttiplicatioB* If A =;$•$; jfeMh 



a ij = [< b iki c kP- (51) 



The properties of closed semirings also hold for the algebra when extended to square matrices. 
That is (* nXn , T, *, 8, 1) also forms a closed semiring, where y nXn denotes die set of nxn matrices 
over J, denotes the matrix whose elements are all 0, and I denotes the matrix with y p '$ on the diagonal 
and O's elsewhere. 

For any closed semiring, the closure operator, denoted , is defined to give the reflexive, transitive 
closure of its argument For example, if AG* 11 * . 

A* = I T A T A-A T A»A?A | ••- •*■ ^oo** (5 ^ 



This operator has me property mat 



A* = I T A*A* (53) 



The closure operator will be useful for analyzing the conductance paths in a network. 
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The operations J and • obey some of the properties of closed semirings even when applied to 
rectangular matrices and to matrices of different dimensions. In a matrix equation, as long as f is applied 
only to conformable arguments and • is applied only to arguments with me number of columns in the 
first equal to the number of rows in the second, dien: 

1. t and* " are associative 

2. t is commutative and idempotent 

3. • distributes over f. 

These properties will allow us to manipulate and transform matrix equations in the strength algebra much 
as in more traditional matrix algebras. 

The set J with the ordering < form a complete lattice {35] with T and J. serving as the least upper 
bound and greatest lower bound operations. Hence the set of vectors f D with the extension of < to 
vectors also forms a complete lattice. For this lattice the pointwise extensions off and i serve as the least 
upper bound and greatest lower bound operations, respectively. As a consequence, both J. and t must be 
monotonk functions, and therefore must be as well. Only the most elementary aspects of lattice theory 
will be used in this presentation. 

One further function over strength values will be required to express the ability of a stronger signal 
to block a weaker when both signals are described by their strength values. This will prove important 
when the equations in the algebra of logic signals are factored into equations in the algebra of signal 
strengths, as will be described later. Define the function Woe* : Sx f -* 3 as follows 

W«rt(a,b) * J. A a ^ b (5.4) 

la a<b 

We can apply the pointwise extension of this function to vector arguments. The function block is 
monotonk in its first argument and antimonotonic in its second. That is, if a < b then for any c 

Woc*(a,c) < block (b,c) 
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block(c,a) > block (c,b). 

As a result, it satisfies the following properties with respect to the least upper bound operation f: 

Woc*(aTb,c) = block (a, c) T block (b,c) 

block (block (a,b),c) = block (a, b T c) 

bTWoc*(a,b) s bfa. 

Furthermore, serves as an identity element for the second argument: 

block (a, 0) = a, 

and block is idempotent 

block (a, a) = a. 

This list of properties is by no means exhaustive. In fact all of the above properties hold for the identity 
function of the first argument. However, mey will allow usio manipulate and solve equations involving 
the function block . 

5.4 The Algebra of Signal States 

The network model allows node states in the set { 0, 1, X }, with and 1 representing the Boolean 
logic states and with X representing an undefined of erroneous state. We have also introduced a new 
value _L representing a null state. This value will only be associated with the null signal which has 
strength 0, indicating that this signal is devoid of state. 



102 



The set of signal states is denoted T*= { _L, 0, 1, X } and the elements are partially ordered as 
follows: 



A 

A 

\/ 



The least upper bound operation for this partial ordering, denotedU, has thefollowing function table 

U _L 1 X 

X X 1 X 

X X 

1 1 X 1 X 
X XX X X 

The rule for combining logic signals of equal strength is that the resulting signal will have state equal to 

the original states if they were equal and state X if they were not Assuming the value JL is only 

associated with signals of strength 0, this rule is expressed by the operation |_l. That is, if two signals with 

equal strength and states x and y are combined, the resulting signal will have state x U y. The operation 

U is termed the consistency operation, because for nonnull arguments it gives a "proper" value (0 or 1) 

only if the arguments are both proper and equal. OfteiwispUgjxe^aaeir^ya^X. 

The set Twith the ordering <C forms a lattice. Hence the operation LI obeys fee following 

properties: 
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1. yUy = y (idempotency) 

2. xUy = yUx (commutativity) 

3. zU(xUy) = (zUx)JJy (associativity) 

4. x<y=*zUx<zUy (monotonkity) 

The set of signal state vectors T n along with the extension of < to vectors also forms a lattice with the 
pointwise extension of U serving as the least upper bound operation. Hence, the pointwise extension of 
U also satisfies the properties listed above. 

The lattice <^ <> has both the structure and interpretation of the flat lattices used in denotational 
semantics of programming languages (35, 48]. The ,f proper" values and 1 are incomparable, while the 
bottom element J_ represents an "underdefmed" or null value and the top element X represents an 
"overdefined" or erroneous value. 

5.5 The Algebra of Signals 

55.1 Signal Values 

A logic signal is represented by a pair of values <s, y> where s € 3 is die strength and y G fis the 
state with the restriction that y = J_ if and only if s = 0. That is, only a null signal can have a null state. 
Logic signals form a set 

JL = {K l ....,K f y l ....,y p }x{Q.l.%} U {<0,_L». 

The null signal <0, _L> is denoted X. The expressions +s, -s, and xs denote signals with strength s and 
states 1, 0, and X, respectively. The symbols +, -, and x can be viewed as denoting unary functions 
mapping strength values to signals of a particular state with the convention that +0 = -0 = xO = X. 
These functions can be extended pointwise to vectors and matrices as welL 
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Signal-valued variables are written with italicized characters, while strength and state-valued 
variables are written with normal characters. 

Let«a> denote the state of signal a and H at denote itsstrength. These symbols denote functions 
<•>: JL -> t: and I . i:A -» S. 

SJS2 Signal Combination 

The binary operation V is defined to describe the effect of combining two signals. This operation is 
defined in terms of the strength and state of the resulting, steal: 

iflVfrl = IfflttH 



{ 



<d>, lal>lftl 



This operation provides a formal statement of me rules mat a stronger signal wil! override a weaker, and 
mat signals of equal strength wfll combine to form a signal with the same strength and with state equal to 
die least upper bound of the two states. Note that this operation has only a distant relation to the Boolean 
"or" operation which is sometimes denoted V. 
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The operation V defines the following partial ordering among signal values: 



■y P + y P 



x Vi 



"Vi + Vi 



X V2 



XYl 



-Yl + Yl 



x Kq 



XKi 



-*1 +K 1 

\/ 

\ 

That is a < b if and only if a V 6 = b. With this partial ordering we must maintain a distinction between 
the terms "greater" and "stronger". The signal xs is greater than -s or +s but not stronger. 
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The set of signal values J. with the partial ordering '< forms a bake with minimum element A, 
maximum element x-jy and with V serving as the least upper bound operation. As a consequence, V 
must be idempotent, commutative, associative, and monotonia 

We have defined partial orderings for the sets of signal strengths f, signal states T, and signals JL In 
many cases the functions from one domain to another preserve these orderings, Le. they are monotonia 
For example, the functions +, -, and x are monotonia because signals with the same state are totally 
ordered by their strengths. Similarly the function I . I is monotonk because signals of different strengths 
are totally ordered by their strengths. The function <>, on the other hand, is not monotonk. For 
example -y x < +y 2 , but «-?]> < <+y 2 > - 

The pointwise extensions of <■>, I • B, and V can be defined and obey all of the properties listed 
thus far. 

5L5.3 Signal Factorization 

As shall be seen, the algebra of logfc signals lacks many of the properties one might desire in a 
mathematical system, such as a total ordering and monotonkity of certain operations. As a consequence, 
we will often factor equations in the algebra of signals into equations in the more tractabk algebra of 
signal strengths to aid the mathematical development 

The domain of logk signals loosely resembles comptex numbers with strength corresponding to 
magnitude and state corresponding to phase. Just as a complex number is characterized by either its 
magnitude and phase or by its real and imaginary parts, a signal is characterized by either its strength and 
state or by its "1" and "0" parts. The factored form of a signal is a pair of strength values (u, d), with u 
indkating the strength with which the signal will pull a node toward 1, and d indicating the strength with 
which the signal will pull a node toward 0. The following tabk shows die two representations of a signal 
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Signal 


Strength - State 


1 part - part 


X 


<o,x> 


(0,0) 


+s 


<S,1> 


(8,0) 


-*- 


<s,o> 


(0,s) 


xs 


<s,X> 


(s,s) 



The factored form Must either have both parts equal or one part equal to zero, corresponding to the rule 
that a stronger signal will override a weaker, and the weaker can be ignored. A signal with state X has 
equal 1 and parts, because X represents a conflict between signals of equal strength pulling a node 
toward 1 and 0. The null signal X, on the other hand, has no ability to pull a node toward 1 or 0. 
Observe that the ordering < between signals is equivalent to the extension of the strength ordering < to 
the factored form representation. The functions PI and L-J are defined to select the 1 and parts of a 
signal, respectively. 

Tal =* {**V «a>*lorX (5.5) 

\ 0, <a> = 0orX 

UJ _ f-tWj ^►atorX (5.6) 

\0, <a> = lorX 

The factored form of a signal can be viewed as describing two signals, one with state 1 (or J_) and one 
with state 6 (or X), which when combined will yield the original signal: 

a = +fVlV-LaJ (5.7) 

The strength of a signal equals the maximum of its two parts: 

lal = ra1|UJ (5.8) 

The state of a signal can also he determined by combining the states of the signals represented by its I 
andO parts: 

<a> = <+ral>U<-laS>. (5.9) 

This identity holds, because at least one part of a signal will equal zero unless the signal state equals X. 
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When two signals are combined with the operation V, the resulting signal has the following factored 
form: 

raVbi = block (r<n t n>\i a y b i) (s.ig) 

LaVbJ = Woc*(LflJTLbJ t laVAI) 

In the above equations the function bhck preserves die restriction that if one part of a signal is less than 
the other (or equivalentiy, it is less than die maximum of the two parts), it must equal zero. For example 

r+Y 2 V- n l = Mx*( T2 T0,Y 2 > = T 2 
L+72V-7JJ .= bhckm Yi-72) = * 

These identities illustrate how the function Mac* entetswaen equations in the signal algebra are factored 
into equations in the strength algebra. 

15.4 Signal Coupling 

To complete the set of operations for manipulating signal values we require a means of expressing 
the effect of a signal coupled through a logkal conductance. Our rule for signal coupling is that a signal a 
win be coupled through a logical conductance of nonzero strength s to form a signal with strength 
tat i s and state <a>. In Chapter 4 we defined die function q>ler.J.Xf-* J. to yield the signal 
resulting from an application of Ms rule. For our forma) dev el opm en t, instead of using a function with 
different types of arguments, we will express signal coupling with the binary operation • over signal- 
values. The conductance s is then represented by the signal xs, indicating that it can conduct and 1 
signals equally well. Define ° as follows: 
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a *b = HrallTbl) V -(UJ|Li>i). (5.H) 

In other words, if two signals are represented in factored form, " is equivalent to the pointwise extension 
of i applied to the corresponding parte of the two signals. The signal a eoupled through a logical 
conductance s.then forms a signal xs«>a. In all cases usetlbeie, one argument of • will have state X (or 
J_.) This formulation takes advantage of the feet B^tfwe|ei*4.^ denote Asset: 

\X X = {xslsGJ}.' 

then the subalgebra (U. x , V, », X, xy_) is isomorphic to the algebra (*, T, 1. 0, y^. Thus when the logical 
conductance s is represented by the signal xs, it has the same algebraic properties. 
The operation ° obeys the following properties: 



1. (J., °, xy J is a monoid: 



i. a,6€^ s *a l, iG^ (closed) 

ii. a*(k*c) = (« p ^°c -n;.;S,(Mo<MfHi-.' 
iii. a ° xy = xy • a = a (xy_ is the identity) 



2. X is an annihilator for •: 
a«X = X. 



Thus the algebraic system (JL, V, °, X, xy_) o/mos/ obeys all of the properties of a closed semiring. In 
general, however, ° is not monolonic. For example +yj < -y2> but 

XYl 0+ Yi = +Yi < "Yi = x Yi°-Y2- 
As a consequence, ° does /io/ distribute over V. For example >h; 



1. Any operation which distributes over the least upper bound operation for a lattice must be monotonia 
For example, suppose a < c. If it were the case that » distributes oyer V then 

b° aV b°c = b°(aV c) - b"c 

which implies that b° a<b°c. The converse statement need not hold, unless the set is totally ordered. 
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x Yl «»(+ ri V -y 2 ) = -f! ^ xti^+Yi V xyf-12 = x l\- < 512 ) 

This lack of distributivky expresses in mathematical terms that signal paths cannot be analyzed 
independently, because weaker signals may be blocked along a path by a stronger signal with a different 
state. This phenomenon was demonstrated with the networks shown m Rgure4.7. In fact me left hand 
side of equation S.12 expresses how die steady state signal on node nj is formed in me first example of 
Figure 4.7, while the right hand side expresses how the steady signal on node nj is formed in die second 
example. This lack of distributivky will cause some difficulty in our mathematical development 

If we restrict our attention to strong transistors in the 1 and state, however, distributivky (and 
hence monotonicky) holds. That is, if 6€ { X, xy } then 

M<iVc) = *•« V b'c, 

and hence the algebraic system ({ k, xy }, V, °, A, xyj does form a closed semiring. In fact, this 
algebraic system is equivalent to die Boolean dosed semiring ({ 9, 1 }, +, •. 0, 1). Tins shows that 
restricted logical conductance networks approach relay networks m their simplicity. 

If one of die arguments to ° represents a conductance value, flien • is monotonic in mis argument 
That is, if b, c € f, and b < c, men 

xb°a < xc«« " 

for any a€ .4. 

The operation • was seen to equal to pointwise extension of J. to the factored representation of 
signals. This leads to the following identities: 

ra»bi = rvum <5.D) 

UfbJ = UJ4LU 
■ a»6I = (roTirtl) T (UJILUJl 

For die common case where <a> or <b> isintheset{X.X}: 
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*«•&* .« I all 1 61 (514) 

5£5 Matrix Operations 

The matrix product operation p is defined to describe the effects of first coupling a set of signals 
through logical conductances and a»en combining them. That is, if c = ^oithen 

Ci = V (^-6|). (5.15) 

The operation o is closed and assbeiiBve. tMerti^^^ebrafe System U nXn , V,o, 0, /) a/raos/ 
obeys the properties of a closed semiring, where the identity matrix / has xy^'s on the diagonal and X's 
ebewhere^^hfle the iefb rtaSrtx consists of aU X*s: ftt genefal SoweVer, 6 does not distribute over V 
and hence is not monotonic. 

In the special case where *€{ X, xy p } nX m , though, 

/te(*V«) = AoiVAoc, 

and tor ttu> restricted case o jsrao^otonic. Tte.&t%te&&m.{&K x^} aX ^ V,o^fti)then^rBW 
aclc^s^ring^ Tbeisfbrftw^ 

A* = lVAVAoAVAoAoAV... =* V W*. (518) 

(KkXoo 

This closure is isomorphic to die Boolean transitive closure. 

If one of die arguments to o represents a conductance matrix, then o is monotonic for this 

argumenL ThatisifB,CG 3 mXn ,andB<C,then 

xBoi < xCo« (5.17) 

fbrany«Gjt a . 
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The following identity follows from the properties of V and •. 

Mo* I = (m»r*T) T (LOLAJ). (5.18) 

For die common case where each element of A is either a null signal or has state X, Le. A € JL x n x m , 

lAoH = M 1*1 J I. (5.19) 

A matrix product can also be factored into its 1 and parts: 

rAokl = bbck(rA>n\tAoH) (5 JO) 

L4oJJ = bhck{lAJ*lU, iA**l) 

The pointwise extension of die function Mock maintains the restriction that each part of a signal must 
equal the strength of die signal or equal 0. 

5.6 Summary 

The concept and properties of logic signals have been formalized into an abstract algebra with a 
domain corresponding to the signal values ;ai^%ra^ opt e r a tkjhs crMsp^g^tolh^ rules for combining 
and coupling signals. Manybfthtpropmierbftogt^ 

smau\ discrete set, and the operations obey many desirable mathematical properties. Even the 
complications arising from signal blocking are reflected by the lack of distributivity in the algebra except 
for a restricted domain. 

The domains and their operations are summarized below. 



Signal Strengths 
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Elements: 



3 = {0,K l ,...,K q ,y l ,...,y p } 



Ordering: 0< kj< . . . < k < y^< . . . <y. 



Operations: 



T 
I 
block 



maximum (least upper bound) 

minimum (greatest lower bound) 

signal blocking 

(f |) matrix product 

closure 



Signal States 



Elements: T = { _J_, 0, 1, X } 



Ordering: 



Operations: 



U 



\/ 



JL 



consistency (least upper bound) 
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Signals 



Elements: J. = {\,+K 1 ,-K 1 ,XK 1 ,...,+y p ,-y xy } 

Ordering: 



Xy P 



- y P + y P 



c Vi 



"Vi + Vi 



x V2 



/\ 



x Kq 



XKi 



-Kl +Ki 
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Operations: 

V signal combination (least upper bound) 

o signal coupling 

o (V °) matrix product 

closure (for restricted domain only) 

Functions from signals to states: 

<■> state of signal 

Functions from signals to strengths: 

II ■ II strength of signal 

m 1 part 

L- J part 

Functions from strengths to signals: 

+ form signal with state 1 (or J_) 

form signal with state (or J_) 

x form signal with state X (or _J_) 
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6. Computation of the Target State 
6.1 Introduction 

The algebra of logic signals along with the related algebras of signal strengths and signal states 
provide a set of mathematical took for farther developing the switch-level abstraction. In this chapter a 
method will be developed for rinding die minimum solution of the steady state signal equation given in 
Chapter 4, which then gives us a method for finding the steady state of a logical conductance network. 
This method is then generalized into one for finding the target state of an arbitrary switch-level network. 
This development will utilize only the logic signal abstraction as expressed by our idgebras, along with 
two equations which were derived from the properties of the electrical model: the definition of the target 
state in terms of the steady states of a set of togkadeoaductance networks given in equation 3.20, and the 
equation for the steady state signals given in equation 4 JS. While we could arrive at the final results more 
directly by utilizing additional properties of the electrical model, this approach demonstrates mat the 
concept of logic signals is quite powerful and self-contained. The earlier work with the electrical model 
was presented only to motivate and justify die more abstract concepts. The issues of implementing these 
techniques with efficient computer algorithms are deferred to Chapter 7. In particular, the matrix 
notation in this chapter is used only for mathematical convenience and need not imply that the simulation 
algorithms involve matrix operations. 

6,2 Toe Target State Equation 

First, let us restate die definition of the target state in terms of the steady states of a set of logical 
conductance networks. The function target (x, y, z) was defined as giving die set of states y which the 
normal nodes would reach if the input nodes were held in state x, die transistors were held in state z, and 
the normal nodes were initialized to state y. With the transistors held in state i die network of transistors 
can be described by a set of logical conductance networks, where a transistor in the 1 state forms a logical 
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conductance equal to its strength, a transistor in the state forms a logical conductance of (i.e. an open 
circuit), and a transistor in the X state forms a logical coiiductance^ither equal to its strength or to 0. We 
have seen that a set of parallel logical conductances can be replaced by a single logical conductance equal 
to the maximum element in the set. Thus Jhe range of ^o^ conductance rjetworjcs corresrwnding to a 
switch-level network in transistor state z can be described by four logical conductance matrices; G" 1 ", 
G"", E rato , and FT". Each element g™"^ of G"*" equals the maximum strength of all transistors in the 1 

' ■■■< V . ?> 

state connecting normal nodes n^ and n-. or equals if no such transistor exists. The elements of G™" 
equal the corresponding values for all transistors in either the 1 or the X state. Each element e mm jj of E 
equals the maximum strength of all transistors in the 1 state connecting normal node «j and input node ^ 
or equals if no such transistor exists. The elements of E™" equal the corresponding values for fran§istors 
in the 1 or X state. More formally, if 7y denotes the followingset of transistors: 

7y = {/ k |(SOURCE(/ k )=n i and DRAlN(/ k )=/ij) or (SOtJRCE(/ k )=nj and DRAIN(/ k )=n i )}, 
then 

g"":: = T SUV, (6.1) 



" ij ^6%^=!^' 



and 



g"":; = T ,str k . (62) 

In both equations, the maximum of an empty set is defined to equal 0. Both G** and G™" are symmetric 
matrices with elements in the set { 0, y % , ...,y p }. In general these matrices are very sparse, because 
each node is connected to only a limited number of other nodes. Similarly, if T '^ denotes the following 
set of transistors: 

7"y = {/ k |(SOURCE(/ k )=/J i and DRAIN(/ k )=»j) or (SOURCE(/ k )*= »j and DRAIN(/ k )=n i )}, 
then 



* str k , 



e^.. = J str^ (6.3) 
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and 

E""" and E 1 "" are n x m matrices with elements in the set { 0, y x v }. These matrices may be less 

sparse than G""" and G"", because some input riodes (e.g. VDD, GND) serve as source and drain nodes 

for many transistors. The network may also contain transistors connecting input nodes to one another, 

but these have no effect on the logical behavior. To compute the target state function, we need only look 

at the network of logical conductances described by these four matrices, without considering the 

transistor configurations or states which give rise to them. 

The presence of transistors in the X state implies that the actual conductance matrices G and E lie 

within the ranges: 

G^ < G < G— (6.5) 

E^ < E < £■". 

Let { G } and { E } denote the following sets of matrices 

{E} = {E|e jk = (^ jk (p) W e jk =:e~ jk (p)} 

{G} = {G|g jk = g- ta jk (p)org jk = g- jk (p), andg jk =g kj }. 

Note that these definitions of { E } and { G } are equivalent to those given in Chapter 4 in terms of the 
sets of order of magnitude network matrices E and G as were defined in equations 3.7 and 3.8. 

For logical conductance matrices G and E, let y<G, E) denote the steady state of the logical 
conductance network with these values of logical conductances. Equation 3.20 defines the target state of 
a node in the switch-level network as 



{ 



1, yj(G,E)= 1 forallGe{G}andEG{E} (3.20) 

y ' * 0, yj(G,E) = forallGG{G}andEG{E} 

x, else. 
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That Is the target state equals 1 or if and only if it has tins unique state regardless of the conductances 
formed by transistors in the X state, and otherwise it equals X. This equation can be expressed more 
concisely using the consistency operation U for logic states: 

y = U 7(G,E). (6.6) 

1 {G},{E} M 

This equation shows that the target state function for an arbitrary switch-level network can be determined 

by computing the steady states for the logical conductance networks represented by all possible matrices 

G and E in the sets { G } and { E } and then combining these states with the operation U. This reduces 

the problem of computing the target state equation for a switch-level network to one of computing the 

steady stele of alog^conductance«etwork. Asshallb^see^ ! y^,«»eme^o^fer computing to i stea«fy 

state of a logical conductance network can be generalized into a method for directly computing ine target 

state of a switch-level network. 

63 The Steady State Signal Equation 

Let us turn our attention to computing the steady state y for a particular set of conductance 
matrices G and E. 

The rate for forming Jogk signals is that an input node forms a logic signal wWi state equal to Hie 
node state and strength equal to y^ As was defined In <3iap*er^^**tor * denotes the set of signals 
forn^ed by me input nodes in «tete*: 

<*i> *= 34 <67) 

Similarly, a normal node forms an initial signal with state equal to the node state and strength equal to the 
node size. As was defined in Chapter 4, the vector y denotes the set of signals formed by the normal 
nodes in state y: 
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<j\> = v i («J) 

1^1 = capj. 

In Chapter 4 it was shown that the vector of steady state signals r must be die minimum vector 
satisfying the constraints: 

n = l.u.b.({^}U {cpfc( J r j ,e ij )Jl<j<in}U " {q*(^gy)|l<i<n}|. (4.5) 

for all i. This equation can be expressed using operations in die algebra of logic signals as 

*i = V\ V VCxey-Xj) V Vfrg^vp. 
J j 

If die matrices J? and Care defined as the matrices of togk signals representing the logical conductance 



matrices E and G.U. 



E = xE (6.9) 

G = xG 



then this set of equations can be expressed by a single matrix equation: 

v = Eo x V y V Go 9. (6.10) 

Unlike equations in other algebras, in which all expressions involving the, dependent variables can be 
move to one side of the equality, we have no inverse operation in the signal algebra to peflr&canceMatian 
of the left hand side. Hence the equation for the steady state signal m*^j^ expressed ^ a recurrence 
relation. It will be shown shortly that equation 6. JO has a^Qique minimum solution. 
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6.4 Solution for Restricted Logical Conductance Networks 

We wish to find the minimum vector v such that v = /(r), where 

f(a) - £orV;V Goo. 

A general technique for solving such recurrence equations is only known for monotonic recurrence 
equations, i.e. ones in which the recurrence is expressed by a monotonic function. In general, the 
nonmonotonicity of the operation o implies that the function / may not be monotonic, and hence this 
technique does not apply. For restricted logical conductance networks, however, » which normal nodes 
may only be interconnected by conductances of strength y_» the function/ is menotonic. That is, a 
restricted logical conductance network has a condu(^aiK*piauTiC» with each element equal to or y« In 
this case each element of G equals X or x y „ This implies mat for any a and b 

Go(aVA) = Gom V Gob, 

and therefore/ is monotonic. The following theorem, a special case of one given by Scott [35], shows 
how to solve such an equation. 

■'■- ■■'■■ - ' ' _ ■'• B 

Theorem 6.1. 

For a monotonic function/: jl n -+ -4." the equation 

a = m 

has a unique minimum solution given by 

k— » oo 

where denotes a vector of all X's, and the superscript k denotes k applications of the function /. 
Furtoermore,ttwlin«tw& 
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Proof of Theorem 6.1: 

Consider the following sequence 

a. m. firm), .... /to, ... 

First we win prove that its limit exists. Clearly 0</(0, and by the monotonicity of/ one can prove by 
induction on k that 

Hence the sequence is nondecreasing. Any striafy increasing sequence in Ji n would have a length of at 
most n-M|. Therefore for some j < n-M I, A* = /* +1 to a«! then for any k > j, fa = /*(0. 
From this we can see mat rf - * must be a solution, because 

RnaDy, suppose for some «,«=/(«). Starting with the basis* > 0, we can prove by induction on k mat 

« = /*« > /*** 
Therefore 

k— » oo 
Thus rf""" is the unique minimum solution of the recunrax eqUationJ 

This theorem follows as a special case of a theorem proved by Scott [35] regarding the least fixed 
point of a continuous function on a continuous lattice, where "continuity" here is only distantly related to 
the continuity of real analysis. In finding the minimum solution of the recurrence equation a = /(«) we 
are computing the least fixed pea* of the fenction/. Any firate lattice » continuous, as is any monotonfc 
function on a finite lattice. Scott shows that a result similar to equation 6.11 holds for any continuous 
function on a continuous lattice. Finite convergence, however, may not be guaranteed for functions on 
infinite lattices. 
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The process outlined in this theorem of computing the sequence /(fl), f(f(0f) f (<9, . . . 

Until it converges corresponds to a straightforward relaxation algorithm. It starts by setting flj to X for 
each node and then performing relaxations of the form a- «—/"}(«) until it converges, i.e. any further 
relaxation steps would not change a. Since we have expressed the method in vector form, each 
application of/ corresponds to applying relaxation computations at all nodes simultaneously. 

For example the logical conductance network shown in Figure 6.1 models an nMOS Nand gate and 
pass transistor with all transistors in the 1 state. The extra self-loop has been added to n$ to intredwee the 
possibility of extraneous solutions. Equation 6.1u"canbe< 



v l = +T1 V x ?2 ° *2 V x n ' "3 

"2 = ~n v *T2**1 

vj = xyy-^ Vxy 2 «»*j. 



The above equations have the charging signals y left out for simplicity. The relaxatkMi method gives toe 
following sequences: ,^ ; "- 



^: X m -Y 2 t -T2 , -T2 

a|: X -«y i:: -y 2 -i 2 -Y2 

ay X X +Ti -T2 ~T2 



ww 



F4 6.1. Restricted Logical Conductance Network Example 



I 



1 — CD— -^ 



,h 



I 



n 



* 
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indicating that all nodes have steady state signals -y 2 and hence have steady states 0. Observe that an 
extraneous solution never arises because at all times each value <£ is less than or equal to the steady state 
signal. 

For this particular equation the minimum solution can be expressed in terms of the closure 
operation. 



CoroBary U.l. 

For a matrix Ge { X, xy p } n x a , and a vector h € U n , the equation 

« = * V Gom. (6.12) 

has a unique minimum solution given by 

rf* = <7*oJl (6J3) 



ft«tf«fCb*oilai!y6.U: 

The fonctk>n/(«) = b V Go * is monotonic, and an expansion of/*(#) gives 

/*(«) = *VCb(...*VGb(4VCbi)...)* V Gh>k = t V cHoJi 

0<j<k I0<j<k J 



Hence 

k-+oo 
This result shows that for a restricted logical conductance network the steady state signal is given by 

* = G*o(EoxVji (6.14) 

As mentioned in Chapter 5, the transitive closure operation in this restricted case is equivalent to a 
Boolean transitive closure. Element ij of thefmatrix G* equate xy j if nodes n| and * are connected by 
some path of logical conductances and equals X if they are not Thus the matrix G* partitions the 
network into a set of equivalence classes, where /ij and «= are in the^ame class if and only if element ij of 
G equals xy^. Since xy^ is the identity element for • and X is an annihilator, computing element i of 
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G* o b in this case simply involves selecting the elements of * for the nodes in the same class as /ij and 
combining them with the operation V. Thus if we let b equal the vector of signals given by Eo x V y, 
and node nj is in equivalence class Cj^, then 

1 *jec k J 

Observe that all nodes in a class will have the same steady state signal. If we were only concerned with 
restricted switch-level networks with no transistors in the X estate, die computation of the target state 
would be quite simple. 

6.5 Solution for General Logical Conductance Networks 

Unfortunately, the technique outlined for computing the, steady state signal in restricted logical 
conductance networks can fail for general networks! because it may converge on some solution other than 
the minimum. Consider how the proof of Theorem 6.1 relies on the monotonicity of the function /. 
First it sis used to show mat the sequence/ 1 ^ is monotone and hence converges. In fact, it appears that 
mis sequence would converge for equations of the form of 6.12 regardless oHhe monotonicity of/, 
although this has not been proved formally. More importantly, however, the monotonkity of/ is used to 
show that a > j\(fy for any a which satisfies the recurrence relation, and hence the limit to this sequence 
must be the minimum solution. This result does not hold when the recurrence function/ is not 
monotonia 

The network shown hi Figure 6.2 resembles the network shtiwhnf Figure 6.1 except that the pass 
transistor has strength y^, and therefore we no longer have a restricted logical conductance network. 
Equation 6.10 can be expanded as: 
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Fig.6.2. General Logical Conductance Network Example 




CD- 5 



n 



u 



*l = -Tr2 v *T 2 -1 

The above equations have the charging agnate y left out for smiplicity. The miliimiim solution of these 
equations is ^ = -y^ ^ = -y^ and »j t= -^ ghring a steady slate of on aM three nodes, just as one 
would expect. The relaxation method gives die sequences: 



«1 : 


A 


*ri 


"T2 


-T2 


-*2 


«2 : 


A 


~n 


"*2 


-T2 


-T2 


•J : 


A 


A 


*fi' 


*n 


**1 



The sequence converges with the signal x-yj on node ny which is greater than the minimum value ^- 
and therefore finds an extraneous solution with state X on node/13. 

This error arises due to an interplay between the effects of signal blocking (giving a nonmonotonic 
recurrence function), and the possibility of extraneous solutions of the equation. The network of Figure 
6.2 has been contrived to cause this interplay. On the second relaxation step we introduce information 
about the path from VDD to /13 in setting a$ to +yj, and this value is not less than or equal to the steady 
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state signal -yj. Due to the presence of the self-loop at n$, this information wift remain during Hirther 
relaxation steps. When the third relaxation step introduces information about the path from GND to /13 
into ay it combines with the old value of 03 to give xyj instead of -yj. Thus our relaxation method does 
not take the effects of signal blocking into account properly and hence may reach an extraneous solution. 
While this example seems rather contrived, similar effects can occur with more realistic (but larger) 
networks. 

The steady state signal for a general logical conductance network can be found by a method of 
conditioned relaxations which first computes the strength of the steady state signal and then uses this 
information while computing successive^ relaxations tq prevent <j node from being set to a nonnull signal 
weaker than the steady state signal. For example, suppose for the network shown in Figure 6.2, we could 
determine that the steady state signals will have strength y2 on nodes «j and /i 2 and will have strength yj 
on node no. In generating the sequences of values on each node, any time our original method would set 
the node to a signal weaker man die steady state signal, we will instead set it to X. Tliis gives die 
following sequences: 

"1= * * ■ -1% -T2 -72 • • ' 

>2 : X "72 '72 '72 "72 ■■• 

n$ \ \ X ; "71 "*7i ••- 

The signal +yi is weaker than the steady state signal for node «j and hence die first relaxation 
computation at this node will give a value \. As a consequence, tfce signal *yf h never propagated to 
node /13 and will never create an extraneous signal This example shows how the conditioned relaxation 
technique prevents extraneous signals from arising by killing weak signals before they can become 
extraneous. We will now prove formally that this method produces the correct results. 
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6.5.1 Factored Equations 

The method of conditioned relaxations can be derived formally by factoring the steady state signal 

equation into its 1 and parts, giving more tractable equations in the algebra of signal strengths. All 

elements of £and (Tare null signals or have state X, which implies by equation 5.19 that 

lEoxt = l£l*lxl = E-IjcI 
IGorl = IGl-IrR = G-Irl 

Therefore 

Brl = E-lxl T i^l T G'lrl, (6J5) 

which shows that I f I satisfies a recurrence relation of the form Irl = h (I r I), where me function h is 
monotonia Recurrence relations for r>l and LfJ can be derived as well: 

m = block (E • rxi t ryi t g • r>-i, i r i) (6.16) 

LpJ = block (E • LxJ T kM T G • LrJ, Irl). (6.17) 

Thus, if we could determine the value of I rl, we would have recurrence relations of the form 
r>l = /j(rvl) and LvJ = / (LkJ), where boA/j and/ n are monotonfc functions. We will show later 
that H v I, T v\ and L v J must be the minimum solutions of their respective recurrence equations. 

6S2 Recurrence Equations io the Strength Algebra 

The minimum solution of a monotonic recurrence equation in the strength algebra can be found by 
a relaxation method similar to the one shown for the signal algebra, as is proved in the following theorem. 
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Theorem 6.1 

For ammotmk function/:* -♦ S a , the equation 

a = m 

has a unique minimum solution given by 

a- = ■ *" AO) (6.18) 

k— ♦ oo 

where denotes a vector of all O's. Furthermore, this limit will be reached for some k < n-| f |. 



The proof of this theorem parallels the proof of Theorem 6.1. For equations of the form of 6.15, the 
minimum solution can be expressed in terms of the closure operation. 



Corollary 6.11. 

For a matrix G G Jf" x n , and a vector b € J 8 , the equation 

a = b f'G«* 

has a unique minimum solution given by 

a"" = G*»k. 



The proof of this corollary parallels the proof of Corollary 6.1.1. Observe that this result holds for 
strength values in unrestricted as well as restricted networks. 

As a conclusion to our study of recurrence equations in the strength algebra, let us took at Ae 
relation between solutions of different equations. Define the relation <. between two functions/ and g 
as/ < g if and only if /(a) < «(a) for all a. The following theorem shows that this relation will then 
hold between the minimum solutions of their respective recurrence equations. This theorem 1 will prove 
valuable in comparing the steady state signals of different logical conductance networks. 
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Theorem 63. 

If monotonic functions /:!?" -* f* and g :3° -♦ 3° areordered/< g and a"* and .V" are the minimum 
solutions to the equations 

a = /(») 



respectively, then 



a - " < F*. 



Proof of Theorem 6.3: 

By induction on k / k (0) < $ k (0), and therefcie 



I — ♦ 00 C — f CO 



This theorem is also a special case of one given by Scott which states that the least fixed point 
operator is monotonic when applied to monotonic functions. 

6JS3 Solution Technique 

The following theorem shows that the minimum solutions of equations 6.15, 6.16, and 6.17 do 
indeed lead to the value of the steady state signal 
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Theorem 6.4 

For a matrix G G J 11 X n , and a vector bEJL n , define Gas G = xG. The unique minimum solution of 

the equation a = /(a), where 

/(«) = bVGom (6.19) 

is given by 

rf* = +«-" V -#*, 

where u™ 1 " and iF* 1 are the minimum solutions of the equations u = ffi&md i = ftfi), respectively, 
and the functions/j and^j are defined as: . > 

f x (u) = Moc*(i>1 T G'% r) (6.2Q) 

/ (d) = Moc*(LW T C*4 '). ' l -■■<« #•&> 



and 

r = G*-I*l. im 



The proof of mis theorem is given in Appendix IL It serves nwinly to confirm one's intuition but 
requires proving niany subtte pqii^ Prwia^^^ '*"" 

satisfies the recurrence relation a = /(a). It is then straightforward to prove that it must be the mini m um 

solution. 

If we let the vector b in equation 6.19 equal Eo x V jr. then we can apply Theorem 6.4 to find the 
steady state signal v. Alternatively, the reSufc of ttife theorem can tfc ixpres&d in a manne* more 
suggestive of a sequence of cohditidned relocations in mi s^ial algebra. Define the function 
kttl:JLxS**JL* 

kill (a, b) = +block(ral,b) V -6/oc*(LoJ,b). (6.23) 

In other words 
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kill(a,b) = f a.iai'2 
\X,lal< 



>b 
b. 



This function expresses our technique of killing any signal weaker than some strength value. The 
following corollary to Theorem 6.4 describes the method of conditioned relaxations. 



Corollary 6.4.1. 

For a matrix G G J x n , and a vector * £ U n , if Cis defined as G = xG, men the minimum solution of 
the equation «=/(«) where 

/(«) = * V Go* 
is given by 

k-*oo 
where 

/(«) = kUKfimliX (6JS) 

mi 

r = G*-I*l. 



The proof of Corollary 6.4.1 is given Appendix II. It proceeds by factoring the function / and 
showing that 

f\o> = */ x H^ v -fifm 

With mis, one can easily see that the sequence win converge toy*. 

The conditioned relaxation method succeeds where the straightforward relaxation method fails, 
because the function/ is monotonic over die domain { a | a < /*, and a - block (a, I f** I) }, whereas 
the function/ may not be. Furthermore,/' is closed over this domain. Thus, when successive values of c 
are formed by repeated applications of/, we will get a monotonic sequence converging to d*. This 
result is given as a corollary to the main theorem, because it does not generalize to switch-level networks 
containing transistors in the X state while the method given in the theorem does. 
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6.6 The Target State 

Returning to the problem of computing the target state y of a switch-level network, where the 
network may have conductance matrices in the sets { G } and { E }, we saw that the target state could be 
found by computing all possible steady states y (G, E) and men combining these possible states with the 
consistency operation U. Such an approach, however, would have exponential complexity if a large 
number of transistors were in state X. Instead, we will derive a method of directly computing y by 
generalizing the method of Theorem 6.4. 

To determine the target state y j of node /ij we need only find the range of values y j(G, E) can 
assume. If some setting of G and E can be round which gives 1 or X for y £G, E), and some other «ttiag 
can be found which gives or X, then y"| equals X; "Hit* the other hand, either of these two attempts 
fiufe, then y i equals the state found by me other attempt Define <G,E) as die vector of Steady state 
signals for the logical conductance meJw^rwimcondttetancemafiFfcerG and E. Then y j(G, E) equals 1 
or X if and only if rv^G, E)1 is greater tharifl. Similarly, y^CE) equals rtrX if and only if Lv^G, E)J 
is greater than 0. This suggests that the target state of a node can be found by performing two 
optimization processes. The first maximizes r vfiG, E)1 for aU possible G and E giving a result u * while 
me second maximizes Uj(G, E)J giving a result d * these values can be combined to give the target 
state: 



{ 



l.u^Oaod 4*^ = (6-26) 

*i = ^ 0, (F^O and u * = 
X, ebe. 



That is, the target state will equal a proper value (0 or 1) if and only if the corresponding optimization 
process succeeds (obtains a nonzero value), while the other fails. This can be expressed by the following 
equation: 

?i = «+if , i > U <-d*j>. (6.27) 
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At first this optimization approach might not seem to improve on the full enumeration technique. 
It seems to call for a separate set of optimizations for each node, with each involving the trial of a number 
of conductance matrices. Fortunately, these difficulties do not arise. Instead, the values I"v(G, E)l will 
be maximized for all nodes with one particular pair of matrices G and E, and a similar result holds for 
Lv^G, E)J. Furthermore, these values can be computed without ever finding the particular values of G 
and E which give rise to mem. Instead, the vectors a * and #"* can be computed directly by slightly 
modifying equations 620 and 6.21, as is shown in the following theorem. 



Theorem 6£. 

The target state y of a switch-level network is gwen by 

y = «Hf>U <-#•>, (6.28) 

where ii * and d * are die minimum solutions of the equations a = £ t («) and i - gM), respectively, 
and the functions gj and §q are defined as 

«!<■) = jfoc*(E—«r*i t ryif cr<-«,i) {%&) 

g^(d) = Woe*(E"--LxJ t k*J T G"" 'i,il (6J0) 

and 

r = G" , '*-(E - --IjcIT l/l). <6J1) 



The full proof of Theorem 6.5 is given in Appendix H. It involves showing that 

"" = {6HE}^ S <" 2 > 

where 0* equals the minimum solution of equation 629, and a"*(G, E) equals the minimum solution of 
the equation a = /^a) for 

/j(u) = block(E ' fxl t Tyl | G • a, G*»(E«I xl | !*■)). 
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That is u min (G, E) equals Vv(G, E)"l. We can see that u 01 " must be greater or equal to any u min (G, E) for 
any G in { G } and E in { E } as follows. For any such G and E, G raia <G<G ra " and E min <E<E max , and 
since block is monotonic in its first argument and antimonotonic in its second, this implies that/j < gj. 
Therefore Theorem 6.3 shows that u "' must be greater than or equal to u min (G, E). To complete the 
proof, we need only find a matrix G G { G } and a matrix E G { E } which give u rain (G, E) = u opt . The 
following matrices satisfy this requirement, although the proof is rather tedious: 

{ g min u opt. _ Q or u opt. _ 
g" 1 "^ u opt i >0andu op, j>0, 

r e^" TjcI = 
y \e max i j, TjCjIX). 

Observe that these matrices are defined in terms of the solutions they lead to. Our solution technique 
bypasses the search for the optimal settings of G and E and yields the optimal solution directly. By 
symmetry, one can see that a similar result holds for d 01 *. 

6.7 Explanation and Example of the Solution Method 

Theorem 6.5 describes an efficient technique for computing the target state of an arbitrary 
switch-level network. First we compute the vector r by applying the relaxation method to the equation 

r = E min • II jc 11 T II y\\ T G mto • r, (6.33) 

which gives the same result as equation 6.31. The elements of r equal the strengths of the steady state 
signals in the logical conductance network formed when only transistors in the 1 state are conducting, 
since the equation involves the matrices E rain and G min . For any allowable values of G and E, any signal on 
node «j with strength less than rj will be blocked regardless of the conductances of transistors in the X 
state and hence can be killed. We then compute the vectors u "' and d opt by applying the conditioned 
relaxation method to equations 6.29 and 6.30. These computations consider transistors in both the 1 and 
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X state to form logical conductances equal to their strengths, giving C^ and E"" in the first argument to 
block in bom equations, but any strength value on node «j less thai rj is Mocked The computations are 
performed separately for the 1 and signals, so that a signal (represented by a strength value) will not be 
killed if it could be the dominant signal for some set of transistor conductances. Signals of state X can be 
considered to have both state 1 and state -ft, and hence they enter into bo* computations. Once these 
strength values have been found, they can be combined to give fee target stales. 

This technique requires no enumeration over possible sets of transistor conductances whatsoever. 
This method cannot be formulated as a more intuitive method in the signal algebra, such as the method 
shown in Corollary 6.4.1, because our optimization technique does not correspond to the operation of any 
single logical conductance network. Furthermorej no simpler method has been obtained for restricted 
networks, because transistors in the X state in some ways, resemble weak transistors. That is, the signal 
paths in a restricted network cannot be analyzed independently when transistors in the X state are present 



fig. 63. Switch-Level Network Example 
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This method can be illustrated by the example shown in Figure 63, representing an mMOS Nand 
gate with both inputs equal to 1 connected to a pass transistor in state X. Assume that node n 3 has size icj 
and initial state 0. The recurrence equation for r can be written as: 

r l = n T (Y2* r 2) 
r 2 = ^2 T (Y2W 
r 3 = "1- 

The minimum solution of this set of equations is rj = r 2 = Y2> mA r 3 = K l- ^ recurrence equation 
for u op ' can be written as 



u x = block (yj T (Y2 1 u 2> t (Y2 1 "3). Y2> 
u 2 = block ty l.ui. Tj) 
u 3 = block ^i^, "I** 



The minimum solution of this set of equations is u^j = u "^ = ^3 ■* °» indicating- that regardless of 

the conductance formed by the pass transistor, no signal with state 1 or X can form on any nodes. Even 

though thepullup transistor provides a signal of strength +yj. our computation correctly recognizes that 

this signal will be blocked by the signal -y 2 - The recurrence equation for d * is 

d x = block ((y 2 i d 2 ) T (Y 2 i <*3>. Y 2 ) 
d 2 = block (y 2 t (Y 2 i dp. Y2) 
d 3 s= block (jcj T ty-lty.. icj). 

The minimum solution of this set of equations is d "^ = d ^ = d ^ = Y2- Thus ' since ^^ values are 
all nonzero, while the values of u 01 * are all 0, all three nodes have a target state 0. 

If the same network has initial state 1 for /13, we would find that u ^ = k^ while all other 
elements of u 01 * and d * have the same values as before. This gives target states y j = y 2 = 0, and 
y 3 = X, indicating that the unknown conductance of the pass transistor creates an ambiguity in the target 
state of node /13. 
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6.8 Properties of the Target State 

Now that we have a mathematical description of the target state, several useful properties can be 
demonstrated. 

6.8.1 Monotonicity 

Our partial ordering of signal states ranks states according to how well defined they are. The state 
_L is underfined, i.e. it represents an absence of information. This state will never appear on a node in a 
switch-level network, because all nodes store information dynamically and hence can never be devoid of 
state. The states and 1 are well defined, Le. they represent a consistent degree of information. The state 
X is overdefined, i.e. it represents conflicting information. The following theorem shows that the target 
state function is monotonic for this ordering. This indicates that setting some node or transistor to X can 
only lead to target states fin- some nodes equal to X which would otherwise equal Boolean values. 



Theorem 6.6. 
Ifx<x',y<y',z<z'then 

target (x, j, z) < target (x'.y'.x'). 



Proof of Theorem 6.6: 

This theorem can be proved by comparing the derivation of die target state for initial values x, y, z 
(which will be shown with unprimed values), with the derivation for initial values x\ y', z' (which will be 
shown with primed values.) Compare the function 

gj(u) = Wodfc(E"« • m T r>l T <?""•■. r) 

with the function 

gj'(u) = block (E-* • fx'i r ryiTG— '•■,0k 
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where 

r = G nU "*«(E mi,, «*xl T *3>*\ 



and 



r' = G mta '*«(E n,to ' • H x H t iy«). 



One can see from equations 6.1, 6.2, 6.3, and 6.4, that if z < i' then 

G mto > Q-to' 

E"*" < E™"' 

/-mux ^ /«■' 

Since the strengths of the initial stimulus signals are determined only by the node types and sizes, 
II x II = II x II, and II y II = II y ». Therefore r > r'. Furthermore, if x < x\ then x < x and therefore 
Txl < I"*'! Similarly, iyi < !>'!. Therefore g t < ^j' and by Theorem 6.3, u * < u " 4 '. By similar 
reasoning, one can see that g Q < g^ and therefore d * < *"*'. Observe mat the function of b whose 
value is <+b> is monotonic, and therefore <+u opt > < ^h-u^'^ and by similar reasoning 
<-&°*> < <-d°*>. Finally, by the monotonicity of U 

target (x,y,z) = ^u^* U <-#*> < <+u° s *> U <-&<*> = target (x'.y'.z). I 

The monotonicity property then extends to the functions step x and phase. 



Corollary 6.6.1. 

Ifx < x' and y<y' then 



stepJy) < step_<y'). 
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Proof of Corollary 6.6.1: 

One can see from the table in Section 2.5 that trans (x, y) < trans (x, y'). Therefore 

step(y) = target (x, y, trans (x,y)) < target (x',y\ trans (x',y')) = step ¥ '(y'). 



Corollary 6.6.1. 

If x < x' and y < y' then 

phase(x,y) < phase{x,y). 



Proof of Corollary 6.6.1: 

By Corollary 6.6.1 and induction on k, 

step x \y) < step x - k (y'). 

Therefore 

phase(x,y) = lim step x Hy) < i Um step- k (y') = phase{x' ,y*). I 

k— +oo A k— »oo * 

These results show that the presence of an X value on a node can only lead to new network states which 

have some nodes set to X which would otherwise be set to Boolean values. 

6.8.2 Stability of the Target State 

The target state is claimed to be the set of states which the normal nodes would eventually reach if 
the input nodes and transistor states were held in states x and z, and the normal nodes were initialized to 
state y. To really prove this, we must show that once the network reaches the target state, it will stay there 
until some input node or transistor changes state. While this stability can readily be seen for the steady 
state of a logical conductance network, it is less clear for the target state of a switch-level network. The 
following theorem eliminates any such doubts. 
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Theorem 6.7. 

If y = target (x, y, z), then y = target (x, y, z). 



The proof of this theorem is also given in Appendix II. It involves showing that the terms r/1 and 
LyJ in equations 6.29 and 6.30 can be replaced with terms fjl and LjU where j? is the vector of signals 
with fth element having state equal to the target state y j and strength equal to the node size capj. 

6.9 Summary 

Our entire development of the switch-level model so far can be summarized by three equations 

r = Ef'l.xl.^ljrl t.Gfv ( 6 - 33 > 

u = blockiE™ -Txl t V|G"'M) (6-29) 

d = block (E— -LxJ f L.>d > G" 4 " '4 r). (6.30) 

By finding the minimum solution of the first equation and then using mis vahie in computing the 
minimum solutions of the other two, we obtain two vectors of strength values tf* and d * from which the 
target state for each node can be computed as 



{ 



1, u^Oandd^O (6.26) 

Vi ~ \ &, €* x > ©and u^sO 
X, else. 



Consider how far we have progressed rrom me electrical eacuit-oriented view of the switch-level model 
provided by the original definition of the target state in Chapter 2. This new method involves only simple 
operations in a discrete algebra, and the equations can be solved by a straightforward iterative method. 
Furthermore, it finds whether nodes are sensitive to the unknown conductances formed by transistors in 
the X state (and hence should have target state X) without enumerating over possible combinations of 
transistor conductances. Thus, it can be implemented by a very efficient computer algorithm as is shown 
in the next chapter. 
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7. Simulation Algorithms 
7.1 Introduction 

Theorem 6.5 defines a straightforward method for computing the target state function, and this 
method can be implemented by an efficient algorithm to serve as the basis of a switch-level simulator. By 
exploiting the locality of both the interconnections and the activities in the network, the program can 
achieve a performance comparable to logic gate simulators. First a unit delay simulation algorithm is 
presented which provides the same functionality as the program MOSSIM [9] for designs which can be 
described in the MOSSIM network model. Next, it is shown that a slight modification yields a simulator 
with a timing model similar to Terman's program [5], although the functionality of the two algorithms 
differ significantly. The new algorithm differs greatly in its style from both of these previous algorithms, 
largely because it is based on solving equations in a well-defined mathematical domain rather than on the 
intuitive ideas of the simulator designers. These simulation algorithms are compared and contrasted 
toward the end of the chapter. Some performance data from MOSSIM are presented to demonstrate the 
performance characteristics of switch-level simulation and how it compares to logic gate simulation. All 
algorithms are presented as "Pidgin Algol" programs as defined in Aho, Hopcroft, and Ullman [1]. 

Before delving into the details of the simulation algorithm, let us consider its intended mode of use. 
Suppose the design to be simulated will operate as a synchronous circuit with a conservative clocking 
scheme. That is, some external set of clock signals will be provided through input nodes which control 
the sequential operation of the circuit such that as long as these clocks run slowly enough, no timing 
errors can occur. Each clock cycle can be subdivided into a set of simulation phases (called "epochs" in 
Mead and Conway [37]) where during each phase, all clock and data inputs remain constant For 
example, a two-phase, nonoverlapping clock contains four such simulation phases: 
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Phase 

12 3 4 

Phil 10 

Phi2 1 

During each phase, the circuit has sufficient time to stabilize. 

To model the functionality of such a circuit, a simulator can simply compute the state in which the 
network would settle for each phase of each clock cycle, setting the clock and data inputs to new values 
between the phases. The function phase, defined in Chapter 2 as phase(x, y) = , '"" stepJy) serves 
this purpose. To the user, this technique provides the effect of a unit delay timing model in which 
transistors switch one time unit (i.e. on the next computation of stepj after their gate nodes change state. 
Such a technique provides only limited information about the speed of die actual circuit, but gives an 
indication of the function computed. The characteristics of this and other timing models are discussed in 
Chapter 8. 

This technique has been applied to simulating self-timed systems [36] as well, in which activities 
may occur independentiy and asynchronously. Each phase then corresponds to a particular setting of the 
input data and control signals, and it is assumed that the circuit will settle before the inputs are changed 
Although actual circuits may not obey these assumptions, almost all can be modeled as if they did. 

7.2 Complexity Model 

In a switch-level network, each node could be connected to every other node by any number of 
transistors, giving an unbounded number of transistors relative to the number of nodes. In practice, 
however, the number of transistors grows only linearly with the number of nodes due to the limited 
connectivity allowed by a two-dimensional integrated circuit chip and to electrical and functional 
considerations. To evaluate simulation algorithms, we should have a model of the complexity of 
networks which more closely matches actual circuits. The following set of assumptions has been observed 
to hold for a variety of designs, although it has not been subjected to a rigorous study. 
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Define the connectivity set of a node as the set of transistors for#nich the node serves as the source 
or drain connection and the connectivity degree as the size of this sit The fanout setaf a node is defined 
as the set of transistors for which the node serves as the gate connection, and the fanout degree is defined 
as the size of this set An informal study of a variety of designs has shown that almost all normal nodes 
have connectivity degree less than 5. Exceptions include large busses and the output nodes of large Nor 
gates. Input nodes, especially VDD and GND, however, may have a high degree of connectivity. A 
similar statistic holds for fanout degree with the exception of nodes providing major control signals such 
as clocks and reset or enabling commands. 

We will assume the network may contain 0(1) (i.e. a constant number) of input nodes each with 
0(n) fanout and connectivity degree, where n is the number of normal nodes. An 0(1) subset of the 
normal nodes may also each have 0(n) fanout and connectivity degree, but the remaining normal nodes 
must each have 0(1) fanout and connectivity degree. From either the fanout or the connectivity 
assumptions, one can see that die network can contain only 0(n) transistors. 

The sparseness of interconnections in a logic design leads to a localization of the activities. When a 
node changes state, generally only a small number of nodes will be directly affected. Furthermore, in 
most synchronous designs, each logic element will be activated only a small number of times during each 
clock cycle. That is, during a single simulation phase, information will only propagate ftooa tfie outputs 
of one set of storage elements through some combinational logic to the inputs of other (or perhaps the 
same) storage elements. Even allowing for a small number of dynamic hazards (transient pulses caused 
by unequal path delays), each node win change state only 0(1) times during each phase. In fact, 
experience has shown that often significantly fewer state changes occur. For example, in a random access 



1. Structures involving many transistors of the same strength and type connected in parallel, such as 
large Nor gates could be simulated more efficiently if mey were modeled by special "multi-transistors" 
which have multiple pte nodes, any of which can activate the switch. A simple count of the number of 
gate nodes in the X and 1 state would indicate the state of such an element 
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memory, only a small percentage of the nodes change state during each clock cycle. While this example 
represents an extreme case, most networks contain only a small number of active elements at any given 
time. 

7.3 Sparse and Incremental Equations 

Our complexity model shows that the connectivity in the network is very sparse, and that changes in 
the network state occur only incrementally, i.e. a small number of nodes at a time. A well-designed 
simulation algorithm can exploit both of these reductions in complexity and thereby achieve considerably 
better performance than would a naive implementation of the matrix equations. These techniques will be 
demonstrated by developing algorithms for solving sparse and incremental equations in the strength 
algebra. 

7.3.1 Sparse Equations 

Suppose we wish to find the minimum solution of the equation a = /(a) where f-.f 1 -> tf* can be 
expressed as 

W>\ = b 4 T I /y(aj) (7.D 

'j c i 

The set P-. is called the adjacency set of node n y and in our application will equal the set of normal nodes 
connected to the node by transistors in the 1 or X state. We will assume that all connections are 
bidirectional, i.e. «j G P-. if and only if n: G Py Furthermore the function / is assumed to be both 
monotonic and passive. The general definition of passive functions is described in Appendix II, but for 
the function/ above implies that for all indices i and j and all strength valujs s,f-Js) < s. 
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As an example the function /j defined in equation 630 as 

^(a) = block (r*1 T G«a, r) (6.20) 

can be expressed in this form with bj = block (rip, q) and/y(aj) = block (gy 1 3j, q). One can see that 
f r is both monotonk and passive. This example also motivates the term "passive"! "which corresponds to 
the property that a signal coupled through a logical conductance can never be increased in strength. This 
example also demonstrates how the locality of interconnections in the network lead to sparse matrix 
equations. Element y of matrix G can only be greater than if a transistor in the 1 or X state connects 
nodes n { and /ij. By our assumptions about network a>nnecuvity, j .Pj j can be 0(n) for only 0(1) values of 
i, and must be 0(1) otherwise. Therefore 2| /* [ must be 0(n). 

The following program solves this equation where S denotes some data structure such as a stack in 
which elements can be inserted (push) and removed (pop) in unit time. The order in which elements are 
removed is unimportant 
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proeedureSOLVEK/): 
begin 

S«-0; 

for i «- 1 until n do 
begin 

push(S, /ij); 
donejt-j&fae 

end; 
whileS 7^0 do 
begin 

n^ «- pop(S); 

if done: = false then 



beght 



done: «— fra«; 
for each «j G P; do 
if/ ij (a j )>a i then 
begin 

»i+- fifty 
doae^*- false, 

push(S, «j) 

end 

end 

end; 

retura(a) 

end 

The procedure SOLVE1 resembles the relaxation method outlined in Chapter 4, except that it tries 
to minimize the amount of computation by computing the effects of a node value on adjacent nodes only 
when the value changes. It starts by setting all nodes to the initial values given by b and placing the nodes 
in the list S. At any time, any node * in the list for which donej equals fake has had a new value assigned 
to a:, the effect of which has not been propagated to neighboring fiode& The procedure performs a series 
of relaxations, each of which starts by selecting a node flj from S such that donej equals false. The effect 
of the value a: on each neighboring node in P y Le. /-(ap is computed, and if this exceeds the previous 
value on the neighboring node, the neighbor is updated and placed in the list with donej set to false. This 
may lead to duplications in the list S, because n- x may already be on it. The flag donej, however, provides 
a means of checking whether the effects of the node value on adjacent node values have already been 
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computed. This technique eliminates the need to test a node for membership in S before adding it to S. 1 
These relaxations continue until all nodes are consistent with one another, i.e. any further relaxations 
would have no effect 

The correctness of the procedure SOLVE1 relies only on the monotonicity of /. This can be proved 
by inductive assertion on the number of relaxation steps (i.e. executions of the while loop body.) Let a k 
equal the value of the array a after k steps. It is claimed that for any k, b <a k < a""" and for any j such 
that donej equals true, /y^) < aj for all i in P y This clearly holds at the start, because 
b = a =/(0) < a™" and donej equals false for all j. Now suppose mis assertion holds after relaxation 
step k. A relaxation step involves propagating values according ttHae Set of functions /y, and therefore 
a k < a k+ * < a k T/(a k ), which gives 

•» < a k < a k+1 < a k T/(Jr k ) < a""" T /(«"") = a"*". 

Furthermore, the second part of the assertion clearly holds, because the program sets done: to false any 
time aj has been changed and only sets it to true once the effects of the new value have been computed. 
By induction the assertion must hold when the procedure terminates with done: equal to true for all j. 
This implies that for the final value of the vector a, aj > fJaJ for all i and j. We also know that aj > bj 
for all i, and hence a >/(a). By the monotonicity of/ and induction on k, a >/*(a) for all k, and 
therefore 

* ^ v *"«/**■> ^ w *L/*to = ■" ,, • 

ic-"»go k— *oo 

Combining this inequality with the inequality of the induction assertion gives a = a - " when SOLVE1 
terminates. 



1. One could test the flag donej before inserting a node in S to avoid duplication in S, but the method 
given here leads more naturally to our further developments. 
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To analyze the complexity of this procedure, observe that a nodefc pashedon S only -when its value 
is increased. Thus each node n { can only be pushed a maximum of \T | times, each time causing at most 
one relaxation step requiring | P { | operations. Therefore the algorithm has coi^exity less than or equal 
to 0(21* 1-1^1) = 0(|y|.2|P i |) = 0(|J|.n). 

The following example shows that this degree of complexity can actually be realized. Define b as a 
vector of decreasing strength values followed by ffs, Le. 



{ 



b i - t *p+q+l-v P*l<*<?*9 
0, p+q+ l<Zi<p. 



Let P l = { n 2 }, P n = { n^ }, and P { = { n^. « i+1 } for 1 <i< n. Define /y as the identity function 
for «: G P. This example corresponds to a linear chain of normal nodes with the nodes having initial 
signals of decreasing strength. Suppose that S is implemented as a stack and that nodes are selected from 
each set P x in order of their subscripts. SOLVE1 will first set all nodes itj such that i > p+ q to * v and 
then it will set all nodes Hj such that i > p+ q-1 ttt«$ andso on thfo*gh all possible strength values until 
finally all nodes are set to y . Thus, the worst case compteiaty of a SOLVE1 equals 0(| 1 |-n). 

For most MOS networks we can assume that t is a very small set For example, the network model 
of MOSSIM can be implemented with J = { 0, K V y v y 2 }• Therefore SOLVE1 provides a linear and 
hence optimal solution. Nonetheless, our worst case example shows that SOLVE1 can waste much effort 
in propagating values which will only be overridden later. A slight refinement, however, leads to an 
algorithm with complexity 0(n). It replaces the list S with and array ef lists B, with one list for each 
possible strength value. 
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procedure SOLVE2</): 
begin 

for each se^dol^sj 4-0; 
for i *- 1 until n do 
Dcgn 

ai^-bj; 
pushCBIaj], «j); 
done:*— false 
end; 
PROPAGATE(B, a,/); 
return(a) 
end 

In this program node n- x is inserted into the stack corresponding to the initial value bj. The procedure 

PROPAGATE, defined below, then spreads these values through toe network. 



procedure PRO 


PAGATE 


(B,a,/): 


begin 






s*"V 






repeat 






begin 






while B{sJ ^ 


0do 




luinia 






V 


-pop(BIsD; 




if done j =^jiwtben 






begin 






done| *- true, 






pnitj 4- false,] 






for each /^ G /j do 






if/ y (a i )>a i t-« 






DcgBi 






<4 *~fd»dl 






auJ 






ead 




cad; 




s«- 


pndto 




end 






until s = 







end 

The function />r«/ in this procedure is die predecessor function for 3. That is, pred(Q) = 0, and for s > 0, 
/va/(s) equals the greatest element of 3 less than s. This function is used to enumerate the strength values 
in descending order. The line enclosed in square brackets is required for our next extension of the 
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algorithm. 

The procedure SOLVE2 operates much like SOLVE1, except that a node is selected for relaxation 
only if it has maximal strength of all nodes for which donej equals false. It does this by keeping the nodes 
on separate lists according to their strengths and going through these lists in decreasing order. It relies on 
the passiveness of/ to assure that if a node of maximal strength is selected for relaxation, no nodes will be 
set to greater strength during the remainder of the computation. Therefore each node is selected for 
relaxation only once, and hence the complexity of SOLVE2 equals 0(| P { |) = O(n). It is unclear whether 
SOLVE2 will actually achieve a better performance than SOLVE1, because the cost of implementing the 
array of stacks might exceed the gain in efficiency. This depends on the details of the programming 
language as well as on the networks to be simulated. Nonetheless, we shall pursue the algorithm for 
SOLVE2 for further development 

7.3.2 Incremental Equations 

As an extension of this technique, suppose that we have computed the minimum solution a m,n of the 
equation a = /(a) with/ defined by equation 7.1 and now wish to find the minimum solution a min ' of the 
equation a = /'(a), where 



[Aa)]i = b'j T n .J^/W 

and /' obeys the same restrictions as /. Furthermore, assume that b'j differs from bj for only a small 
number of values of i, that /'y differs from /y for only a small number of values of i and j, and that 
P- ^ P 1 - for only a small number of values of i. 

For the example of the function /j given earlier, the differences between bj and b'j reflect changing 
input node states or changing connections to input nodes. The differences between /jj and /'y and 
between P-. and P"-. reflect changing connections between normal nodes. 
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The differences between the functions/ and/' can be regarded as perturbations of the function/. 

That is node n { is perturbed if b; ^ b' v P { ^ ^,/y ^/.., <*/.. ^/.j. Such a perturbation can only 

affect nodes in the vicinity of node n- v where «j is defined to be in the vicinity of n { if there exists some 

path from «j to n { in the graph with edges defined by the connectivity sets / y k . Thus, the new equation 

can be solved "incrementally", by only computing new values for those nodes in the vicinity of a 

perturbed node. 

procedure INCR_SOLVE(/",/, a): 
begin 

E«-0; 

for each i such that b { ^ b\ or P { ^ f [ or f.. jLf or /.. j&f.. do 
begin j j j j 

E-EUlnj}; 

donej«— false 

end; 

for each sG *do B[s]<- 0; 
for each n^EE such that donej = false do 
begin 

INITIALIZER, b', a, B); 

PROPAGATE(B, a,/) 
end; 

return(a) 
end 

The procedure INITIALIZE sets all nodes in the vicinity of ^ to the initial values given by b', using a 
depth-first search technique for finding all connected nodes in a graph as is described in Aho, Hopcroft 
and Ullman [1]. The flag initj is used to indicate whether the node has already been initialized. It is 
assumed to equal false at the start of the program and is also set to false by the procedure PROPAGATE 
once the relaxations begin. The procedure also places each initialized node into the appropriate list in the 
array B. 
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procedure INITIALIZER, b', a, B): 
if initj = false 
then 
begin 

initj «— true, 

push(B[aj], /ij); 
for each n-. £ P- y do 

INITIALIZE(nj, b', a, B) 
end 

The procedure PROPAGATE has already been given. It will take the initial values set by INITIALIZE 

and spread them through the network in the vicinity of the perturbation. Observe that the flag donej is 

used for two purposes. It is used by the procedure INCR.SOLVE to avoid redundant computations 

when two perturbed nodes are in the same vicinity. It is also used by the procedure PROPAGATE to 

take care of the case where a node has been moved to a new list but has not been deleted from the old. 

The complexity of the procedure INCR_SOLVE is proportional to the number of nodes in the 

vicinity of a perturbed node. This can range from 0(1) to 0(n). In any case, the procedure is close to 

optimal, because it only looks at nodes in the vicinity of perturbations. A truly optimal algorithm would 

only look at a node if its value will change, but this is hard to achieve. 

7.4 Unit Delay Simulation Algorithm 

Theorem 6.5 shows that the target state can be computed by solving equations for r, u op \ and d ** in 
the strength algebra. We will now drop the superscripts on u and d. As the phase simulation progresses, 
these values change only incrementally. Thus the techniques developed in the previous section lead 
directly to an algorithm for a switch-level simulator. In the following programs it is assumed that the 
network structure and state information as well as the vectors r, u and d are available as global variables 
and need not be passed as arguments. Only an outline of the algorithm will be given here, because it 
involves many details and requires further development of the proper set of data structures. 
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The program PHASE shown below computes the function 

phase(x,y) = lim stepjHj). (2.3) 

k— +00 * 

It takes advantage of the fact that each computation of step x involves only incremental changes to the 

network state. The procedure is given a list of node-state JMirs .as. anar^ument. Each element of this list 

is of the form <jj, x>, indicating a new setting for art input node, or <n-, y>, indicating a new setting for a 

normal node. The variable "newval" represents some element of this list In general, only input nodes 

will be changed, because in an actual circuit only these nodes are accessible externally. 

procedure PHASE(A): 
begin 

E«-0; 

for each newval E A do 
begin 

S£T_NODE(newval); 
E 4- E U PERTURBJ^QDECnewval); 
SETjmANS«?ltMlS(aewv4; 
E i- E U PERTURB^TRANSISTORS(newvaO 
end^ 

while E^0 do 
E«-STEP(E) 
end 

The procedure SET_NODE updates the node state, and the procedure PERTURJkWDE places those 
normal nodes perturbed by mis change in the list E as well as sets die flags donej to^sefor these nodes. 
That is, a changing input node will perturb all normal nodes connected by conducting transistors, white a 
changing normal node wfll perturb only itself. The procedure SETJTRANSBTORS updates the state of 
every transistor in the fanout set or the node, while the procedure PERTORBLTRANSISTORS places 
each normal node for which a transistor in its connectivity set has changed state in the Bst E and sets die 
flags donej to false for these nodes. The procedure STEPsimulates the effects of the perturbations on die 
nodes in E which then creates a new set of perturbations to be simulated. This process continues until a 
stable state is reached, i.e. no perturbations remain. 
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The procedure STEP is shown below. 

procedure STEP(E): 
begin 

A<-0; 

for each «j G E such that donej = false do 

A «- A U UPDATE^); 

E*-0; 

for each newval £ A do 
begin 

SET_TRANSISTORS(newval); 
E^-EU PERTURB_TRANSISTORS(newval) 
end; 
return(E) 
end 

The procedure STEP selects a node from the list of perturbed nodes and calls the procedure UPDATE to 
compute the new states of all nodes in the vicinity of the selected node. Those nodes which change state 
during this process are accumulated in a list A. Once the effects of all perturbations have been simulated, 
the transistors in the fanout sets of nodes in A are set to their new states. This will cause new 
perturbations which are accumulated in the new list E. Observe that the procedure PERTURB_NODE 
need not be called, because by Theorem 6.7. the target state will remain stable unless either an input node 
or transistor changes state. By changing the transistor states only after all perturbations have been 
simulated, the procedure STEP creates the effect of all node states changing simultaneously and then all 
transistor states changing simultaneously. This implements a timing model in which transistors switch 
one time unit after their gate nodes change state. 

The procedure UPDATE is shown below. It applies the technique shown in the procedure 
INCR_SOLVE to solve the equations 

r = E mln • II x II T M' T G mia • r (6.33) 

u = block (ET^'Txl T r/1 t G max • u, r) (6.29) 

d = block (E max 'LjcJ T L>d T G max »d,r). (6.30) 

From these the value of the target state is computed for each node as 
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{1, Uj>Oanddp= 
0, d i >Oandu j = 
X, dsc 



= (6.26) 





Each node which changes state is placed in a list, along with its new vatae. 

procedure UPDATE^) 
begin 

for each s € B(sJ do B|sJ *- 0; 

INITIALIZE-Rdij, B); 

PROPAGATE_R(B); 

INITIALIZE-U^, B); 
PROPAGATE_U(B); 

INITIAIJZEJXnj, B); 
PROPAGATEJXB); 

A <- UPDATEJSTATE(/ii); 

retum(A) 
end 

The remaining procedures will not be given. 

The speed of the procedure UPDATE could be improved for the case where no paths of conducting 

transistors from n^ to input nodes or other normal nodes contain transistors in the X state. Suppose 

furthermore that all paths from nj to other normal nodes contain only transistors of strength ?-. Then 

Corollary 6.1.1 shows that the steady state signal for every node connected by some path to it can be 

found by combining the initial signals on these nodes using the operation V, and the state of this signal 

equals the target state of all of these nodes. The procedure INITIALIZEJR could check whether this 

particular condition exists, and if so a procedure could be called to perform this computation, and men 

the nodes could be set to the state of tbis signal. This cempujatipn involves considerably less effort than 

computing r, u, and d. Considering that in most simulations the X state arises only rarely once the X's 

present at the start of the simulation have been forced away, such an optimization could provide a 

substantial speed gain. For the case where some normal nodes are connected by transistors with strength 
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less than y (i.e. an unrestricted network) we could employ the method of Corollary 6.4.1 to find the 
steady state signals and consequently the target states. This method, however, may not provide a 
significant savings over the original method. 

Our network complexity gives a rather weak bound on the complexity of PHASE. If we assume 
that each node changes state only 0(1) times, then the total number of perturbations to the network will 
be proportional to the sum of the fanout degrees of all nodes, which is 0(n). However, a perturbation 
may require 0(n) operation to simulate its effects, although such cases are rare. This gives a total 
complexity of 0(n 2 ). Such complexity is achieved only by highly contrived examples, however. 
Experience has shown that typical networks require at most 0(n) operations per phase. 

7.5 Pseudo Unit Delay Simulation Algorithm 

The algorithm presented in the previous section carefully holds all transistors Fixed until all 
perturbations have been simulated, thereby creating the effect of all nodes changing state simultaneously 
and then all transistors switching simultaneously one time unit later. If we instead switch transistors 
immediately after their gate nodes change state, we obtain an algorithm with many characteristics of a 
unit delay simulator but in which all events are ordered. The characteristics of this timing model are 
discussed further in Chapter 8. 
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procedure PHASE2(A): 
begin 
E*-0; 

for each newval £ A do 
begin 

SET_NODE(newval); 
E^EU PERTURB_NODE(newval); 
SET_TRANSISTORS(newval); 
Ef-EU PERTURB_TRANSISTORS(newval) 
end; 

while E ^ do 
begin 

«i «— dequeue(E); 
if donej = false then 
begin 

A «- UPDAT^/ij); 

for each newval £ A do 
begin 

SET_TRANSISTORS(newval); 
E^EU PERTURB_TRANSISTORS(newval) 
end 
end 
end 
end 

In this procedure, elements are removed from the list E in first-in, first-out order so that the effects of 
simultaneous perturbations will be simulated before any subsequent effects are simulated. As a result, the 
algorithm provides a similar timing model to the unit delay algorithm, even though the list E evolves 
continuously during the simulation phase rather than being repeatedly filled and emptied. 

7.6 Comparison to Other Switch-Level Simulators 

Both MOSSIM [9] and the simulator developed by Terman [5] are designed along similar lines to 
the ones described here. A comparison between these three simulators serves to highlight some issues in 
simulator design. 
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7.6.1 MOSSIM 

For networks which can be described in the MOSSIM network model, the unit delay simulation 
algorithm presented here will produce the same results as MOSSIM. The two programs, however, differ 
greatly in their internal structure. 

MOSSIM precedes the simulation by a relatively complex network analysis, primarily to partition 
the network in a way which allows selective updating. At the start of this analysis each input node is 
replicated to give a separate copy for each transistor in its connectivity set. Then the network is 
partitioned into transistor groups with each group corresponding to a connected component in the 
undirected graph with a vertex for each node and an edge between each pair of vertices corresponding to 
the source and drain nodes of a transistor. This partitioning divides the network into components which 
interact only through fanout connections, i.e. from a node in one group to a transistor gate in another. 



Fig. 7.1. A MOSSIM Network Partitioned into Transistor Groups 




pullup node 

normal node 
input node 
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Since we assume the gate node is electrically isolated from the source and drain, such a connection & 
purely unidirectional and depends only on the state of the node signal and not on its strength. Thus each 
transistor group can be yiewed as a logic block with inputs, outputs, and internal state and which 
communicates with other blocks only with logic values 0, 1, and X. An example of a MOSSIM network 
partitioned into transistor groups is shown in Figure 7.1. A similar partitioning method has been used by 
researchers at the Nippon Electric Company in an analog simulator to achieve a localization of activities 
[42]. Although our new algorithm does not require this partitioning, such a technique provides a way of 
introducing some degree of structure into die network. This structuring has several applications that will 
be described shortly. 

During the simulation a transistor group need only be simulated when one of its inputs has 
changed. This tends to restrict the simulation to me active portions of the network, thereby achieving 
some of the effect of the incremental updating technique used in the new algorithm. MOSSIM, however, 
can only take advantage of the static locality in the network. U. thai .which can be detected without 
considering transistor states. The new algwitlim ato tiuraadwrntag* o/^ 
source and drain nodes of a transistor in the state are also l considered electrically isolated. Thus, only 
nodes connected by paths of conducting transistors to a perturbed node heed be updated. This added 
selectively should yield some gain in speed. 

MOSSIM also uses a significantly different method jrf updating a $et of nodes following a 
perturbation. It exploits the feet that in a restricted network a set 0f normal nodes connected by 
transistors in the 1 state will attain the same target state as was shown in Corollary 6.1.1. MOSSIM 
simulates a transistor group by partitioning the normal nodes in a gronp into equivalence classes and 
computing the steady state signal for each class, ignoring all transistors in A* X state. Then the effects of 
transistors in the X state are simulated by first forming a "supergraph" widi a vertex for each equivalence 
class and an edge between two vertices if a transistor in the X state has its source node in one class and its 
drain node in the other. This supergraph is inspected to see which classes should be set to X because of 
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possible connections to classes with different state and greater or equal strength. 

The MOSSIM algorithm has several drawbacks which motivated the new approach. First, the 
initial network analysis, as well as the computation of equivalence classes, supergraphs, and so on involves 
a great deal of dynamic storage allocation. TMs requires much computational effort and cannot be 
programmed easily in languages which lack automatic heap storage management The new algorithm, in 
contrast, utilizes only recursive procedure calls and data structures such as sets and a small array of stacks. 
Furthermore, the MOSSIM algorithm cannot be extended to unrestricted networks easily. Although a set 
of nodes connected by transistors with strength y p and state 1 in an unrestricted network will form an 
equivalence class with the same target state, the computation of this state becomes much more difficult 
While almost all MOS designs can be described as restricted networks, the greater generality of the new 
algorithm gives added flexibility. 

7.6.2 Terman's Simulator 

The algorithm used in Terman's simulator provides a timing model similar to the pseudo unit delay 
algorithm presented here. These two algorithms differ significantiy in their functionality and internal 

structure, however. 

Terman's program deals only with restricted networks. For nodes connected by paths of transistors 
in the 1 state, it combines the initial signals on the nodes with an operation similar to the operation V, 
just as was suggested to improve the procedure UPDATE. As mentioned in Section 2.9, however, chafge 
sharing is simulated by real-valued capacitances. For nodes conaeriedby transistors in the X state it 
attempts to encode information about the network condition into additional "states" and propagate this 
information much as it does the other states (which can be likened to signals.) However, the small 
number of additional states provides insufficient detail about the network condition/and this forces a 
rather inaccurate simulation. Furthermore, simply adding more states would not correct this problem, 
because it seems as if an accurate algorithm requires two passes over the set of nodes: the first to perform 
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a pre-conditioning step and the second to compute the new node sates. Ik functionality of Terman's 
algorithm can be expressed in the signal algebra (except ibr charge sharing^ 



y 



■ I 



1 *k vj^.i^^^r). 



That is, a node is set to X unless it has the same steady state signal when all transistors in the X state are 
fuDy conducting as when they are nonconducting. Furthermore, when charged nodes in different states 
are connected by transistors in the X state, no attempt is made to take their relative capacitances into 
account Instead these nodes are set to X. As was shown in Section 19, one has little choice in this case, 
because real-vahied capacitances and transistor in die X state (apparently) cannot be dealt with 
consistendy. It can be shown that Termans algorithm is more conservative than the one given here, 
except when charge sharing is involved. That is, whenever it sets a node to or i ours sets the node to 
the same value, but in some cases it may set a node to X when ours sets the node to or 1. To see this, 
recall that for a restricted network, both G"* 1 and G"" must be elements of { 0, y.} nXn , and therefore 
any matrix in the set { G } must also obey mis property. Let E equal any matrix in the set { E } and G 
equal any matrix in the set { G }. Equation 5.17 shows that o is monotonic for arguments representing 
conductance matrices. Therefore for any m G JL n 

xE^ox V jr V xG-^oa < xEoxVjrV xGo« < xB-ojr V/V xG""o* 

ami furtheraKire, these funcuoiis are oioaotonkmc AUieoremsiiml^toTrieoremU thenshowsaiat 

KG-M"*) < i(G,£) S «^-.P^ 

If v i (G" ta ,E- fc )= vj(G--,E-»), then 7j(G,E) = j£G" t P*) for any GG{G} and E<E{E}. 
Therefore y t = y ^G - ', E"*). Thus, for the cases in which Terman's simulator sets a node to V or i, 
ours sets ft to the same value, except when charge sharing is involved. 
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Fig. 7.2. Inaccuracies in Terman's Simulator 



Initial Correct Simulation 
Value Result Result 



IK* 



II- * 



X y x X X 

J_ y 2 1 1 X 

(^=1.0 ^=4.0 



As the examples in Figure 7.2 show, however, Terman's algorithm may set a node to X when it 
clearly should set it to or 1. These examples are shown in the MOSSIM network model in which the 
pullup resistor corresponds to a d-type transistor of strength -y^ while all other transistors are n-type 
transistors of strength y 2 . In me first example the node is being driven to 1 by a transistor of strength y± 
in the 1 state and a transistor of strength y 2 in the X state. The node will have a different steady state 
signal if the second transistor is conducting or nonconducting, and hence Terman's program will set it to 
X. The state of the steady state signal will be 1 in either case, and hence the node should be set to 1. The 
second example shows a similar result when a normal node is initially charged to 1, and then connected 
by a transistor in the X state to VDD. In the third example, if the transistor had zero conductance, node 
« 2 would stay charged at 1, while if it has nonzero conductance, the two nodes will share charge but n^ 
will remain above the positive logic threshold. Therefore « 2 has a target state of 1, but Terman's program 
sets it to X. In many other cases, however, such as in simulating logic gates with some inputs in the X 
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Fig. 7.1 Inaccuracies in Terman's Simulator 
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As aw examples in Figure 7.2 show, however, TettahV algorii^ a note to X when it 
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signal if die second transistor is conducting or h^^ will set kto 
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by a transistor in me X statelo VDD. In me third example, if mi nimsfetof had *ero conductance, node 
n 2 would stay charged at 1, while if it has nonzero conductance, me two nodes Wfll Share charge but n^ 
will remain above the positive logic threshold. Therefore itj has a target state of 1, but Terman's program 
sets it to X. In many other cases, however, such as in simulating logic gates with some inputs in the X 



165- 



As a further application of mixed-level simulation, in most MOS designs certain transistor 
configurations arise very often and can be replaced by their functional representations to improve the 
simulator speed. For example, MOSSIM recognizes the transistor configurations corresponding to. 6 
different nMOS logic elements: inverter, Nand, Nor, and each of these logic gates with a single pass 
transistor on its output. MOSSIM performs this optimization only when the configuration comprises an 
entire transistor group. Since transistor groups interact with one another only through fanout 
connections, each group has a clearly defined set of inputs and outputs. Hence the functionality of the 
transistor configuration can be guaranteed to match the functionality of the logic gate. For example, 
group G 2 in Figure 7.1 can only behave as an inverter, and group G 4 can only behave as an inverter with 
pass transistor. These optimizations affect only the speed of the simulation and not its functionality. This 
cautious approach overlooks other possible optimizations but involves no guesswork. With MOSSIM 
these optimizations are performed during the network analysis prior to simulation and entail little 
additional effort because the network must be partitioned into transistor groups anyway. With the new 
algorithm, no such partitioning is required unless the optimizations are to be performed, and hence the 
added cost of applying them becomes much higher. However, for networks which will be simulated over 
long sequences of inputs, the net savings can be significant 

Implementing a mixed-level simulator requires small extensions of the procedures PHASE and 
STEP given earlier. In addition to maintaining the list of perturbed nodes E, we must also maintain a list 
of perturbed function blocks F, i.e. those blocks for which some input has changed since the most recent 
updating. When the procedure PERTURB.TRANSISTORS is called to find which nodes are perturbed 
when a changing node state causes a changing transistor state, we should also call the procedure 
PERTURB_BLOCKS which will add any block to the list F which is perturbed by the changing node 
state. In the procedure STEP, blocks in the list F should be simulated and the nodes at their outputs 
should be set to their new states. Any node which changes state is added to the list A. With this 
implementation function blocks will be simulated much as they are in traditional event-driven logic gate 
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shnulators. 

7.8 Performance of MOSSIM 

Although the algorithms presented in this chapter have not been implemented, their performance 
should be comparable to MOSSIM. Thus we can use the performance of MOSSIM as a measureof me 
speed of switch-level simulation. Furthermore, since MOSSIM ean replace transistor configurations 
corresponding to certain logic gates by a gate-level representation, we can compare the relative speeds of 
switch-level and logic gate simulation. 

MOSSIM is written ia the language, CLU [26J, All times were measured on a DEC 20/60. White 
the program, programming language, and computer system have been designed to provide reasonable 
performance, there is room lor speed improvements in aH three areas. 

Table 7.1 lists the characteristics of six different binary counter circuits, three each of two bask 
designs. Bom designs have me circuit shown in Figure 73 for each bit position. Data is stored 
dynamically, and no initiahzation circuitry has been included. The chain of half adders forms a 
cany-ripple adder. The two designs differ in how the half adders are implemented. The first, called 
CNTR, utilizes four Nand gates and an inverter to implement the sum and ^arry logic in a conventional 
way. The second, called MCNTR utilizes a pre-charged Manchester carry chain as shown in Mead and 
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Logic Gj 


CNTR10-OPT 





70 


CNTRIO-UNOPT 


200 





CNTR16-OPT 





112 


MCNTR10-OPT 


124 


52 


MCNTRIO-UNOPT 


258 





MCNTR16-OPT 


200 


84 



* includes depletion mode transistors 



-167- 



Fig. 7.3. One Bit of Counter Circuit 
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Conway [28, p.150] to implement the carry logic and a selector logic block [28, p.152] to implement the 
sum logic. This design achieves high speed by pre-charging each bit position in the carry chain on each 
clock cycle, so that the carry lines are never driven through load transistors. 

Both designs were tried in 10 and 16 bit versions (e.g. CNTR10 and CNTR16). The suffix "OPT" 
for the entries in Table 7.1 indicates that MOSSIM replaced the transistors for whatever logic gates it 
could find with the logic gate representation. The suffix "UNOPT" indicates that the network was 
simulated entirely at a transistor level. As can be seen, the conventional counter design can be replaced 
entirely by a gate level representation, while in the other design only 50% of the network (measured by 
the number of transistors) could be replaced. Experience has shown that between 50% and 80% of typical 
designs can be represented at a gate level. 

Table 7.II gives performance data for the six circuits. All times are measured in CPU milliseconds 
per clock cycle. Best and worst case times were measured by finding which cases minimized or 
maximized the time, while average times were measured by averaging 1024 clock cycles. The best and 
worst case times could not be measured with complete accuracy for the faster circuits. 
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Program 7.II. Performance Statistics 

CPU milliseconds / clock cycle simulated 

WOISt 

70 

220 

110 

320 
440 
530 



Circuit 


average 


best 


CNTR10OPT 


30 


30 


CNTRKRJNOPT 


95 


170 


CNTR16-OPT 


34 


30 


MCNTRIO-OPT 


300 


260 


MCNTRIO-UNOPT 


400 


360 


MCNTR16-OPT 


490 


450 



First let us consider the data for CNTR. The transistor level simulation requires 3 times longer than 
the gate level simulation. This provides a measure of speed of a simulator which operates at a transistor 
level relative to one which operates at a logic gate level. Furthermore, if the simulator were designed only 
to simulate logic gate circuits, its speed could be improved further. Nonetheless, it shows mat the speed 
of a switch-level simulator can approach that of a logic gate simulator. Observe that the best and worst 
case times differ greatly. This indicates that unless the carry propagates a long way during a clock cycle, 
large portions of the network remain inactive. This is also seen in the 16 bit version where the average 
time is only slightly higher than for the 10 bit version, but the worst case time grows in proportion to the 
network size. 

The data for MCNTR indicate much different performance characteristics for this design than for a 
conventional gate implementation. The circuit requires up to 13 tunes longer (for the 16 bit version) than 
CNTR. This difference is due largely to the greater amount of activity in MCNTR. On every clock cycle, 
all carry lines are pre-charged, and since the sum logic depends on the carry values, the simulator will first 
compute the sum based on the pre-charged value and men on the final carry value. Furthermore, 
transient X states arise due to sneak paths in the push-pull drivers for the carry chain and the sum logic 
causing many more activities to be simulated. Thus, the amount of activity in the network is almost 
independent of the length of the carry propagate. This is indicated by the closeness of the best and worst 
case times and by the fact that the simulation time grows in direct proportion to the network size. Note 
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also that in replacing 50% of the network with logic gate representations, we improve the speed by 25%, 
showing that these optimizations do not provide a major performance gain. 

These measurements indicate that the performance of a switch-level simulator can vary widely 
depending on the nature of the circuit to be simulated. For cimuitsimpleiwated mostly with logic gates, 
the activity is highly localized and much of the design can be simulated at a gate level. For designs using 
more exotic tecbnicpes such as prercharged and pass transistor logic, activity occurs throughout the 
circuit and a larger portion must be simulated at the transistor level to either case, however, the 
simulation time gf ows at most linearly with the network size. 

7.9 Summary 

Unlike previous switch-level simulators which were baaed solely on the intuitions of the designers, 
the algorithms presented here are based on a mathematical theory- TWs provides a framework much Uke 
numerical analysts have in which problems are formulated asa set of equations, and the goal becomes to 
find efficient algorithms to solve them. In our case, the Mutation algorithm relies mainly on a technique 
for solving recurrence equations in the strength algebra. By studying this simplified and more abstract 
problem, one can see more clearly the trade-ofBs between algorithmic complexity and simplicity of 
implementation. Furthermore, the algorithms can be proved correct The characteristics of equations 
which arise in simulating MOS networks such as sparseness and locality of changes can be exploited to 
obtain an algorithm which is particularly efficient for this application. 

Once the basic unit delay algorithm has been devetoped, it can be altered to provide a slightly 
different timing algorithm or extended to improve the efficiency and generality of the simulator. It can 
be seen that switch-level simulation can be combined with logic gate and functional simulation so that the 
best features of each may be utilized. 
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8. Timing Models 

8.1 Introduction 

The unit delay algorithm presented in Chapter 7 simulates die network as moving through a 
sequence of target states. To the external viewer, this provides a timing model in which a transistor 
switches onetime unit alter its gale node changes state, but in which signals propagate along pads of 
conducting transistors instantaneously. For common implementations of inverters, Nand gates, and Nor 
gates, such as those shown in Figure 2.2, this algorithm atee simulates logic gates as having unit delay. In 
this chapter, the modeling of time in MOS circuits will be investigated in terms of both simulating the 
functional behavior of a design and detecting timing errors. Example networks will be given containing 
logic gates. It is assumed that these gates arc implemented m one of trie styles shown in Figure 2.2, all of 
which behave identically ftom a logical point of view. 

As was mentioned in Chapter 2, switch-level networks share many commonalities with relay and 
logic gate networks when viewed as systems computing logical runctions. TT» target state provides the 
base characterization of die logical function computed by a switch-level network. R gives the node states 
created by die network in response to die current stales. Thus it describes the excitation of the network 
much as die excitation of a Boolean gate network [20] is defined as the output values of all logic gates as 
functions of their current input values. Many of the theoretical techniques and algorithms developed for 
logic gate networks (and relay networks) can be adapted to the switch-level model In doing so, however, 
several characteristics of MOS systems must be kept in mind. 

First, die sheer size of MOS LSI systems imposes cons&atets on the practicality and usefulness of 
many techniques. Such tools as Karnaugh maps [20J and flow tables (22, 43) require a complete 
enumeration of all possible network states, which would grow exponential with die number of nodes. 
Such techniques can only be applied to small sections of a design. In fact, any algorithm with nonlinear 
time complexity must be viewed with skepticism for networks in which the elements number in the 
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thousands. Similarly, any tool which requires "hints" from the user such as the location of feedback 
paths, state variables, or delays becomes unwieldy for large networks, unless this information can be 
derived algorithmically from some source such as the layout specification. This problem becomes even 
greater as LSI designs are generated by automated or partially automated systems, because the designer 
may not have the intimate knowledge of the network required to provide such hints. Thus the desire to 
handle large networks and to implement any technique as a computer algorithm changes the standards by 
which techniques are measured. 

As a further point, switch-level networks differ from networks in other logic models in that the X 
level can arise during normal network operation due either to short circuits or improper charge sharing. 
For example, sneak paths between input nodes in different states often arise in pass transistor logic 
circuits such as the push-pull drivers used in output pads [28, p. 165]. Generally these error conditions 
occur transiently and have no effect on the ultimate network behavior. The presence of X states may not 
indicate a badly designed circuit but only a temporary ambiguity in the network operation which must be 
scrutinized to see how far it propagates. This contrasts with logic gate models in which either Boolean 
behavior is assumed at all times or at least that an X value can only arise as the result of other X's. 
Techniques developed for other logic models often require modification to handle the X state. 

Traditional assumptions about sources of delay in digital systems do not apply to MOS circuits, 
either. Most analytic techniques for asynchronous circuits [22, 29, 43], assume that delays occur only in 
logic elements (gate delay) and not in wires (line delay.) Furthermore, they assume that a logic element 
will respond to all inputs simultaneously. This would require in MOS implementations of logic gates, 
such as those shown in Figure 2.2, that all transistors respond to changing input signals at the same rate. 
Such an assumption may not hold, because p-type and n-type devices may have different timing 
characteristics, and even transistors of the same type may not behave identically due to variations in 
geometry, fabrication details, and stray capacitances. A timing model for MOS circuits should allow each 
transistor to have an independent switching time. With this degree of generality, however, we can also 
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model line delay by incorporating it into the switching times of the transistors with their gates attached to 
the wire. Furthermore, transit delay (the time required for a signal to travel through a conducting 
transistor) can usually be modeled by incorporating it into tile switching time of die transistor. Thus, a 
timing model which allows arbitrary delay times for ea* transistor cah inodel all forms of delays in an 
MOS circuit Given that wiring delays are predicted to dominate in future^tSl circuits {36J, the abiHty 
to model wiring delays may play an important rote. Many of the traditional analytic took, however, 
cannot deal with this degree of generality. 

As a final, and somewhat more optimistic note, much of the concern about timing in traditional 
logic design need not concern the MOS designer. Most MOS systems are designed to operate 
synchronously with conservative docking schemes. For estampte, in a properly designed circuit with a 
two-phase, nonoverlapping clocking scheme {28J, no malfunctions due to timing can arise as long as the 
clocks run slowly enough for the circuit to settle duringeaeh time epoch, but fast enough to avoid the loss 
of stored charge. These methodologies have been adopted, in fact, to compensate for the difficulty in 
accurately predicting the exact time behavior of MOS circuits, especially if the design is to operate 
correctly despite variations in fabrication and despite the inability to fine tune circuits once they have 
been fabricated. For such systems, ahnost any timing model would provide sufficient accuracy to test the 
functionality of a design. However, timing still remains an issue for those designs in which relative path 
delays can affect the logical behavior and for ascertaining that a design can operate at a particular 
clocking rate. 



-173 



8.2 Simulation Timing Models 

Much attention has been focused on timing models for logic gate simulators. Some of these 
techniques will be described and their suitability for switch-level simulation will be discussed. 

8.2.1 Unit Delay 

To implement a unit delay model the simulator computes the excitation of the network and then 
sets the node states to these values. This model has been used successfully in both logic gate and 
switch-level [9] simulators. It provides the same level of accuracy as logic designers use when they analyze 
timing by counting logic gate delays. This suffices to model the functionality of a a wide variety of 
designs, because very few circuits rely on a path with fewer logic gates having a longer delay than a path 
with more. A unit delay model, however, may deceive the logic designer who finds that a design can be 
made to simulate correctly if extra delays (e.g. inverter pairs) are inserted along some paths. Often the 
actual circuit cannot be corrected so easily because of factors such as the assymetric rise and fall times of 
ratioed logic gates and the inexact behavior of the circuit during transitions. 

Unit delay simulators can fail to terminate, both in cases where the actual circuit would run 
indefinitely, and in cases where the actual circuit would settle. For example, a simulation of an inverter 
ring would not terminate, such as the one formed when the input to the network shown in Figure 8.1a is 
set to 1, because the circuit would run indefinitely. More importantly, the circuit shown in Figure 8.1b 
has a critical race if the input is changed to 1, but eventually the slight differences in the two Nand gate 
delays would cause the conflict to be resolved (although not with a predictable outcome.) With a unit 
delay simulator, however, both logic gates are simulated as having the exact same delay and the 
oscillations continue indefinitely. As a practical matter, this problem has arisen only a few times out of 
many simulation runs by many users. The conditions leading to these oscillations seldom appear in real 
designs. This nontermination due to perfectly matched delays can create major difficulties, however, if 
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Fig. 8.1. Networks for Which Simulation May Not Terminate 
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the simulator tries to initialize the network nodes to states other than X, because a naive choice could well 
create effects similar to the example in Figure 8.1b. To prevent nonterminating simulations, the 
simulation program can be designed to halt after a maximum number of unit steps, with this limit set 
according to the network size. 

8.2.2 Pseudo Unit Delay Simulation 

A pseudo unit delay simulator proceeds by computing the excitation of a set of nodes resulting from 
an event selected from an event list, where each event indicates a perturbation in the network state. These 
nodes are then set to their excitations and any resulting perturbations are added as events to the end of 
the event list. This process continues until the event list is empty, indicating that the network is in a stable 
state. As long as the event list is maintained as a first-in, first-out queue, this simulator resembles a unit 
delay simulator in that if two nodes are perturbed simultaneously, the effects of both will be simulated 
before any perturbations they cause are simulated. However, activities which are modeled as occurring 
simultaneously (and hence independently) in a unit delay simulator will be ordered, and hence one may 
affect the other. This ordering will depend on the internal details of the simulation program and to the 
user will appear unpredictable. Such a simulator would terminate for the example shown in Figure 8.1b, 
although not in a predictable way, because the simulator would arbitrarily select one logic gate to simulate 
before the other. The network shown in Figure 8.1c, however, has a similar form of critical race, but a 
pseudo unit delay simulator would fail to terminate, because in following a FIFO discipline the simulator 
would alternate between the upper and lower chains, giving the effect of perfectly matched delays. One 
should note, however, that this example is very contrived, whereas networks like that of Figure 8.1b are 
quite common. It is not known whether less pathological examples would fail to terminate with a pseudo 
unit delay simulation in cases where the actual circuit would reach a stable state. Circuits with 
nonterminating behavior such as the inverter ring shown in Figure 8.1a will have nonterminating 
simulations. Thus, a pseudo unit delay simulator partially solves the problem of unbounded oscillations, 
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but at the expense of introducing some degree of unpredfctabflity. The technique has been used 
successfully in a switeh-tevel simulator [5]. 

8.2.3 Zero Delay 

To implement a zero delay model [12] all feedback paths must be broken such that the system 
becomes a combinational network, computing a function from the inputs and current values of die state 
variables to the outputs and new values of the state varabtes. Each pass through the network appears to 
occur instantaneously, and hence the term "zero delay". This modd assumes the circuit is ftee of critical 
races. With this approach an ordering can be placed on the network elements such that each element is 
simulated at most once during a pass, thereby achieving greater efixiency than either a unit delay or 
pseudo unit delay simulator. Simulations of networks such as those shown in Figures 8.1b and&le would 
terminate, assuming that the paths in both cases are broken in «riy one place, although a simulation of 
the network shown in Figure 8.1a would not Such a technique would apply to MOS networks only if the 
feedback loops could be identified automatically. Finding a miaimum set of feedback loops is known to 
be an NP-comptete problem [171 but for most applications a set which fe not minimum would suffice. If 
an algorithm chooses to break the paths in Figures 8.1b and 8.1c in two places, however, a nonterminating 
simulation could result Furthermore, the user would have litueundemand^ of uu simulation timing 
unless informed of the points at which feedback paths are broken. A zero delay model has apparently not 
been tried ma switch-level simulator. 
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8.2.4 Continuous Time 

Many logic gate simulators introduce a continuous time measure by allowing the user to assign 
delay times to the logic elements and using a time-ordered event list scheduler. Some even allow logic 
gates to have different rise and fall times [44}, to model logic gates in ratioed circuits. Such a technique 
has been proposed for switch-level simulation with the network parameters which determine the delay 
computed by a layout analyzer [5]. With MOS circuits, however, delays areaifeetetfby many details 
which the switch-level model ignores, such as both the linear and nonlinear effects of transistors, the 
threshold voltages, the capacitive loadings (which may change during operation), and the exact voltage 
waveforms on the nodes. As one tries to take such details into account, it becomes difficult to achieve a 
consistent -.level of detail without resorting to a fuO scale analog simulation. An inaccurate simulation 
would create more problems titan it would solve, becaHse users tend to place great faith in numerical 
results even if they have no validity. 

Perhaps the most promising approach is to find lower and upper bounds on die circuit delay by 
applying only simplifications of the analog behavior which can be guaranteed either conservative or 
optimistic. The user can then tighten the bounds by increasinglhe level of detail at which the circuit is 
simulated, until it can shown that die required timing condiltoss wdl b» met Recent work by Gtasser 
[18] takes important steps in finding these kinds of simplifications, but much more wor* is required 
before k becomes a practical simulation tooL Tins form of simulation would provide the most reliable 
verification of the ability of a circuit to operate at a partkularcloekkig rate. 
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8.2.5 Summary 

None of the models listed above will satisfy alt uses at all times. However, the unit delay and 
pseudo unit delay models stand out for thek simplicity, understandabihty, a«l insistency with the level 
of detail which the switch-level model is intended to provide. To choose between these two, one must 
compare the value of greater predictability against the value of * greater (but not total) nmriunity to 
nooterminating simulations. 

8.3 Analysis of Timing by Ternary Simulation 

Rather than introducing a continuous time model into a switch-level simulator, one might best 
leave such a level of detail to analog simulators where it can be handled accurately. Instead, the 
switch-level simulator could be used to identify those portions of a system winch warrant closer timing 
analysis by analog simulation or some other technique. Tins would allow the more powerful (but more 
expensive) tools to be applied just where they are needed. For example, the speed of a synchronous 
circuit is often limited by a single critical path, such as the carry chain of an adder. A unit delay simulator 
could generally find tins path by finding which nodes changed state during the last unit step of a 
simulation phase and then working backward. 

A method known as ternary simulation has been developed to detect possible hazards and critical 
races in. logic gate circuits without introducing a continuous tine model {H, 23, 46). This technique uses 
die X state to represent the ambiguity caused by a node in transition from to 1 or 1 to 0. Ternary 
simulation techniques can also be applied to switch-level networks to detect possible sources of timing 
errors. These parts of the circuit can then either be redesigned or analyzed by more detailed methods 
such as circuit simulation. 
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8.3.1 Algorithm 

The algorithm for a ternary simulation of MOS networks can be described in terms of the function 

phase. Suppose a switch-level network is in a stable state y, i.e. y = step x (y), and we v> >sh to simulate the 

effects of changing the input nodes to state x\ The resulting state y' is computed as 

q *— phase{\ U x', y) 
y' 4— phase(x',q). 

That is, first all nodes v. for which X: ^ x': are set to X, and the network is simulated until it settles. Then 
the input nodes are set to their final values x' and the network is again simulated until it settles. For logic 
gate networks implemented with transistor circuits such as those shown in Figure 2.2, this algorithm 
reduces to the one proposed by Brzozowski and Yoeli [11] for logic gate networks. Ternary simulation 
requires only slightly more effort than the unit delay simulator described earlier. Observe, however, that 
the X state will now become the most prevalent state during the simulation, because every time a node 
makes a transition, it must go through the state X. For the method to provide useful information, the X 
state must be simulated accurately and efficiendy. The algorithms presented in Chapter 7 satisfy these 
requirements. 

It is claimed that each node «j will have a new state y'j equal to or 1 if and only if it would have 
this unique state regardless of the magnitude of any delays in the circuit. Any nodes sensitive to network 
delays will have y'j equal to X. This claim has not been proved formally, although it can be motivated 
informally, and weaker statements have been proved. 

The ternary simulation method will first be motivated informally. In setting the input nodes to 
x U x', all inputs which may be changing simultaneously are set to X, indicating that their exact values 
cannot be ascertained. In the first computation of phase, these X's are propagated to all nodes which are 
sensitive to the input node transitions, including into the feedback paths corresponding to unstable state 
variables. In the second computation of phase, the final values of the input nodes are propagated to any 
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nodes which are combinational functions of the inputs and of any stable state variables. If a feedback 
loop has developed containing all X's, however, these nodes and any nodes dependent on them will 
remain in the X state, because the behavior of the actual circuit would depend on the exact delays in the 
feedback paths. The assumption made in our logic model that a set of transistors with the same gate node 
may behave independently when the gate node is in the X state gives the effect of each transistor 
responding to a transition at an independent rate. Thus a node will be set to or 1 if and only if it has 
this unique state regardless of the transistor switching delays in the circuit 

For example, the network shown in Figure 8.2 contains a 2-out-of-3 majority with its output fed 
directly to one input and through an inverter to another. An MOS implementation of the majority gate is 
shown by Seitz in Mead and Conway [37, p.255]. The output of this gate equals 1 or if at least two 
inputs equal 1 or 0, respectively, and equals X otherwise. The reader can also verify that an MOS 
inverter in our model will have output X if its input equals X. Suppose initially that x l = 0, y 1 = 0, 
y 2 = 1, and that we change Xj from to 1. Assuming unbounded transistor switching delays (or 
equivalentiy unbounded line delays), there will be a critical race depending on the relative delays in the 
two feedback paths. A ternary simulation would first set Xj to X, and then compute the new state X for 



Fig. 8.2. Ternary Simulation Example 
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both nodes «j and » 2 - Thus, feedback loops have developed containing only X's. When xj is then set to 
1, both nj and B2 will remain at X, indicating a critical race. This example demonstrates that a MuHer 
C-element (consisting of a majority gate with its output fed back to One of its inputs> implemented in 
MOS may malfunction if its feedback path contains as excessive delay. Iri feet, Dennis and Patil have 
shown [14] that a C-element cannot be constructed which is guaranteed to function despite arbitrafy line 
delays. 

8.3.2 Theoretical Results 

Brzozowski and Yoeli [11] have proved that an algorithm similar to the one shown here will set a 
node in a logic gate network to 8 or 1 ^^jfithasftbuniijaes^ereg»fdtes«of fcgfe &tfe delays. Some 
networks which would function properly assuming only unbounded togic gate delays, such as the 
example^ Figure 8.2, however, may have nodesset to X. lteseau«Wisc«njecatr€,Jbut do not prove, 
that their algorithm will seta node to or 1 If and only (fit has this unique^state regardless of both gate 
and Out delays. Then* proof reMes primarily on the moaotoaicity of tfw *«cifeition function for the partial 
ordering < X and 1 < X. Their proof also assuni» that aklgfcnetworkfiioimaBy contains Only nodes 
with Boolean logic states, but this requirement can be relaxed. The excitation function for the 
switch-tevel model is step % , and we have already seen by theoiem 46 that it is monotonia Thus, 
Brzozowski and Yoeli's proof can be applied to show mat for y' resulting from the ternary simulation, 
each element y' l will equal or 1 only if it has dm unk|i»«afe»regHadlB» of the transator swindling 
delays, with the restriction that transistors with the Same fate node must have Htfr* same delay. 
Furthermore, a proof of their conjecture could be applied to show that each element y'jwal equal or 1 
if and only if it has this unique state for arbitrary transistor switching delays. 
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Brzozowski and Yoeli's proof also shows that the ternary algorithm will always terminate after a 
linear number of simulation steps (in the number of nodes.) For example, the simulations of all three 
networks shown in Figure S.i will terminate, j(but with the nodes in fte X stated) Even though such 
circuits as the inverter ring shown in Figure &la will am indefinitely, their ternary simulations wiH 
terminate. Thus, the worst case performance of ternary simulation isijetter man that of any of the timing 
models described previously. In practice, this potential gain is offset by the larger number of node 
transitions and the greater effort required to simulate the effects of the X state. 

8.33 Application 

A ternary simulation would provide useful results fi^ synchronous systam, because a well designed 
synchronous system should contain no critical races under any delay nxidel.lBsing errors can only arse 
as a result of insufficient clock intervals, and these should be checked by critical path analysis. Far 
asynchronous systems, on die other hand, Dennis and Pattt £14 have shown that no nontrivial sequential 
circuits can be built which will function correctly despite arbitrary line delays. For example, me Mufler 
C-element forms the cornerstone of most speed independent circuit designs [29], and we have already 
seen how it can fail with ternary simulation. A ternary simulation of nu^ asynchronous circuits would 
overwhelm me user with X's and provide Btfle information about the design. Unger (43] describes how 
delay elements could be introduced into a ternary simulation. By mseftun^ deky elements mto carefully 
selected feedback paths* circuits can be made to simulate in a reasonable manner. This style of 
asynchronous design, however, has limited application ia MOS LSI, because it requires too much 
fine-tuning of die circuit. 
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8.4 Toward a Simulator for Self-Timed Systems 

In some design disciplines, timing constraints are assumed to hold within small regions of the 
circuit, but as long as these constraints are satisfied the system will function correctly regardless of the 
delays in other parts of the circuit. For example, many speed-independent designs require only that the 
C-elements operate properly. Seitz has proposed a class of systems called self-timed systems [36] in which 
particular timing constraints can only be assumed to hold within equipolenlial regions, while arbitrary 
delays may be incurred by any communication between regions. This methodology allows much 
flexibility in the design style, because a subsystem contained within a single equipotential region can be 
implemented with a synchronous or asynchronous circuit, or recursively as a self-timed system. 
Self-timed systems may also have overlapping equipotential regions, but we will not consider this 
possibility. Ternary simulation would provide a useful tool for testing self-timed systems if it were 
applied only to those portions which are to function correctly despite arbitrary network delays. The 
simulator should only model the functionality of those portions for which particular timing conditions are 
assumed to hold, thereby providing the appropriate stimulus to the portions simulated by the ternary 

algorithm. 

Such a simulator would have many advantages over detecting timing errors with a continuous time 
model simulator. With a continuous time simulator, numerous cases must be simulated with slight 
changes in operating conditions and network parameters. Even a poorly designed circuit may simulate 
acceptably due to a chance combination of input conditions, element parameters, and simulation model. 
With ternary simulation all possible forms of uncertainty due to timing are condensed into a single state 
and hence a single simulation run analyzes many cases simultaneously, always finding the worst case 
behavior. With a self-timed system simulator, the user could isolate small regions of the design for which 
particular time behavior is assumed to hold. These regions could be tested extensively be a very accurate 
circuit simulation. Then the simulator tests the hypothesis that as long as these regions function correctly, 
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the entire design will be insensitive to other circuit delays. 

A system could be simulated in this manner if the regions of the circuit for which ternary simulation 
is to be applied were separated from those regions to be modeled with a pseudo unit delay simulation by 
special network elements as will be described informally. These elements serve to convert between the 
two modes of simulation. Regions may only be connected through fanout connections, i.e. from a node 
in one region to the gate of a transistor in another. 

Ramp elements convert the 0-1 and 1-0 transitions resulting from a pseudo unit delay simulation to 
0-X-1 and 1-X-O transitions for the ternary simulation. When the input to a ramp element changes, the 
output is first set to X, initiating a transition event, an indication of which is placed on a special event list 
Any time a perturbation resulting from a transition event causes an X to propagate to a node, this change 
also becomes a transition event. Events are selected from the normal event list only when the transition 
event list is empty. The simulation is continued until the network settles (i.e. both event lists are 
emptied), at which time the X's will have propagated as far as possible. Then the simulator sets the 
outputs of those ramp elements with output equal to X to their input values, and the effects of these 
events are simulated until the network settles. It can be seen that the combination of ramp elements and 
this scheduling algorithm leads to a ternary simulation of those regions with ramp elements at their 
inputs. 

Trigger Elements convert the 0-X-1 and 1-X-O transitions resulting from a ternary simulation into 
0-1 and 1-0 transitions. That is, if a transition event causes the input to a trigger element to be set to X, 
the output is held in its previous state. When the input is set to a new state by a normal event, the output 
is set to this state. In this way, the X states arising from the ternary simulation will be blocked from 
entering regions for which the simulator is only to model the functionality. 
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While this simulation technique has not yet been implemented, it shows great promise as a tool for 
locating possible timing errors in self-timed systems. 

8.5 Summary 

The relatively simple unit delay and pseudo unit delay taning models have proved adequate for 
modeling the functional behavior of both synchronous and self-timed systems. These models provide 
little information about possible timing errors, but for eireuitft designed with conservative docking 
schemes, such information is not required. Acontinaous time model simulator would prove most useful 
for verifying that a circuit can sustain a particular doddafcHtt; Such a simulator should try to place 
lower and upper bounds on circuit performance by applying only simpMcations which can be 
guaranteed conservative or optimistic. It should allow the user tettgWen^hese bounds by increasing the 
level of detail at which the circuit is modeled. Ternary and self timed system simulators provide a 
powerful tool for detecting critical races. Unlike analog simulators, winch can only verify a design for a 
particular set of circuit delays, these simulators can verify a design for aft possible sets of circuit delays, 
except fw small regies m which partkmlaraom^ 
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9. Conclusions 
9.1 Final Thoughts 

The value of the switch-level model has already been proved by extensive experience with 
switch-level simulators. These simulators have shown peat versatility and accuracy in modeling the 
logical behavior of systems constructed using many different architectural and circuit design techniques. 
Furthermore, they utilize only information about the actual structure of the design such as can be 
obtained from the mask specifications. By operating at a logical level, switch-level simulators achieve 
performances approaching those of logic gate simulators. 

This thesis has demonstrated that tfie switch-level model can also be developed into a mathematical 
description of the behavior of MOS logic circuits. This allows a rigorous derivation of equations to 
express the operation of a network from which simulation algorithms can be derived and proved correct 
These new algorithms improve on previous algorithms which were based en the programmer's intuition. 

The degree of generality allowed by the switch-level model does not come without its costs. As 
compared to relay and logic gate models, our mathematics model seems much more complex and 
intractable. A long and difficult process was required to derive the basic simulation technique. This 
difficulty stems partly from the novelty of the model. The switch-level model cannot be described 
adequately by such well-developed concepts as Boolean and linear algebras, and hence we were forced to 
develop a new abstraction and justify it in terms of an electrical model. Furthermore, the switch-level 
model supports a much richer variety of circuit and architectural design techniques than do traditional 
logic models. The Mead and Conway approach to custom LSI design differs from other approaches such 
as polycetls [33] and gate arrays {19] in that it allows the designer to select from die entire variety of 
different design techniques and even tailor the individual devices to provide many trade-offs between 
speed, power, density, and ease of design. To provide sufficient expressive power for this level of 
generality, a logic model must be more complex than logic gate models, but as has been shown tilts 
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g^eralityneed«)tconMttirou^iarfAoc«ttensioM. 

Unlike commercial LSI design* however,. Mead and Conway encourage me use of structured 
design methodologies and simplified design technique* This pennies a heavy reliance on computerized 
design and analysis tools to replace the many hours of highly skilled labor used to produce commercial 
LSI designs, The tools can be designed along the hues of mis memodotogy, thereby achieving better 
performance and encouraging me user to produce well structured designs. With the emphasis shifting 
away from techniques appropriate for bumansitoB^oseapprop^is ^ scomputer iraptementlitiott, 
techniques should be judged primarily on their ability to be implemented as computerized tools. In this 
regard switch-level simulators compare favorably with logic gate simulators and greatly outperform 
analog simulators. 

9.2 Suggestions for Further Research 

Thus tar, the switch-level model has only been implemented hi die form of logic simulators with 
sunpte timing models. White these applications have demons&ated the vataeoff me switch-level model, 
they represent only the begianmgs of an entire class of tools for MOS design. The switch-level model 
helps bridge the gap between the views of as LSI design as an electrical shcuit or as a piece of artwork 
and the view as a system which computes a logical ftmctton. 

9.2.1 Simulation 

Simulation will always play an important rote in the LSI design process. Humans have a 
remarkabte abaay to synthesttc designs, often applying §reat cleverness m selecting from a seemingly 
endless range of possible approaches. However, many hours of tedious work and much self-discipline are 
required to generate a totaHy correct design by hand. Computers, in contrast; lack the kind of cleverness 
and intuition which goes into novel and original designs but will willingly perform tasks which humans 
find impossibly dull Thus a natural complement is formed with computers verifying the designs and 



188 



ideas produced by humans. Simulation provides just that form of verification. It directly tests die 
functionality of a design in much less time and at Mich less cost than the production of prototypes. 
Furthermore, it can often provide more information ten can reasonably be measured from an actual 
implementation of a design. 

Future logic simulators could provide more rigorous checks of the correctness of a design. Unlike 
current simulators, which try to provide as much generality as possible, the simulator could be tailored 
toward a particular design methodology to check whether Hw design adheres to this methodology. Many 
of these checks can be added with little additional complexity. 

A slight modification of the algorithm given in Chapter 7 would yield a simulator which provides a 
more rigorous test of CMOS designs. In CMOS, node voltages correspond to valid logic Wets only if 
they very nearly equal 0.0 or V dd . Even the threshold voltage drop resulting from a signal with state 1 
passing through an n-type transistor or a signal with state passing through a p-type transistor is 
considered excessive. Thus, a 1 signal passing through a conducting n-type transistor should become an 
X and similarly foraO uiiough a p-type, unk^ the coinplemeittar^ 

simulator can achieve this effect by modeling a turned-on n-type transistor as baling conductance equal 
to its strength for signals of state and having an unreliable conductance for signal* of state 1, and 
conversely for p-type transistors. The techniques developed for monetae transistors in the X state can be 
extended to describe the behavior of die network for this model with only slightly increased effort 

The self-timed system simulator described in Chapter 8 would test whether a circuit really fulfills 
the requirements of self-timed system design. It does so in a more rigorous and reliable way than would 
an analog simulator and with much less computational co^ Ihealgoritrtm for this simulalor must still be 
developed in detail, however, and the design of a suitable user interface will present many challenges. A 
more fbnrwl understanding of the meoiyoehuxith^ Brzozowski and 

Yoeli's conjecture that a ternary simulation tests whether a design wiH function properly for arbitrary line 
delays remains an open problem. The proposed combination of ternary and pseudo unit delay 
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techniques has also not been studied formally. This style of ^ntilatten could greatiy assist die design of 
self-timed systems. 

Simulators could ate© be designed to provide new kinds of tafbrmation. For example, whereas 
existing simulators only tell die user the state of the network, a wore sophisticated program could also tell 
how it got there. This would greatly aid fte user m debugging a design. 

Simulators could also provide more information about circuit timing. Even with methodologies 
such as conservative clocking schemes and tools such astenM^stenibtion, for some applications theuser 
must know the time behavior of a circuit For example, tf a system must function at a dock rate 
determined by other chips or by the particular application, the designer must verify that this clock rate 
can be sustained. Hopefully, a fu» scale analog simuktiem ear* be avoided by identifying the criticat patii 
and simulating only this part with '-a simplified model It is beieved tint a considerable amount of detail 
must be added to tire switch-level model before it can model circuit timing with the same generality and 
accuracy with which it models the functionality. While the "Order of magnitude" electrical model may 
describe the logical behavior adequately, it provides limited information about circuit speed. Ideally, any 
simplifications of the electrical behavior should be guaranteed cotiservative or optimistic so that (he 
simulator provides bounds oft the performance. The simulator should be able to apply different tevelsof 
detail so that the bounds can be tightened until they are acceptable for die application. Unlike existing 
hybrid simulators in which increased accuracy is achieved only by going to a totally different (and often 
incompatible) model, the simulator should allow a smooth and reliable transition between the levels of 

detail. 

A switch-level simulator eould also provide information about the effects of circuit faults. Unlike 
the traditional mcmods of fault analysis, which assume mat a logic gate wffl fail by becoming "stuck-aT 
some level, the switch-level model can describe the failure of individual transistors or connections. For 
example, a faulty transistor can be modeled as a transistor in me X state, Le. an unknown and unreliable 
conductance. Similarly, a faulty wire can be modeled as two nodes connected by a strong transistor in the 
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X state. This application of the X state has not been studied in detail 

Finally, we must consider what form of simulation will work for VLSI design. The current 
implementations of switch-level simulators model the entire design as a network of transistors or at best 
combine some of these transistors into simple logic gates. While the efficiency of this technique has 
proved adequate for modeling LSI systems containing up to 10,000 transistors, it will dearly fail for VLSI 
systems containing 100,000 or perhaps 1 million transistors. As with all other design tools, the size and 
complexity of VLSI systems forces a rethinking about logic simulation. Systems of this size must be 
designed hierarchically, where each level of the design involves combining subsystems designed at die 
next lower level. A simulator must utilize this hierarchy to achieve the required performance and to assist 
the user in maintaining the hierarchy. Switch-level simulation can provide valuable assistance at the 
lower levels of this hierarchy. It can be used to verify that a transistor circuit implements a particular 
logical function either by automatically comparing k to a functional specification or simply by allowing 
the user to test the design. Once the functionality has been verified, die simulator can replace the 
switch-level representation with a functional representation for simulation at the next level of die 
hierarchy. In some cases the user may also wish to simulate some portions of the network at a functional 
level and other portions at a switch level. As was seen earlier, a simulator can be designed to combine 
simulations at these different levels. Thus switch-level simulation wtt serve an important role in VLSI 



9.2.2 Other Design Took 

Most existing LSI design tools treat a design as a piece of artwork. For example, most layout 
systems and design rule checkers view the design as a collection of geometric primitives with no 
understanding of the actual or intended functionality. As a consequence, they provide only limited 
assistance to the designer. The switch-level model can help introduce a greater understanding of die 
functionality of a design. 
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Current design rule checkers, for example, check only the geometric features in the layout 
Experience indicates that they catch only a limited class of errors and also report many extraneous errors. 
For example, if the designer omits a contact cut, no error is reported. If a wire contains a gap, it will be 
reported only if the gap is smaller than the interwire spacing tolerance. A more useful checker would 
perform both the geometric analysis and the extraction of the switch-level logic network. It would 
produce two results: the switch-level network it believes was intended, and the design rules violated by 
the layout with respect to the network. By simulating the network and checking the list of violations, a 
designer would have a much greater chance of detecting most errors. Furthermore, since the checker 
would have a better understanding of the electrical connectivity in the layout, it could greatly reduce the 
number of extraneous errors reported. Experience with C. Baker's layout extraction program [4, 5] has 
already demonstrated that many errors in the artwork can be detected by deriving the switch-level 
network and applying several simple checks. 

The switch-level logic model could also be incorporated into a design "analyzer" which detects such 
features in an MOS logic design as feedback paths, sites of dynamic storage, potential race conditions and 
short circuits. Many of the theoretical techniques which have been developed for the Boolean gate model 
could be profitably adapted to the switch-level model. 

The ultimate goal of computer-aided design is to automatically synthesize an LSI or VLSI design 
from a high level functional description of the system. Current automated design systems such as Bristle 
Blocks [24] still require the user to design large portions of the layout by hand and have only an 
incomplete understanding of the circuit's functionality. These tools would benefit to some degree by 
adopting the uniform representation of a logic design provided by the switch- level logic model. 
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9.2.3 Theoretical Developments 

The mathematical developments in this thesis provide a firm foundation for switch-level simulation, 
but much more theoretical work is required before the full potential of the switch-level model can be 
realized. Several areas have already been mentioned in connection with their possible applications. For 
example, the exact relation between circuit timing models and both ternary and self-timed system 
simulation techniques remains to be established. Similarly, research has only begun on methods for 
introducing more detailed timing information in ways which can be guaranteed conservative or 
optimistic. Finally, developing computer algorithms to identify feedback paths, sites of dynamic storage, 
and other features of a design will require a much better theoretical understanding of the structure of 
MOS logic networks. Many of these analysis problems appear to be NP-complete, but heuristic 
techniques can probably be devised which will work efficiently for most real designs. 

The simulation model developed in this thesis provides only an operational model of a logic 
network. That is, given a particular initial state, it describes the state in which the network will ultimately 
stabilize. A more powerful model would describe the function computed by the logic network. For 
example, Shannon showed how the Boolean function computed by a combinational relay network can be 
derived from the structure of the network. The Boolean function for a logic gate network can generally 
be derived by inspection. Work on denotational semantics of programming languages is also directed 
toward deriving the function computed by a computer program. Much more work will be required to 
develop a corresponding functional model of switch-level networks. At present, the conversions between 
the logic states, logic signals, and signal strengths tend to obscure the function computed by the network. 
Furthermore, expressing the ternary function computed by a network (i.e. including the state X) involves 
considerably more effort than expressing the Boolean function. 
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A functional model of switch*Ievel networks would have many applications. For example, it could 
be used to verify that a network implements a particular function. Unlike programming language models 
in which this task is in general uncomputable, for logic networks we can always construct truth tables by 
simulation, and hence any design can be verified in at most exponential time. Furthermore, various 
heuristic techniques could reduce the complexity considerably for most problems. Therefore, automatic 
logic design verification may be practical. A functional model could also help make the transition from a 
switch-level simulation to a functional level simulation; That is, ifme 5 f8ricttsh computed by a section of 
a design could be derived automatically, the simulator could replace the swifeh-level representation by a 
functional representation. This seems especially practical if it is applied to small sections of the design, 
such as to the transistor groups used in MOSSIM. Since most designs- contain many instances of a 
particular transistor configuration, only a small number of different amfiguratiotts would need to be 
analyzed. Thus, developing functional models for switch-level networks poses a great challenge but 
should yield many applications. 
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Appendix I • Multi-Port Networks 

This appendix contains derivations of several equations and proofs of some properties regarding 
passive linear resistor networks. 

1.1 Introduction 

Suppose a time-invariant network N contains only voltage sources and passive, linear resistors in 
which each voltage source has its negative terminal connected to the reference node GND and is set to a 
voltage 0.0<v<V < y. Assume furthermore that each node is connected by some path to a voltage source 
and that voltage sources are never connected in parallel, and hence aU node voltages are well-defined. 
The corresponding zero- slate network Nq is defined as the network formed when all voltage sources in N 
are set to 0.0, ie. they are short-circuited. The network Ng contains only passive, linear resistors. 

A port in the network consists of two terminals with any node servingas the positive terminal, 1 and 
the reference node GND serving as the negative terminal. We will be interested in networks with at most 
three ports, labeled 1 through 3. The indices i and j are used to denote any of these ports. We require 
only information about N and Nq which can be observed at the ports. Otherwise the networks can be 
considered "black boxes". 



1. Nodes which are the positive terminals of voltage sources in N correspond to input nodes in the 
switch-level network, while other nodes correspond to normal nodes. Either kind of node can be a 
positive port terminal 
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1.2 Port Parameters 



The open-circuit impedance parameters are defined in terms of the voltages measured at the ports of 
the zero-state network N Q when a current source is conned across one of the ports. That is, if current 
source I is connected across port i, and a voltage v '■ is measured atport j, then Zjj is defined as 



v 



z ji ■"" 1 * 



These parameters form a matrix Z which in our case has dimension 3X3/ Most presentations of 
multi-port networks [15] describe the analogous 2x2 matrix for two-port networks. The diagonal 
elements of Z are the "selFimpBdance" terms, Le. the tmpedarice measufed across each port of N Q when 
afl other ports are left open^ircuhed. The off-diagorial elements are trie "cross-impedance" terms 
indicating the degree to which the positive tenninalS of the two pom are connected; That is, Zy equals 
0.0 if there is no path in N Q between the positive terminate of ports i< and j, 1 or if one of these two 
terminals is connected directly GND mrdogh a zero-state voltage source. At the other extreme, Zy equals 
z if every path in N from the positive terminal of port i to GND passes through the positive terminal of 

portj. 

For a passive resistor network, the matrix Z obeys several important properties. First all elements 
are nonnegative and finite. Second, since the network is reciprocal, the matrix is symmetric: 

zy = z«. foralliandj. 

This follows directly from the Reciprocity Theorem R5]. Third, for z u >&D, when a voltage source with 
a positive voltage v is connected directly across port i of Ng, while all other ports are left open-circuited, it 
will produce the same effect as a current source with current v/zy. Therefore, the voltage at any port j 



1. As the term "path" is used here, a path cannot contain the reference node GND or any node 
connected directly to a voltage source as an intermediate node. 
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under these conditions will be 

Tliis voltage cannot be less than 0.0 or greater than v, and therefore 

0.0 < zjj < zy, for an i and j. (All) 

This result also holds when zy = 0.0, because thiscase occurs only when the positive terminal of port i in 
Nq is connected directly to OND by a zero-state voltage source* and hence za . ~ OJD. 

The open-circuit voltage parameters of N are defined as the voltages vj, *%, and v^ measured at 
ports 1, 2, and 3 of N with all three ports open-circuited. Parameters of this nature have notbeen found 
in any other presentations on multi-port networks. In general, these voltages wiQ not be completely 
independent of one another when the network contains paths between fte positive port terminals. To 
derive some of their relations, suppose that a voltage source of voltage v^Js connected across port i of N. 
Then no voltages in N wiH be changed. Now if afl voltage sources ia N are set to 0.0, die voltage at any 
port j will be given by 

V J Zii v i- 

when zjj ^ 0.0. Furthermore, since this voltage was obtained by changing the settings on some voltage 
sources from nonnegative values to 0.0, it must be less than or equal to die original open-circuit voltage of 
port j in N, and therefore 



Zy 

0.0 < -*vr < V:, for alii such that 2=: ^ 0.0 and all j. (A1.2) 

Zj; J «• 



The Thevenin equivalent of the network N at port i contains a voltage source set to the open-circuit 
port voltage and a resistor of conductance equal to the net conductance across the port with all sources set 
to 0.0. As long as the other ports are left open-circuited, the Thevenin equivalent of N at port i is 
described by the parameters: 
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v thev = *i 
Sthev = L0/2 ii- 



Suppose the resistor conductances are given by rational functions of p with degree greater than 0.0. 
Then the open-circuit parameters wiftateo be given by rational ftactkaas of p with degree less ttiaa 0. 
Furthermore* since O-O^ijCpi^zgCp) fifr all vaUiesofp: 

ffegCzy) < JegCzjj) < foraUiandj. 

Similarly, the open-circuit voltage parameters will be given by .rational functions of p, and since 
O.O^v^p^V^ for all values of p 

<fe«Cvj) < foraUL 

1.3 The Effect of a Variable Resistor 

A network containing a single variable resistor can be described by -.& fixed network N with a 
variable resistor connected across pom 1 and 2, as shown in Figure LI. jSince the positive terminal of 
port 3 can be any node in the network, the port voltage can represent the voltage on any node in the 
network. Suppose the resistor has conductance b> and die resulting port voltages equal w' lf v' 2 , and v' 3 . 



Fig. 1.1. Multi-port Model of Network with Variable Resistor 
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This resistor will have the effect of injecting a current I into port 1 of N and removing the current I from 
port 2, where 

I = IKv'j - V'j). 

By superposition, the voltage at any port equals the sum oftte voltage created by the sources of N, the 
voltage created by the current injected into port I, and tve voltag* created by the<Mrrent removed from 
port2: 

V'j = Vj + IZy - IZjj (A1J) 

y'i = v i + H^i ~ V^u - r^ 

For ports 1 and 2, this gives two equations in the two unknowns v.\ and v' 2 : 

v'j = vj + hW x - y' 2 Xz u - z 12 ) 
v' 2 = v 2 + iKv'i - v'jXz^ - zjjV 

This set of equations can be sobed to give 

» » v l ~ v 2 



»i - v*, = 



1"'2 - I.0 + h(z 11 -2z 12+ z 22 ) 
and this result can be substituted into equation A1.3 to give 



V l V l + 1.0 + IKZ U - 2Z 12 + Zjg) iA1 ^ 



In particular, the voltage at port 3 win be given by 



3 ~ v 3 + 1.0 + h(z 11 -2z 12 +z 22 )- (A1 ^ 



To gain some insight into this equation, observe that for smaB values of h, v'j - v' 2 pwvj - v^ 
and therefore 

v' 3 ~ v 3 + My i - v 2 Xz 13 - Z23* (A1.6) 

In other words, the denominator in equation A1.5 approximately equals 1.0. As h becomes very large, die 



m 



increased conductance between the two nodes will be offset by the decreased voltage difference, and 
hence each node voltage in the network approaches an asymptote: 

Urn v = + ^- v 2^ 1 3- 2 23> , (A1 ,7) 

h-»oo * * z ll _2z 12 +z 22 

Thus Vj is approximately a linear function of h fe^JmaftVaMAllF h and levels off to a constant value for 

large h. 

The general form of equation AL5 is 

v<tO = v(W + 1.0 +\ h - < 34 > 

Furttiermore, since z^2 <2ji and 2j2 <^22 

b = z n - 2z 12 + z 22 = (z u - z u ) + (Z22 - z 12 ) > 0.0. (A1.8) 

Equations 3.4 and A1.8 are sufficient to prove Lemma 3.1. 

We can derive further properties of equation A1.5. When a positive current J. js injected into port 1 
and removed from port 2 of Nq, the positive terminal of port 1 must have die maximum node voltage in 
the network and the positive terminal of port 2 must have die minimum. Therefore 

< z 12 ~ z 22> ^ < z 13 ~ z 23> ^ ( Z U ~ z 12>- 

The second inequality holds when I is negative as well. From this one can see that 

I z 13 — z 23 1 ^ Z U ~ 2z 12 + z 22- (AL9) 

If the resistor conductances are given by rational functions of p, then the voltages will also be a 
rational function of p and will have a general form 

V(p,h(p)) * V(p,0.0) + w^f^y 

Furthermore, since the impedance parameters all have degrees less than Q, degHb) < 0, and since equation 
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A1.9 holds for all values of p: 

de^Y - v 2 Xzi3 - z 23 )) < de& u - z£j) < de&z n ~ 2z 12 + z 22^' 
Therefore 

These results are sufficient to prove Lemma 3.3. 

As a historical note, although equation A1.S would seem to have otter applications in electrical 
network theory, no such result has been found in the literature. Some presentations on sensitivity analysis 
(such as [15, p.678]) derive the approximate equation A1.6 for small values of h, but die more general 
result is not given. The technique of viewing the variable resistor as a current source of unknown current 
and then later solving simultaneous equations to find this current is generally credited to Kron [25] which 
he used in a method called "diakoptks". More recently, mis technique has been applied to sotviflg 
systenisofsparseequatk)nsbyameftodcalkd"tearm;j"I7,32l. 

1.4 The Effect of Connecting Two Ports 

We can determine the effects of connecting together ports 1 and 2 in network N by letting the 
conductance h in the previous development approach infinity. We will assume that z^>0.0 and 
Z22 > 0.0, i.e. there is no voltage source connected directly across ports 1 or 2 of N. In the limit die 
voltages on the two ports will approach a single voltage v where 



y _ Ibn y 

h-»<» 



, = v , < v l - v 2 Xz U - z 21 > 
1 1 + -*u- 2 «U + ?» " 



By rearranging terms we get 

= < z 2? - z 1 2 > v l + ^ 1 1- ^ 2, 
E ll _2z 12 + 2 22 



v = •»■ • 



1. This restriction can be assumed, because the combination rule is never applied directly to die logic 
signals which describe input nodes. 
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If we define kj and k 2 as 



k l = Z 12 /2 11 
k 2 = *12^ z 22 



then 



v thev ~ gjCMJ-k^ + fl^^i) v ' 

where §j and g 2 are the Thevenin conductances at ports I aad2, te. die lecipWcals of z^j and z 22 , 
respectively, and from our assumptions both values must be positive and finite. The following properties 
of die factors kj and k 2 can be derived from their definitions and from equations A1.1 and AL2: 

92 k l ~ 9l k 2 
0.0 < k r < LQ 

0.0 < k 2 < 1.0 

k r v i ^ v 2 

k 2' v 2 ^ v l- 

Equation 4.4 and the above properties of die factors kj and k 2 are suffice* to prove the validity of die 
combination rule for cyclic connections, as is done in Chapter 4. 
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Appendix II ■ Proofs of Results in Chapter 6 

This appendix contains proofs of Theorem 6M, Corollary 6.4.1, and Theorems 6.S and 6J~ 
n.l Passive Networks 

The method of conditioned relaxations relies on me fact that logical conductances act as passive 
elements. That is, when a signal is coupled through a logical conductance, the resulting signal has 
strength less than or equal to the original signal strength. Furthermore, signal combination always 
ignores the weaker signal. This allows us to kill any signal on a node wim strength less than the steady 
state signal strength, knowing that no possible action of the network could amplify this signal into one 
critical to the formation of some steady state signal. 

The passiveness of logical conductances has important implications on the recurrence equations 
describing a logical conductance network. Before me desired theorems can be proved, some general 
properties of "passive" recurrence equations must be derived. For a value s £ 5 let s^ denote the vector 
bbck (&,[&]), where [s] denotes a vector wkh each element equal to s. Similarly, for a function 
f:fi -» 3°, let/ s denote the functk>n/ s (a) = block (/(*), fs J* 

A function/:/ 1 — » J° is said to be passive if for any s and any a, / s (a) — f s (t^}. In other words, if 
h = /(a) then any element of b greater than or equal to s can depend only on elements of a greater man 
or equal to s. No element of a less than s will be amplified into a value greater than or equal to s. A 
function of more than one argument is passive if it is passive for each argument. The functions t. i, \ 
Mock and constant functions are all passive, as is any composition of passive functions. The successor 
function sue, on the other hand, is not, where sue is defined as suc(yj = y ff and for any a < y ff suc(&) 
equals the least element of T greater than a. 
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The following lemma shows the relation between the minimum solutions of the equations a = /(a) 
and a = / s (a). 



Lemma All. 

For a monotonic and passive function f:f* -♦ i°, if a"" is the minimum solution of the equation 
a = /(a) then a"^ is the minimum solution of the equation a = / s (a). 



Proof of Lemma A2.1: 

We wffl show by induction on k that 



This clearly holds for k = 0, so assume it holds for k-L 
The minimum solution of a = / s {a) is given by 

k iToc/>> = k^TooM = *V ■ 

Observe that this property may not hold if the function/ is not passive. For example, the minimum 
solution of the equation a = sucii) equals y„ but for any s > 0, the minimum solution of the equation 
a = sucJa) equals 0, even though the function sue is monotonic. 

We can now prove a lemma mat expresses in a very distilled form that the method of conditioned 
relaxations will not lead to a "runt" solution with some signals weaker than the steady state signals. 



-20*- 



LemmaA2.2. 

For a monotonic and passive function f-.f 1 -*^, if a"* 1 B the minimum solution of the equation 
a = /(a) then a""" is also the minimum solution of me equation a = /"(a) where 

/(a) = bbck{m,**% 



Proof of Lemma A22: 

Let a"*" equal the minimum solution of the equation a = /(a). Clearly a"* 1 is also a solution of this 
equation and therefore a B,i,1 > a"*". We will prove that the two vectors are ia fact equal by induction on 
decreasing strength values. By Lemma A2.1, a mh, s and a"*" s are the minimum solutions of the equations 
a = / s (a), and a = /^a), respectively for any value of s. The function Nock obeys the following 
identities 



Woe* v fcd) = c „ 
blocked) = NockiCyd^^. 



That is, unless the second argument is greater man me first, Nock behaves as an identity function of its 
first argument Therefore 

This implies mat jT* t ^Jf*^ . N» suppose a^a^ = V^w^- We will show that 

We know that »"*" aj ^ s ) < """"s. aod therefore the left hand argument of Nock in the above equation 
must be greater than or equal to the right hand argument This means mat this application of Woe* wffl 
behave as an identity function of its first argument giving 
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which shows that a"* 1 ^ = a""^. The set 5 is totally ordered and finite. Hence by induction on decreasing 
strength values a" to ' = a - ". I 

This result also may not hold for a function / which is not passive. For example, the equation 
a = sudd) has a minimum solution y ff but the equation a = block (suc(a), y J has a minimum solution 0. 

The result of Lemma A2.2 can be expressed in a form closer to what is required to prove Theorem 
6.4. 



Lemma A2J. 

For a matrix G E S° x n , and a vector b € JL n , define the functions /j and/ as: 

jr^u) = b1ock(fr\ fffii r) 

/ (d) = «ocA(UU 4„jG-i,xX 



where 



r a G*"l*fc 



If a"* 1 and d"* 1 are the minimum solutions of ttie equations v = f^a) and d = /j)(d), respectively, then 

r = u^r** 



Proof of Lemma A23: 

Define the function g as 



g(a) = Woc*(l*ITG-a, r). 



Lemma A2.2, combined with Corollary 6.2.1 shows that r must be the minimum solution of the equation 
a = g (a). We will show by induction on k mat 



g k (0) = /^WT/o^O). 
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Clearly this holds for k = 0. Now suppose that it holds for k-L 

g k (0) = block $ A BTG'g^OXr) = Woe* (I All T G-C/^'^O) T f Q l 'Ho)), r) 
g\0) = block(rniG'ff\Q\r)Xbbck(LiJ1G%*-h9\i) = f^O) T f H<». 

Taking die limit of this equation gives the desired result. 

II.2 Theorem H4 

Lemma A2.3 takes care of the most difficult part of the proof of Theorem 6.4. 



Theorem 6.4. 

For a matrix G E f 1 x n , and a vector A £ JL n , define G as G = xG. The unique minimum solution of 
the equation a = /(«), where 

/(«) = * V Go a (6.19) 

is given by 

rf* = +rf- V -d"*\ 

where h"*" and d™*" are the minimum solutions of the equations u - /j(u) and d = / (d), respectively, 
and the functions /j andyjj are defined as: 

f t (u) = block (r AT f G-*r) (6 JO) 

/ (d) = Atoc*(lAJTG-4rX (6.21) 



asd 

r = G*-IAi. (6.22) 
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Proof of Theorem 6.4: 

Define the vector * as * = +u"*' V -d* to . We will first show that A satisfies die recurrence relation 
A = /(A). Observe that Tiki = W** and LAJ = il"*". Lemma A2.3 then shows that 

r = u" to T^* = ■*! = liltG'fll = lAVGoAl. 

Substituting this into the definition of/j gives 

f^rm) = bbck(Tbi t G»m, ia v <7oa«) = r* v <7oA1. 

Similarly, 

/ (LAJ) = LAV GoAJ. 

Since I"A1 = a"", m = /^rJUand similarly LAJ » /^LAJ). This shows that 

A = +r*1 V -LAJ = +r* V (7oAl V -LA V GoAJ = A V Co*. 

Therefore A satisfies the recurrence relation A = /(AX which implies that A > rf" 1 ". By the monotonieity 

ofl i,n,andLJ: 

IAB > H*l (A2.10) 

u- 1 " = TA1 > f/*1 (A2.H) 

d"" = LAJ > LfM. (A2.12) 

We can now show that A must equal rf*" by factoring the equation rf Bfa = /(,£*"). First we can find the 
strength of flf as 

»V rfB « = iiVGo^I = IA1 T G'l^*'!. 

Therefore il a"* 1 II must be greater than or equal to the minimum solution of this recurrence relation, Le. 
H ^ ste B > r = HAN. Combining this result with the inequality of equation A2.10 shows that 
N rf""" N = BAH. Now we can see that 

rv ,a i = rA v Goj**-\ = block (r m t G«Trf Bi H,ii<r ,l,, i) 



"i;?;'- 



• 208' 



Fa"" 1 ! = Uock{Tbl T G • ^1, r) = /^ra*-!). 

Therefore ra"*"l > n"*', whk* when combined wiiv the inequality of equation A2.ll gives 
IV"! = u"*. A similar derivation gives Mr*" J = f". We can now complete |he proof with 

a-" = +rrf*i v -Lrt = +B-* y -f*. I 

II.3 Corollary 6.4.1 



Corollary 6.4.1. 

For a matrix G € 5° x n , and a vector * 6 Mi if*risuieiiied«s tim *«£ thin the miriimum solutionof 
the equation a = /(a) where 



is given by 

where 

and 



/(a) = AVGoa. 



/(a) * M/(fla), lX (625) 



r = G*«l»l. 



Proof of Corollary 6.4J: 

Expanding the equation for/ gives 



/(a) = *///(* V Go a, r) = +Wbc*(r* V Go ml, r) V -4foc*(l# V Go ml, r). 

In the proof of Theorem 6.4 it is shown that r = B d* I and therefore tor any a < a"* 

r = la""! = IM T C'lV*! > 1*1 t G«lal = II V Goal. 

Fo^anya<a ,,, ■ 

block (Vi V Co«l, r) * ttockiUmckim^ H*U\ 1* V Goal), r) 
Moc*(rA V Goal, r) = block (Til T G • fal, r) = ^(faTX 
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where /j is defined in equation 6.20, and a similar result holds for the part with respect to the function 
/q defined in equation 6.21. Then for « <^" 

/(«) = +/i(r«i) V -/ (L«J). (A24*) 

We will show by induction on k that 

/*m- = +/ t k (o) v viNw < A114) 

This clearly holds for k = 0, and assuming it holds for k implk*4fe* 



Fotu"*' and d""" defined in the statement of Theertm 6Afft® < «""* and/ k (0) < d"*'. Therefore, by 
the monotonicity of +, -, and V: 

f\o> = *f l k m V -/<>) < +■- V -d- = ^", 

and equation A2.13 applies when a = f\®k This allows to ejqj«a4/ k+ 1 as: 
f k+ \0) = WHO = +/ 1 (r/ k (01) V -/qO/*(0J> 

which proves the induction assertion. Combining the result of Theorem 6.4 with equation A2.14 gives 

^ = + k Too'^ V "k ^oo^ = k ^co ^** V "'<>>. 
and therefore 

^ = v ^Z*** ■ 

k-»oo 



-210- 



IL4 Theorem 65 



Theorem 6.5. 

The target state y of a switch-level network is given by 

where n** and d°* are the minimum solutions of die equations i = gj(n) and d = g Q (i\ respectively, 
and the functions gj and gg are defined as 

g t (u) = Afart^r^-rxit r/i-tfl^'M)" (6^9) 

= block &?.• LxJ t kjU T G— • 4, rX (6J0) 



and 

r = cK'OF* •■-*!? ■*!* (631) 



Proof of Theorem 6.5: 

First, the idea behind the optimization method can be derived fbnnally% algebraic mariipulatlon. 

The last step above follows tram the identity shownin equation 5A It relies on the fact that unless the 
state of a signal equals X, either its part or its 1 part must equal 0. Furthermore observe that the 
function of b whose value is <+b> is monotonk, and therefore 

«+a>U<H>> = <+(atb)> 
and that a similar identity holds with -. Hence 

' = ^tGMii 1 "* 4 ^ u ^c/w*^ van 
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Since | denotes the pointwise maximum operation, this equation shows that the target state for any node 
nj can be computed by rinding the maximum values of rvj(G, E)1 and Lvj(G, E)J for all G G { G } and 
E€{E}. 

Define ^"(G, E) and (P^G, E) as the minimum solutions of the equation u = /j(u) and d = / (d), 

respectively, where 

/j(u) = Woc*(E-r*-iT ryi t g«u, g*«(e-«xiitnj'>))- 

/„(d) = block (E • LxJ T Lj^J T G • d, G*-(E • II x II T II y »))• 



Theorem 6.4 shows that 



u^G.E) = rKG,E)l 
d^G,E) = LKG,E)J, 



Therefore, if we can prove that 



H * = t u" to (G,E) (A2.16) 

{G},{E} 

d°* = T d^G.E), (A2.17) 

{G},{E} 



then comparing equation 6.28 to equation A2.15, it can be seen that the theorem will be proved. 

Only the proof of equation A2.16 will be shown. Equation A2.17 follows identically. For any 
GG{G},EG{E}anduG/ 1 

E«r*i r ryi r g»u < Er"-rxi t ryi r <?-•■ 

G*-(E - II jell T Hj'H)) > G min *-(E min • H jc II T II J' «))- 

Therefore, since block is monotonic in its first argument and antimonotonic in its second, f± < g^. 
Applying Theorem 6.3 gives u mln (G, E) < u 01 * for any G and E in the allowed range. Therefore 

T u min (G,E) < u 01 *. 
{G},{E} 
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Tocomplete the proof, we need only find some setting of Gand E which^es*"*^, E) « a 8 * 
Define toe matrices G' and E' as follows: 



g'.. = rs^ij. u o, * i = 0oru* j = 

,J I sT^ yT t >o and^ >m- :j 



e a = ( e y. r T = « 



We will show that u^G', E*) = a * From these definitions, wecan see that 

E'Txl = E-Txl, (A2.18) 

and can show that for any u < a** 

Woe*(G'*M) = #£**(&*"•«, r). (A1I9) 

To see this, observe that for a < a *, 

block (G' • a, r) < block (G- • a, r) < Woe* (G« • a*, r) < block (a *, r). 
Therefore, for u * =0, 

< Woc*([G'«Bl i ,r i ) < block (vP v r$ = Q, 
and for «"*>(), 

WaciOG'.aJj.ri) = N <« k (jJ >Q <*'%l «jX fj) 
W^aG'-aJj.ri) = ^^^^(g-yiupirp 
W«r*aC';al i ,rp = block (fG— • ajj. r$. 

Equations A2.18 and A2.19 imply mat a* is a solution to the equation 

« = block®. rvir iyif £•■.!)"' (A2.20) 

and furthermore any solution less than or equal to a " must also satisfy the equation a = ^(uX and 
therefore a * is the minimum solution of equation A2.20. Lemma A2.2 then shows that a* must also be 
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the minimum solution of the equation 

u * Uock(E' • TjcI t r>1 T G' • n, r T a "*). 

Lets = G'*»(E' • II x i T B jl)). It isctaimed that s = r T uf. Cteatiy 

r = G^^E^-lxltijil)) < C*-(E'«*j:ITIj>D) = «• 

We know that sis toe miaimum solution of the equation a ^$a>*diere 

Ma) = E'«HjciT»^iTG'«a. 

This shows that s > n"* because the function in the recurrence equation A2.20 is less than or equal to h. 
Therefore s>rtir*. Next we wttt snow that r T w 8 " = A(rT o^, which wffl prove matfT if*> Vand 
when combined with the previous inequality, shows that r t u**. = s. First observe that 

G'*(r frf*) = G*t t G'»i°* 

This follows because if q > u *, then u * = 0, which implies tiiat the fth column of G' will equal the Ah 
column of G"*\ Similarly 

e'.i x i = e^-ii x « t E'-rxi. 

We also know mat r satisfies the recurrence relation 

r = E^-lx! T ■>■ t <£***• 

Combining these facts with equation A2.20 gives 

„<*t r = E' • Txl T Tjrl T G' • U** t ■» • 

u^T r = E' Txl T E""" ' ■ *« T T^ T '*>■ T G' • u* T G^ • r 

u^fr = ES»xM l/l T G'-fr^Ttf 
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Putting these results together, we have shown that u* is the minimum solution of the equation 

b = Woc*^'.rjciTrjTTG'.«,<3'*Kr'i*BT*ji)X 

and therefore u " 1 = rKG',E')1. TMs completes ourpwbfltat 

A similar result can be proved for d*i whWia^*^ 
IL5 Theorem^ 



Theorem 6.7.. 

If y = target (x, y, z\ then y = target (x, y , *). 



Proof of Theorem 6.7.: 

Define the vector y as the set ofsfenafc formed by the norttial nodes in state y;Le. 

<y> = y 

■/■ = cap. 

Observe that I / 1 = I y I = can, and therefore for r defined in equation 6J1 

r = G- h *-(E-. lxlTI> | ) = G**^!-.!,!!,];^ 

Compare die functions 

Let u«* and S* be the minimum solutions of the equatioiisn = ^(oX and. == ^(uX respectively. We 
win show mat these two vectors are equal. First, r j?,l = u^capj for all i Furthermore. 
u*j > block (f^l, rp, and tyl < cap; for all i. which implies that r~y{\ > block (l^l, ifl. Therefore. 
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since block is monotonic in its first argument 

block (Ty\i) > block (block (Vyl,!), r) = block (r/l.r), 

which shows that g x > g 1? and by Theorem 6.3, u " 1 > u 01 ". We will now show that these two - lues are 
in fact equal by showing if* = gjOi " 1 ). Since u " 1 >Tyl> block (Vyl, r), 

g 1 (u opt ) = u " 1 = u opt TWocJt(rJl,r) = i^O t War* (r>l,r), 

and furthermore g 1 (u opt ) > g^u "') = u m > block (r/1, r) which gives 

Thus u "* = u " 1 , and by a similar argument we can show that d 01 " = H 01 *, where d* "' is the minimum 
solution of the equation 

d = block (E max *LxJ | I-3U | G max • d, r). 

Theorem 6.5 shows that target (x, y , z) is given by 

target (x,y,z) = <+u opt > U <-3 0pt > = <+u opt > U <-i opt > = target (x, y, z). I 
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